TITLE: 

  ABSTRACT:
Despite its established inter-individual variability, sildenafil has been the subject of only a few pharmacogenetic investigations, with limited data regarding the genetic modulators of its pharmacokinetics. We conducted a pharmacogenetic substudy of patients randomized to sildenafil (n = 85) in the RELAX trial, which investigated the impact of high-dose sildenafil in patients with heart failure with preserved left ventricular ejection fraction (HFpEF). In the overall population, the CYP3A4 inferred phenotype appeared associated with the dose-adjusted peak concentrations of sildenafil at week 12 and week 24 (adjusted P=0.045 for repeated measures analysis), although this P value did not meet our corrected significance threshold of 0.0167. In the more homogeneous Caucasian subgroup, this association was significant (adjusted P=0.0165 for repeated measures). Hence, CYP3A4 inferred phenotype is associated with peak sildenafil dose-adjusted concentrations in patients with HFpEF receiving high doses of sildenafil. The clinical impact of this association requires further investigation.

BODY.INTRODUCTION:
The phosphodiesterase 5 inhibitor (PDE5) sildenafil (Viagra®) has been one of the most widely prescribed drugs in the world for nearly two decades. Despite the established inter-individual variability of its pharmacokinetics and efficacy,1–3 sildenafil has been the subject of only a few pharmacogenomic (PGx) investigations. Limited and inconsistent data from small to moderate size studies have suggested that genetic variants related to a variety of biological pathways could influence sildenafil's hemodynamic effects in patients with pulmonary hypertension or its efficacy in treating erectile dysfunction.4–9  Considering that a majority of the markers judged to be clinically useful by the National Institutes of Health's Clinical PGx Implementation Consortium are absorption, distribution, metabolism and elimination (ADME) genes,10 it is surprising that few studies have focused on variants of genes which could modulate sildenafil concentrations and dosing requirements. In particular, given that sildenafil is extensively metabolized in the liver by the cytochrome P450 3A (CYP3A) isoenzymes (major metabolizing route; 79%) and CYP2C9 (20%),11–13 genetic variants coding for these isoenzymes would biologically appear to be likely genetic modulators of the effects of sildenafil. Although small studies (n < 25) have investigated whether CYP3A5 and CYP2C9 were associated with sildenafil concentrations and pharmacokinetics,2,14 the small sample size from these studies limits their statistical power to identify significant associations.  The RELAX (Phosphodiesterase-5 Inhibition to Improve Clinical Status and Exercise Capacity in Heart Failure with Preserved Ejection Fraction) trial, which investigated the impact of high-dose sildenafil on the exercise capacity and clinical status of patients with heart failure with preserved left ventricular ejection fraction (HFpEF), represents a unique opportunity to further explore the genetic determinants of serum concentrations of sildenafil in a larger population than previous studies. Indeed, as part of RELAX, peak sildenafil concentrations were measured after twelve and twenty-four weeks of treatment.15 Because no beneficial hemodynamic or remodelling effects were observed in RELAX, the aim of this pharmacogenetic sub-study was to identify predictors of sildenafil peak concentrations. Based on existing evidence,2,14,16 the primary goal of this ancillary study was to investigate the impact of variants in the CYP2C9, CYP3A4 and CYP3A5 genes on dose-adjusted peak concentrations of sildenafil measured after 12 and 24 weeks of treatment. We hypothesized that patients carrying genotypes associated with a greater metabolizing capacity for these isoenzymes, such as extensive metabolizers (EM), would present lower dose-adjusted peak concentrations than carriers of genetic variations associated with a lower metabolizing capacity, such as intermediate metabolizers (IM) or poor metabolizers (PM).

BODY.METHODS.OVERVIEW OF STUDY DESIGN:
The methods and results of the RELAX trial (clinicaltrials.gov Identifier: NCT00763867) have been reported previously.15,17 Briefly, RELAX was a multicenter randomized placebo-controlled trial which investigated the impact of high-dose sildenafil on exercise tolerance in patients with HFpEF (LVEF ≥ 50% in the last 12 months presenting New York Heart Association functional class II through IV) whose symptoms were stable while receiving medical therapy.15 The use of significant CYP3A4 inhibitors (e.g. ketoconazole, erythromycin), as well as a current or anticipated future need for nitrate therapy, were exclusion criteria.  Sildenafil was administered orally at 20 mg three times a day for 12 weeks. After 12 weeks, study endpoints were measured as previously described15,17 including peak sildenafil concentrations, which were obtained through phlebotomy 45 to 120 minutes after the scheduled dose. Following this, if the 20 mg three times a day dose was well tolerated, it was then increased to 60 mg three times a day for 12 weeks; otherwise the dose was maintained at 20 mg three times a day. Study endpoint measurements were repeated at week 24, including peak sildenafil concentrations. Sildenafil concentrations were measured by liquid chromatography tandem mass spectrometry (LC-MS/MS) as previously reported.15 Plasma cyclic GMP (cGMP) levels were measured as previously reported at week 24.15 For the current sub–study, we limited our investigations to the 85 patients in the sildenafil group who provided informed consent to participate in the genomic/pharmacogenomic sub-study and who provided a genetic sample.

