Abstract

MRP2 encoded by ABCC2 gene is involved in the secretion of numerous drugs and endogenous substrates. Patients with Dubin-Johnson syndrome due to mutation in ABCC2 gene have elevated urinary coproporphyrin ratio (UCP I/(I + III)). Here we investigated whether this ratio could serve as a biomarker of MRP2 function. Phenotype-genotype relationships were studied in 74 healthy subjects by measuring individual UCP I/(I + III) ratio obtained on 24-hour urine and by analyzing five common SNPs in ABCC2 gene. The UCP I/(I + III) ratio varied from 14.7% to 46.0% in our population. Subjects with 3972TT genotype had a higher ratio ( ) than those carrying the C allele. This higher UCP I/(I + III) ratio was correlated with a higher level of isomer I excretion. This study provides a proof of concept that UCP I/(I + III) ratio can be used as a biomarker of MRP2 function in clinical studies as it provides quantitative information about the in vivo activity of MRP2 in a given patient.

1. Introduction

Excessive accumulation of drugs and consequent drug-induced toxicity may arise from a defect in the ability to secrete drugs into the bile or urine. Various membrane transporters, including those of the ATP-binding cassette (ABC) family, are involved in these processes in the kidneys and liver. One such transporter, multidrug resistance protein 2 (MRP2), encoded by the ABCC2 (ABC subfamily C2) gene, plays a key role in the secretion of numerous drugs and drug metabolites [1]. Studies in Mrp2-deficient rodents have confirmed the relevance of this transporter in drug pharmacokinetics, with exposure levels reaching 140% to 615% of those in normal animals, depending on the drug considered [25].

MRP2 is also involved in the elimination of various endogenous compounds. Uroporphyrins and coproporphyrins, natural intermediates of heme biosynthesis, are the main porphyrins in excess in the plasma. Coproporphyrins exist in two isomeric forms (isomer I and isomer III), both physiologically excreted by the liver and the kidney. Although no direct lines of in vitro evidence indicated a role for MRP2 in coproporphyrins transport, animal studies recently confirmed the major role of Mrp2 in the elimination of these compounds from the plasma [6]. The authors used a bile fistula rat model to study the biliary and urinary excretion of coproporphyrin isomers in wild-type and Mrp2-deficient rats receiving a continuous infusion of equimolar amounts of the coproporphyrin isomers. Mrp2-deficient rats excreted larger amounts of coproporphyrins in urine, with very low levels in bile compared to wild-type rats. The proportion of each isomer in urine also differed between the two groups of animals. Isomer III predominated in the urine of wild-type animals, whereas the proportion of isomer I was higher in the urine of Mrp2-deficient rats, accounting for 70% of the coproporphyrins present. These results confirmed that MRP2 is involved in the excretion of coproporphyrins and that the impairment of MRP2 activity leads to a switch in excretion pattern, from excretion in bile to excretion in urine.

Dubin-Johnson syndrome (DJS) is an autosomal recessive disorder in which the production or function of MRP2 protein is impaired [7, 8]. Missense, nonsense, or splice-site mutations in the ABCC2 gene are responsible for this syndrome. It is characterized by conjugated hyperbilirubinemia and a deposition of a melanin-like pigment in hepatocytes, with liver function otherwise normal. DJS patients also present a modified pattern of coproporphyrin isomers excretion in urine, characterized by a high urinary coproporphyrin I-to-coproporphyrin I + III ratio (UCP I/(I + III) ratio) [9, 10]. Subjects carrying a homozygous mutation in the ABCC2 gene have a UCP I/(I + III) ratio greater than 80%, whereas in normal subjects it is only about 30%. Subjects with heterozygous mutations have moderately high UCP I/(I + III) ratios (about 60%) [9]. Based on results from studies on the distribution of transporters in liver and kidney [11] and results of Moriondo et al. [6], Frank et al. [9], and Kaplowitz et al. [12], we propose a possible physiopathologic pathway for coproporphyrin elimination (Figure 1).