BODY.METHODS.GENETIC ANALYSES:
The methods of the Heart Failure Network (HFN) genomics/pharmacogenomics sub-studies, including DNA extraction, genotyping and quality control have been previously reported.18 The genotyping strategy included multiple commercial and custom platforms. Among the commercial platforms, we used Sequenom's iPLEX® ADME PGx Panel (Sequenom [now Agena Bioscience], San Diego, CA, USA) to genotype functional SNPs related to the absorption, distribution, metabolism and excretion (ADME). For the current report, we limited our investigations to the variants included on Sequenom's iPLEX® ADME PGx Panel which contains 192 genetic variants, including 183 SNPs which passed genotyping quality control. Following additional data clean-up specific to the group of participants and SNPs of the current analysis (see supplementary information), 75 SNPs with a minor allele frequency (MAF) of ≥0.01 (from 32 genes; see supplementary information) were included.

BODY.METHODS.METABOLIZER STATUS INFERENCE:
Genetic variants genotyped on the Sequenom ADME panel were used to infer the metabolizer status of the primary genes of interest (CYP2C9, CYP3A4, CYP3A5) according to published and recognized scientific evidence in the field (see supplementary information). This approach has the advantage of combining the information of multiple SNPs from a given gene into one variable, thus reducing the number of statistical tests and the severity adjustment made to control for multiple comparisons. This approach was also used for other selected secondary genes (CYP2D6, CYP2C19, CYP2E1, CYP2B6, CYP2A6, CYP2C8, DPYD, SLCO1B1, UGT1A1, SUL1A1, NAT1; see supplementary information), although all secondary genes were also tested as individual SNPs, given that, with the exception of CYP2D6 and CYP2C19, more limited consensus data is available regarding the functional consequences of variants in these genes. Based on previous evidence, and depending on the gene, the metabolizer status was defined, for example, as poor metabolizer (PM), intermediate metabolizer (IM), extensive metabolizer (EM), or ultra metabolizer (UM).

BODY.METHODS.STUDY OBJECTIVES AND ENDPOINTS:
The primary objective of this RELAX pharmacogenetic sub-study was to test the association between the predicted phenotype for CYP2C9, CYP3A4 and CYP3A5 and the concentration:dose ratio of sildenafil measured at 12 and 24 weeks. Similar ratios have been previously used to take into account the impact of different doses on concentrations measured.19–22 Another important factor that motivated our selection of this phenotype as a primary outcome is that it was measured on 2 separate occasions in RELAX. Indeed, our group has used such repeated measures for analyses in previous pharmacogenomic studies21,23 to improve statistical power.  The main secondary objective was to evaluate the association between these three genes and concentrations of cGMP at week 24 in patients still receiving sildenafil. As part of the primary RELAX trial, sildenafil was shown to significantly increase cGMP concentrations compared to baseline, although this difference was not statistically significant compared with placebo.  Other secondary objectives included the evaluation of the association between these endpoints and other ADME genes included on the Sequenom iPlex ADME panel. Finally, as an exploratory analysis, we planned to investigate the relationship of variants significantly associated with the concentration: dose ratio of sildenafil with the occurrence of hypotension and flushing, which were more frequent in the sildenafil group than the placebo group in RELAX.

BODY.METHODS.STATISTICAL ANALYSIS:
For continuous variables, descriptive statistics are presented as mean ± standard deviation, while categorical variables are presented as counts and percentages of each category. Prior to the statistical analyses, the normality test was performed on each endpoint using Shapiro-Wilk test at the 0.01 significance level. As this test appeared to be significant, the normality assumption was rejected and the data was transformed using the log10, square root, inverse and the square functions. The square root of sildenafil peak concentration, the log10 of concentration/dose ratio of sildenafil, and the log10 of cGMP were selected to be the best functions to be statistically analyzed. Data are presented as non-transformed.