Apart from mutations, several common single nucleotide polymorphisms (SNPs) in the ABCC2 gene have also been identified [13, 14]. The phenotypic consequences of these SNPs have been evaluated indirectly, by studying the pharmacokinetics of several MRP2 substrates. These studies showed an impact of particular SNPs (−24C/T, 1249G/A, 3563T/A, 3972C/T, and 4544G/A) on methotrexate [15], mycophenolic acid [16], irinotecan [17], pravastatin [18], and doxorubicin [19] elimination. Identifying patients with low levels of MRP2 activity before the administration of a drug that might act as a substrate for MRP2 could therefore potentially lower the risk of overexposure. However, no phenotypic marker of MRP2 activity has yet been identified. We hypothesized that the pattern of coproporphyrin elimination in urine would constitute a potential biomarker for MRP2 activity. We thus analyzed phenotype-genotype relationships, by studying the excretion of coproporphyrin isomers in urine and by genotyping several ABCC2 SNPs in a group of healthy volunteers.

2. Methods

2.1. Subjects and Study Procedures

Eighty healthy human volunteers were recruited for the COVOL study (Clinicaltrial.gov number NCT00746044). Caucasian male and female subjects were eligible for participation, provided they met the following criteria: at least 18 years of age, no clinical abnormality (investigator’s clinical judgment), no porphyria (according to the urine analysis results), and normal renal and liver functions. Subjects had no chronic treatment other than oral contraception in female subjects (comprising ethynylestradiol, 20 to 30 μg) and took no medication in the two weeks preceding the study. The COVOL protocol was approved by the ethics committee of Tours University Hospital, France. All subjects gave written informed consent for genetic testing and urine analysis. A blood sample was drawn from each subject for SNP analysis. For coproporphyrin determinations, urinary samples were obtained according to the following timetable: day 1, 5 mL of the morning micturition sample was dispensed into aliquots, and all subsequent urine over the next 24 hours was collected. We also collected 24-hour urine samples on two other occasions, within a maximum of two months. Sodium carbonate was immediately added to samples as an oxidizing agent, and samples were then kept frozen until analysis. Procedures for urine handling and storage were established before this study [20]. Urine volumes were precisely measured.

2.2. HPLC Analysis of Coproporphyrins and UCP I/(I + III) Ratio Determination

Urinary coproporphyrin levels were determined by HPLC (high-performance liquid chromatography), using a method described in more detail elsewhere [20]. Briefly, analyses were carried out with a chromatographic system composed of a pump and an Ultimate 3000 autosampler (Dionex, Sunnyvale, Calif, USA). The 474 fluorescence detector (Waters, Saint-Quentin-en-Yvelines, France) was set at 365 nm for excitation and 624 nm for emission. Chromatographic separation was achieved on a reverse-phase C18 Symmetry (Waters, Saint-Quentin-en-Yvelines, France) column (4.8 × 250 mm) with a precolumn. All chemicals were HPLC grade. The mobile phase was a gradient of acetonitrile (Carlo-Erba, Val de Reuil, France) and 0.015 M sodium acetate (Merck, Darmstadt, Germany) buffer, adjusted to pH 4 with glacial acetic acid (Carlo-Erba, Val de Reuil, France).

Urine samples (400 μL) were acidified by adding 20 μL 10 N hydrochloric acid (Carlo-Erba, Val de Reuil, France), vortexing for 20 seconds and centrifuging the mixture for 10 minutes at 18000 × g. We then injected 50 μL of the supernatant into the chromatography system. The limits of quantification were 7 nmol/L for isomer I and 10 nmol/L for isomer III. The UCP I/(I + III) ratio was determined with a precision of about 99%.

We quantified both coproporphyrin isomers in urine from standard curves obtained with commercially available calibrators (Recipe, Munich, Germany). The UCP I/(I + III) ratio was obtained by dividing the peak height of isomer I by the sum of peak heights of isomers I and III.

2.3. DNA Analysis

Genomic DNA was extracted from peripheral blood leukocytes according to the kit manufacturer’s protocol (FlexiGene DNA kit, QIAGEN, Hilden, Germany). Five SNPs within the ABCC2 gene were selected for genotyping on the basis of previous studies evaluating their impact on MPR2 activity [13, 19, 2123]. The five SNPs (−24C/T (rs717620), 1249G/A (rs2273697), 3563T/A (rs17216324), 3972C/T (rs3740066), and 4544G/A (rs8187710)) were genotyped by direct sequencing. PCR was carried out in a total volume of 50 μL containing 100 ng of total DNA, 10 pmol of each primer (Table 1) [24], 0.2 mmol/L dNTPs, 1.5 mmol/L MgCl2, and 1.5 units of Taq polymerase (Invitrogen, Carlsbad, Calif, USA). The amplification protocol was as follows: initial denaturation at 94°C for 5 minutes, 35 cycles (1 minute at 94°C, 30 s at 60°C (rs17216324, rs2273697, and rs8187710), 56°C (rs717620), or 65°C (rs3740066), and 1 minute at 72°C), and a final extension at 72°C for 10 minutes. Then, PCR products were purified with a DNA gel extraction kit (Millipore, Billerica, Mass, USA) and sequenced with the Big Dye Terminator (version 3.1) cycle sequencing kit (Applied Biosystems, Carlsbad, Calif, USA). Sequencing analysis was carried out on an Applied 3130XL automated DNA capillary sequencer (Applied Biosystems, Carlsbad, Calif, USA).

2.4. Statistical Analysis

The normality of the UCP I/(I + III) ratios’s distribution was tested using the Shapiro-Wilk test. A Mann-Whitney test was performed to compare UCP I/(I + III) ratio between men and women. UCP I/(I + III) ratios were grouped into 4 categories based on distribution according to quartiles. The lower quartile (Q1) included the data of the 25% of patients with the lowest ratio; quartiles Q2 and Q3 extended 25% below and above the median ratio, respectively; the upper quartile (Q4) included the 25% of patients with the highest ratio. Student’s test was used to compare the concentrations of the two isomers in patients with the highest ratio (Q4) versus others.

We used a Bland and Altman plot and its intraclass correlation coefficient (ICC) to compare the UCP I/(I + III) ratios obtained for the morning sample and the 24-hour urine samples. The same test was used to evaluate the stability of the UCP I/(I + III) ratio (obtained for the first, second, and third 24-hour urine samples) over time in a given subject (intraindividual variability on three separate occasions).

All SNPs were tested for Hardy-Weinberg equilibrium. A search for a relationship between UCP I/(I + III) ratio (obtained for the first 24-hour urine samples) and each SNP was first performed using a Kruskal-Wallis rank sum test and the post hoc Siegel-Castellan test. Then, a Mann-Whitney test was used to compare UCP I/(I + III) ratio as a function of each ABCC2 polymorphism (wild-type allele carriers versus homozygous variant for each SNP). Statistical analysis was performed with version 2.8.0, an open source program developed by the Foundation for Statistical Computing (Vienna, Austria). We used THESIAS software (Paris, France) [25] to construct ABCC2 haplotypes and to analyze associations between haplotypes and UCP I/(I + III) ratio.

3. Results

3.1. Distribution of UCP I/(I + III) Ratio in Healthy Volunteers

Eighty subjects were recruited (38 men and 42 women). Mean age was 25.2 (±6.5) years, mean weight was 63.1 (±9.9) kg, and mean height was 169.0 (±15.4) cm. Five subjects were excluded from the analysis because we were unable to obtain the totality of samples (urine or DNA) and one was excluded because he took a nonsteroidal anti-inflammatory drug (NSAID) during the study.

UCP I/(I + III) ratio varied from 14.7% to 46.0% in our group (as determined from the first 24-hour urine sample) with a median at 28.0% (Q1: 25.1%; Q4: 33.6%; Figure 2). The distribution of UCP I/(I + III) in our population was not Gaussian ( ). No difference in UCP I/(I + III) ratio was found between men (median (interquartile range)) (30.0 (7.8)%) and women (28.3 (9.1)%; ). In parallel, we did not observe any correlation ( , ) between the UCP I/(I + III) ratio and the age of the subjects. Subjects with a high UCP I/(I + III) ratio (greater than 33.6%) had higher isomer I concentrations in urine (mean ± sd; 25.6 ± 12.0 nmol/L) than did the other subjects (16.5 ± 8.6 nmol/L) ( ), whereas no difference in the concentration of isomer III was observed between these two groups (41.3 ± 24.5 nmol/L and 48.5 ± 29.7 nmol/L, resp.; ).