BODY.METHODS.GENETIC ASSOCIATION ANALYSES:
The concentration: dose ratio was statistically analyzed using the MIXED procedure of SAS, first including only non-genetic covariates into the analysis model. The repeating factor was the study visit. For cGMP, we used a GLM regression model that, again, was first fitted including non–genetic covariates into the analysis. For the risk of hypotension or flushing, a similar approach was utilized using a logistic regression. For all models, the non-genetic variables investigated were age, sex, body mass index, estimated glomerular filtration rate (eGFR) and amiodarone use. A univariate step selection was conducted between each endpoint and each covariate at 0.05 significance level. When more than one covariate was found to be significantly associated to the endpoint, a stepwise selection method was then performed between each endpoint and these clinical covariates with a P-value of less than 0.1 (entry P-value=0.1, stay P–value≤0.15).  The most important principal components (PC) of the population structure (components 1–3) were included in all genetic analyses to control for genetic ancestry and to avoid confounding by population structure,24 whether or not they were significant at the stepwise level. Thus, all the multivariate genetic models for each of the endpoints were adjusted for the non-genetic covariate meeting the aforementioned criteria and genetic ancestry (PC).  For genes, the inferred metabolizer status was defined, depending on the gene, as poor metabolizer (PM), intermediate metabolizer (IM), extensive metabolizer (EM), or ultra metabolizer (UM) and coded respectively as 0, 1, 2, and 3. For the individual SNP analyses, genotype was coded as 0, 1, or 2 according to the number of copies of minor alleles.  As sensitivity analyses, primary genes and any other genes or SNPs meeting our pre–defined statistical threshold were then analyzed in Caucasians only. Statistical analyses were performed using SAS, PLINK, and existing scripts.

BODY.METHODS.SIGNIFICANCE THRESHOLD:
Statistical tests performed were two-sided and adjusted to account for the multiple testing of genes. For the three genes for the primary analysis (CYP2C9, CYP3A4 and CYP3A5) we used a Bonferroni-adjusted p value threshold of 0.05/3=0.0167. For all other variants, the significance threshold was adjusted using a Bonferroni correction based on the effective number of independent tests (Meff) of all candidate genes, including those of objective 1. The Meff was computed using the method of Gao et al25 and was set at 49. Then, the adjusted significance threshold was calculated as 0.05/49 = 0.001.

BODY.METHODS.ETHICAL CONSIDERATION:
The RELAX study was approved by all local Institutional Review Boards (IRBs) and all patients provided a written consent to participate in the study. Moreover, in centers electing to participate in the genomics/pharmacogenomics sub-study of the HFN clinical trials, the sub-study was approved by the IRB at each side. Patients taking part in the RELAX trial in these participating centers were offered the possibility to participate in the genetics sub-study. All participants willing to participate in the study provided written informed consent.

BODY.RESULTS:
A total of 85 patients were included as part of this sub-study (Table 1). Study subjects were more likely to be males (54%), Caucasian (92%) and were 67.9±10.5 years old. The predicted phenotypes and related genotypes of the study participants for CYP2C9, CYP3A4, CYP3A5 are detailed in Table 2.

BODY.RESULTS.PREDICTORS OF DOSE-CONCENTRATION RATIO OF SILDENAFIL:
Of the 85 patients enrolled in this sub-study, 77 individuals had at least one sildenafil concentration measured (week 12: n = 71; week 24: n = 66). In univariate analysis, both age and BMI were significantly associated with concentration: dose ratios of sildenafil (all P < 0.0001).  In regards to the primary genes of interest, in the overall population, consistent with their predicted lower metabolising capacity, the IM for the CYP3A4 inferred phenotype presented a higher ratio of sildenafil than EM at both week 12 and week 24 (Figure 1, P = 0.045 for repeated measures analysis, after adjusting for age, BMI and three PCs for genetic ancestry). All IM were CYP3A4*22 carriers. Although this association p value was below the nominal threshold of 0.05, it did not meet our corrected significance threshold of 0.0167. Nevertheless, in the more homogeneous Caucasian subgroup, this association was found to be significant after adjusting for multiple comparisons (P = 0.0165 for repeated measures analysis, after adjusting for age, BMI and three PCs for genetic ancestry [Figure 1]).  No association was observed between the CYP2C9 inferred phenotype and the sildenafil concentration:dose ratio (Figure 2). In regards to the CYP3A5 inferred phenotype, although PM tended to have higher concentrations, particularly in the Caucasian subgroup, P values did not reach our significance corrected threshold in the overall population nor in the specific Caucasian subgroup (all P values > 0.14, Figure 3). In an exploratory analysis, we grouped Caucasian patients based on their inferred CYP3A4 and CYP3A5 phenotypes. We observed that increasing levels of predicted metabolizing capacity appeared to be associated with decreasing levels of dose-adjusted sildenafil concentrations (F test P = 0.0238; Figure 4). None of the pairwise comparisons between specific groups reached the pre-established threshold.