3.2. Stability of UCP I/(I + III) Ratio over Time in Healthy Volunteers

Good concordance (ICC = 0.67 with a confidence interval of [0.42; 0.96]) between the UCP I/(I + III) ratios measured for the morning urine sample and those obtained for the 24-hour urine samples was observed in our population (Figure 3). The exclusion of the one subject with discordant UCP I/(I + III) ratios (black circle) increased the ICC to 0.90, with a confidence interval of [0.77; 0.97]. It is therefore possible to use morning urine samples in place of 24-h urine samples.

Comparisons of the UCP I/(I + III) ratios obtained for the same subject on three separate occasions (24-h urine samples) showed that the UCP I/(I + III) ratio was stable over time, with a median coefficient of variation (CV) of 7.8%, except in three subjects for whom CV values of 21 to 27% were obtained. The ICC values of these subjects were in the normal range, at 0.80 [0.64; 0.87]. Thus, UCP I/(I + III) ratio remains stable over time in a given subject.

3.3. Genotype-Phenotype Correlation in Healthy Subjects

The results of the genotypic analysis for the 74 volunteers are summarized in Table 2. The “minor” allele frequencies (MAFs) of the −24C/T, 1249G/A, 3563T/A, 3972C/T, and 4544G/A polymorphisms of ABCC2 were 0.22, 0.15, 0.08, 0.34, and 0.08. All polymorphisms were in Hardy-Weinberg equilibrium. We found a linkage disequilibrium between −24C/T and 3972C/T ( , ) and complete linkage disequilibrium between 3563T/A, and 4544G/A ( , ). Haplotype analysis was then performed for four SNPs (−24C/T, 1249G/A, 3563T/A, and 3972C/T). Five haplotypes were observed: CGTC ( [frequency] = 0.42), CGTT ( ), CGAC ( ), CATC ( ), and TGTT ( ).

In a preliminary trend analysis, we observed a nonsignificant difference in UCP I/(I + III) ratio according to 3972C/T genotype ( ). The same was true for the −24C/T SNP ( ). When we grouped together CC and CT subjects, we found that subjects who were homozygous TT for the 3972C/T SNP had higher UCP I/(I + III) ratios ( ) than allele C carriers (Figure 4(a)). The same was true ( ) for −24TT subjects (Figure 4(b)). We found no relationship with any of the other three ABCC2 SNPs studied. Similarly, no association was found between haplotypes and UCP I/(I + III) ratio.

4. Discussion

This is the first study exploring the role of MRP2 polymorphisms in coproporphyrin elimination in humans, with the ultimate aim of evaluating UCP I/(I + III) ratio as a biomarker of MRP2 function in vivo.

Measurements of the UCP I/(I + III) ratio in 74 healthy subjects effectively demonstrated wide interindividual variability, with a median (range) of 28.0 (14.7–46.0)%, consistent with values previously reported for the general population [9]. By contrast, UCP I/(I + III) ratio remained stable over three different measurements made in individual subjects taking no drug that might interfere with MRP2 at any time during the study period, indicating that this ratio is indeed an individual characteristic. Note that some female subjects took an oral contraception ( ) comprising ethinylestradiol (EE), known to decrease bile flow and to decrease expression and activity of Mrp2 in rats when administered subcutaneously at 5 mg/kg/day [28]. Consistent with the very low EE dosage used, we found no difference in UCP I/(I + III) ratio between women who took EE and those who did not (Mann-Whitney test, ).

Genotype analysis was conducted with the five SNPs of the ABCC2 gene previously shown to be involved in the pharmacokinetics or toxicity of MRP2 substrates. This genetic strategy was specifically selected to study the phenotypic consequences of these frequent common SNPs, and not of rare mutations leading to DJS (incidence = 0.05%). The minor allele frequencies (MAFs) for the ABCC2 SNPs in our subjects were consistent with those reported in the HAPMAP project on a Caucasian population (residents of Utah descended from immigrants from northern and western Europe) [29] and none of these SNPs deviated from Hardy-Weinberg equilibrium. As 4544G/A was in total linkage disequilibrium with 3563T/A ( , ), it was not included in the haplotype analysis.

Our phenotype-genotype analysis revealed that homozygous 3972T and −24T subjects had higher UCP I/(I + III) ratios than wild-type allele carriers. We found no statistical relationship with the other polymorphisms of ABCC2.

These results are concordant with experimental results. The 3972C/T polymorphism leads to a decreased MRP2 protein expression associated with an impaired transport activity [30]. As the 3972C/T polymorphism is synonymous, its influence on MRP2 activity is probably mediated by the variation of another SNP in very strong linkage disequilibrium [16]. We effectively found a linkage disequilibrium between this SNP and the −24C/T polymorphism, although not complete ( , ), explaining why the trend for an association between −24C/T polymorphism and UCP I/(I + III) was not significant contrary to 3972C/T. Another hypothesis is that 3972C/T SNP may affect protein conformation and substrate specificity [31], modulate protein production by altering mRNA secondary structure and stability [32, 33], or lower translation efficacy [30]. The −24C/T polymorphism is located in the 5′-UTR, a region that appears essential for ABCC2 transcription in liver cells [34]. Similarly, the −24T allele was associated with reduced ABCC2 mRNA levels, [21] and reduced MRP2 protein expression albeit without significant difference in transport activity [30]. This SNP could affect ABCC2 mRNA stability and/or translational activity [30].

Laechelt et al. [30] showed in cells transfected with −24T or 3972T ABCC2 cDNA that 3 haplotypes, H9 (−24C/1249G/3972T), H10 (−24T/1249G/3972C), and H12 (−24T/1249G/3972T), were associated with altered MRP2 protein expression (33%, 73%, and 44%, resp., as compared to wild-type haplotypes) and altered efflux rates (significant for 3972T-containing haplotypes H9 and H12). This confirms that both polymorphisms influence the function of ABCC2 transporter, the effect of 3972C/T being predominant. We did not find any relationship between haplotypes in our population and the UCP I/(I + III) ratio. Our results show a SNP-dependent coproporphyrin elimination and confirmatory studies in a larger population are required to analyze haplotypes and to increase the study power.

Several authors also assessed the functional consequences of ABCC2 SNPs for the pharmacokinetics of drugs acting as substrates of MRP2. In their study on pediatric cancer patients, Rau et al. showed that children carrying the −24T allele had lower levels of MTX clearance than other children, this difference being most marked (tripling of the area under the curve AUC) in female patients [15]. Naesens et al. showed an increase of the mycophenolic acid AUC (17% at day 42 after transplantation, ) in renal allograft recipients carriers of the −24T allele and a similar relationship in patients carriers of the 3972T allele [16]. Recently, Laechelt et al. found that subjects bearing the 3972T allele with or without the −24T presented an increased AUCoral of talinolol. This indicates that screening for ABCC2 polymorphisms may help to better understand the interindividual differences in the pharmacokinetics of some drugs.

Beyond the genetic factors, the activity of MRP2 transporter can also be modulated by exogenous factors, and particularly drug interactions, by phenotypical modification of its activity [1]. Personal data, not yet published, show that the basal value of UCP I/(I + III) ratio in a given subject does indeed increase following the administration of an MRP2 inhibitor, suggesting that UCP I/(I + III) ratio reflects MRP2 activity at a given time in a given patient. In subjects of the present study taking no drug, comparable UCP I/(I + III) ratios were obtained for morning micturition and 24-hour urine. So, based on a simple urine sample, UCP I/(I + III) ratio is a practicable test for future prospective studies evaluating its predictive value for the pharmacokinetics of MRP2 substrates.