BODY.RESULTS.CYCLIC GMP:
After adjusting for age, BMI, baseline cGMP and three PC of genetic ancestry, none of the primary genes were associated with cGMP at week 24.

BODY.RESULTS.SECONDARY GENES:
None of the secondary genetic variants or inferred phenotypes reached our pre-established threshold after multiple comparisons (P = 0.001).

BODY.RESULTS.VASCULAR ADVERSE EFFECTS:
Given the association between the CYP3A4 inferred phenotype and the dose-adjusted concentrations of sildenafil, we explored the association of this genotype with the risk of hypotension and flushing. We observed no significant association, neither in the overall population (IM: 25% vs EM: 17%; P = 0.15), nor with the Caucasian subgroup (IM: 27% vs EM: 18%; P = 0.82).

BODY.DISCUSSION:
In the current study, we investigated the association between three candidate drug metabolism genes and peak dose-adjusted concentrations of sildenafil in HFpEF patients receiving sildenafil. To the best of our knowledge, this pharmacogenetic sub-study is the largest to have investigated the association between candidate genes and concentrations of sildenafil. Furthermore, we believe we are reporting for the first time that the inferred CYP3A4 phenotype, which in this population corresponded to carriers of the CYP3A4*22 variant, may be associated with dose–adjusted peak concentrations of sildenafil. While this association was only numerically apparent in our entire population, it reached our predetermined significance threshold in the Caucasian subgroup. Although we saw numerically higher dose-adjusted concentrations in CYP3A5 PM, particularly in Caucasians, our results did not reach statistical significance. Moreover we found no indication that genetic variations in the CYP2C9-coding gene (CYP2C9), a minor elimination pathway of sildenafil, influenced sildenafil peak concentrations.  CYP3A4 is the most abundant CYP in the liver.26 It contributes to the metabolism of approximately 30% of all drugs, including sildenafil,13,26,27 through oxidative biotransformation.26 Other substrates of CYP3A4 include the immunosuppressant drug tacrolimus, the lipid-lowering drug atorvastatin and endoxifen, the principal active metabolite of the anti-neoplastic agent tamoxifen.27,28 Existing data suggest that sildenafil undergoes significant first-pass metabolism in the liver and that CYP3A4 is the isoenzyme responsible for a majority of this metabolism,12,13 although CYP3A5 may also contribute to the metabolism of sildenafil by the CYP3A subfamily.16
CYP3A4 is also expressed in the gut wall and this pre-systemic metabolism is also thought to limit the bioavailability (40%) of sildenafil.13  The CYP3A4-coding gene, CYP3A4, is located in a cluster on chromosome 7 along with CYP3A5, CYP3A7, and CYP3A4326,27 Few genetic variants in CYP3A4 have been shown to modulate CYP3A4 activity. The intronic variant rs35599367 (C>T), which is located in intron 6 and which was allocated the *22 allele name by the Human Cytochrome P450 (CYP) Allele Nomenclature Database, is, to our knowledge, the only "common" (MAF > 0.01) SNP convincingly shown to influence CYP3A4 expression.27,29 Existing data suggest that the hepatic CYP3A4 mRNA level and enzyme activity are respectively 1.7- and 2.5-fold greater in individuals with the CC genotype than in carriers of the CT and TT carriers.29
CYP3A4*22's MAF is reported to be of approximately 4% in Asians and African-Americans and 8% in Caucasians.27 The frequency in our Caucasian subgroup (7%) is thus consistent with existing data. Our observation that CYP3A4*22 carriers present high dose-adjusted concentration of sildenafil is consistent with other reports showing that CYP3A4*22 carriers present higher concentrations of CYP3A4 substrates, including atorvastatin and tacrolimus.30–32  Despite the fact that we observed a significant association between CYP3A4 and sildenafil concentrations, this did not translate into an association with its downstream biomarker cGMP. The hemodynamic effect of sildenafil and other PDE5 inhibitors in the treatment of erectile dysfunction and pulmonary arterial hypertension is secondary to the inhibition of cGMP degradation by PDE5, which leads to smooth muscle relaxation.11,13 Although this lack of association may appear surprising, it is consistent with the findings from the RELAX trial where sildenafil increased cGMP when compared to baseline but was not statistically significant compared to placebo.15 Moreover, the correlation between sildenafil peak concentrations and cGMP was weak in RELAX, perhaps reflecting that multiple other PDE are implicated in the metabolism of cGMP.15 Nevertheless, this lack of association underlines the uncertainty as to whether the observations in the current and other reports30–32 with CYPA4*22 can be extended to the clinical response to drugs.  Given the high degree of similarity in DNA sequence between CYP3A4 and CYP3A5, the known substrate overlap which exist between these isoenzymes,26 and existing data regarding the implication of CYP3A5 in the metabolism of sildenafil,16 we hypothesised that CYP3A5 could also influence the dose-adjusted peak sildenafil concentrations. In a similar fashion to a previous report in 21 healthy Korean males, we observed higher peak dose-adjusted concentrations in Caucasians participants who were CYP3A5 PM in the present study.2 Nevertheless, as with that report, these trends were not statistically significant, which could be the result of our small sample size (see supplementary material).28 Whether combining information from CYP3A4 and CYP3A5 further allows do decipher inter-individual differences in sildenafil's pharmacokinetics requires additional investigation. If our preliminary observations are confirmed, given the apparent 2.5 fold difference in dose-adjusted concentrations between patients presenting with a CYP3A4 IM / CYP3A5 PM profile and those with a CYP3A4 EM / CYP3A5 IM profile, genotyping could be useful to personalise sildenafil dosing in the clinic.  We did not observe any association between sildenafil concentrations and the CYP2C9 inferred phenotype, a minor route of sildenafil metabolism (20%).11,12 This is consistent with a prior study of 23 healthy males.14 Given the limited metabolism of sildenafil by CYP2C9 but the know functional impact of CYP2C9 variants on the pharmacokinetics of multiple drugs, including warfarin,33 we cannot exclude the possibility that CYP2C9 could exert a modest impact on the pharmacokinetics of sildenafil. Demonstrating such an association would likely necessitate a much larger sample size than the one of the current investigation.

BODY.DISCUSSION.LIMITATIONS:
Our results should be interpreted while taking into consideration some inherent study limitations. First, our results are issued from the RELAX clinical trial which relied on three times per day dosing of sildenafil. The generalizability of our results to dosing schemes commonly used to treat erectile dysfunction requires further investigation.13 Nevertheless, the RELAX dosing regimen is similar to those investigated in pulmonary artery hypertension11 and thus, the potential clinical consequences of our findings may be of particular interest to investigate in these patients. Second, there was no measure or control for the types of food that may have been taken with sildenafil. Ingesting sildenafil with a high-fat meal is known to reduce its absorption rate and peak concentration,12 while foods that can modify CYP3A4 activity such as grapefruit or pomelo juice could increase concentrations.11,12 Despite this, and the inherent heterogeneity in multicenter studies, we were able to identify a novel genetic association between the predicted CYP3A4 phenotype, ultimately CYP3A4*22, and sildenafil dose-adjusted peak concentrations, the association would be expected to be even stronger in a study controlling for nutritional modifiers of the effect. Third, one could argue that because CYP2C9 inhibitors were allowed as part of the study, this could potentially have contributed to our inability to detect genetic association with CYP2C9 variants. As part of RELAX, restriction on CYP2C9 inhibitors were not required given evidence suggesting that CYP2C9 inhibitors exert no clinically significant impact on the pharmacokinetics of sildenafil.11,12 Consistent with this premise, the use of amiodarone, a known CYP2C9 inhibitor,26 was not associated with our primary endpoint. A fourth limitation stems from the single concentration measure available at each visit, which limited our ability to conduct a more detailed pharmacokinetic study of sildenafil to assess other pharmacokinetic parameters of interest such as bioavailability, the elimination half-life or the area under the curve are influenced by CYP3A4. Fifth, despite constituting the largest pharmacogenetic study investigating the pharmacokinetics of sildenafil, our sample size was not sufficiently large to identify variants or inferred phenotypes with a more modest effect. Finally, it should be noted that there is a possibility that some patients may have carried rarer genetic variants that were not tested as part of our assay which could have contributed to misclassification of their metabolizer status. This may even be more likely for non-Caucasian patients.  In conclusion, we observed that the CYP3A4 inferred phenotype, which reflected the CYP3A4*22 carrier status, was associated with peak sildenafil dose adjusted concentrations in patients with HFpEF receiving high doses of sildenafil. Whether these results can be replicated in other populations for which sildenafil is currently indicated and influence its clinical response requires further investigation. Should this observation be validated, genotype-guided prescription of sildenafil in clinical practice or clinical trials could lead to a more predictable clinical response from this agent.

BODY.SUPPLEMENTARY MATERIAL:
1

Question: Compared to CYP2C9 Extensive metabolizers  what was the result of CYP2C9 Poor metabolizers  on Sildenafil levels?
Evidence: No association was observed between the CYP2C9 inferred phenotype and the sildenafil concentration:dose ratio
Answer: no significant difference