The precise mechanism controlling the secretion of coproporphyrin isomers into urine has yet to be fully described. Evidence for the involvement of MRP2 in this process is provided by the higher UCP I/(I + III) ratio observed in DJS patients and Mrp2-deficient rats and is supported by the results of our study showing that a moderate impairment of MRP2 function due to polymorphisms increases this ratio, although to a lesser extent than in subjects completely lacking the protein. Frank et al. suggested that the higher UCP I/(I + III) ratio in DJS patients resulted from the excretion of larger amounts of isomer I in urine [9], but these authors did not quantify each product separately. We showed, with 24-hour urine samples, that a high UCP I/(I + III) ratio (subjects from Q4) effectively resulted from an increase in isomer I excretion in urine. We are currently conducting transport studies on MDCKII (Madin Darby Canine Kidney II) cells to verify if one or both coproporphyrin isomers are substrates of MRP2. Increased excretion of isomer I in urine, associated with an inverted proportion of the two isomers in case of MRP2 defect (increase of the UCP I/(I + III) ratio), suggests that coproporphyrin elimination may involve at least a second transporter in the kidney. The proposed physiopathologic pathway in Figure 2, accounting for the inverted proportions of isomers in Dubin-Johnson patients (ratio ~80%) compared to healthy subjects (ratio ~30%), is consistent with the intermediate ratio (about 40%) observed in subjects with lower MRP2 activity, such as subjects with 24TT or 3972TT genotypes. In this case, only part of the coproporphyrins are excreted into the bile via MRP2, again with a higher proportion of isomer I (biliary ratio ~70%) than isomer III, due to the higher affinity of MRP2 for isomer I. The remaining coproporphyrins are excreted in urine via MRP2 or other transporters (such as BCRP, MRP4). Knowing that there is less of isomer I than of isomer III which arrives at the kidney, no competition occurs and the UCP I/(I + III) ratio is about 35 to 45%.

In turn, UCP I/(I + III) ratio provides quantitative information about the in vivo activity of MRP2. This would account for the variation of UCP I/(I + III) ratio with ABCC2 genotype, being high in patients with DJS mutations, moderately high in subjects with relevant polymorphisms, and normal in wild-type subjects.

In our opinion, this study provides a proof of concept that UCP I/(I + III) ratio can be used as a biomarker of MRP2 function in clinical studies. Indeed, it is stable over time in a given subject, provided that no environmental factor (e.g., drug interaction) modifies MRP2 function. It depends on genetic characteristics, and the administration of an MRP2 inhibitor increases the ratio (personal data). Thus, UCP I/(I + III) ratio may reflect the functionality of MRP2 in a given patient, at a given time, taking into account individual genotype and concurrent environmental factors. We now need to carry out a prospective study to demonstrate that the UCP I/(I + III) ratio measured at the time of drug intake is effectively correlated with the pharmacokinetics of the drug. We are currently carrying out such a study with methotrexate, a well-known substrate of MRP2.

5. Conclusion

MRP2 is involved in the elimination of several drugs. No biomarker of its activity is available so far. Previous studies have shown that the urinary coproporphyrin ratio UCP I/(I + III) is elevated in patients with Dubin-Johnson syndrome, as a consequence of mutations in the ABCC2 gene. In healthy volunteers, the ratio was highly variable between subjects and tended to be elevated in subjects homozygous for two distinct SNPs, 3972C/T and −24C/T. This study provides a proof of concept that UCP I/(I + III) ratio can be used as a biomarker of MRP2 function in clinical studies as it provides quantitative information about the in vivo activity of MRP2 in a given patient. Based on these results, we propose the use of UCP I/(I + III) ratio to evaluate MRP2 activity in future pharmacokinetic studies.

Abbreviations

MRP1:Multidrug resistance-associated protein 1
MRP2:Multidrug resistance-associated protein 2
MRP3:Multidrug resistance-associated protein 3
MRP4:Multidrug resistance-associated protein 4
BCRP:Breast cancer resistance protein.

Conflict of Interest

The authors declare no conflict of interest.

Acknowledgments

The authors thank all the volunteers who contributed to the study. They would also like to thank the staff of the Clinical Investigation Center of Tours University Hospital (research nurses, clinical research fellows, project manager, and biostatisticians) and the Centre de Ressources Biologiques de Touraine (CRBT). Additional thanks go to Professor Philippe Lechat for his scientific contribution. This study was funded by grants from the Ligue Contre le Cancer, the Entreprises du Médicament (LEEM), and the CANCEN association (Cancérologie Centre).