Liver eQTL meta-analysis illuminates potential molecular mechanisms of cardiometabolic traits.

Understanding the molecular mechanisms of complex traits is essential for developing targeted interventions. We analyzed liver expression quantitative-trait locus (eQTL) meta-analysis data on 1,183 participants to identify conditionally distinct signals. We found 9,013 eQTL signals for 6,564 genes; 23% of eGenes had two signals, and 6% had three or more signals. We then integrated the eQTL results with data from 29 cardiometabolic genome-wide association study (GWAS) traits and identified 1,582 GWAS-eQTL colocalizations for 747 eGenes. Non-primary eQTL signals accounted for 17% of all colocalizations. Isolating signals by conditional analysis prior to coloc resulted in 37% more colocalizations than using marginal eQTL and GWAS data, highlighting the importance of signal isolation. Isolating signals also led to stronger evidence of colocalization: among 343 eQTL-GWAS signal pairs in multi-signal regions, analyses that isolated the signals of interest resulted in higher posterior probability of colocalization for 41% of tests. Leveraging allelic heterogeneity, we predicted causal effects of gene expression on liver traits for four genes. To predict functional variants and regulatory elements, we colocalized eQTL with liver chromatin accessibility QTL (caQTL) and found 391 colocalizations, including 73 with non-primary eQTL signals and 60 eQTL signals that colocalized with both a caQTL and a GWAS signal. Finally, we used publicly available massively parallel reporter assays in HepG2 to highlight 14 eQTL signals that include at least one expression-modulating variant. This multi-faceted approach to unraveling the genetic underpinnings of liver-related traits could lead to therapeutic development.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD).

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Germline cis variant determines epigenetic regulation of the anti-cancer drug metabolism gene dihydropyrimidine dehydrogenase (DPYD).

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Sequencing of genes of drug response in tumor DNA and implications for precision medicine in cancer patients.

Tumor DNA sequencing is becoming standard-of-care for patient treatment decisions. We evaluated genotype concordance between tumor DNA and genomic DNA from blood and catalogued functional effects of somatic mutations in 21 drug response genes in 752 solid tumor patients. Using a threshold of 10% difference between tumor and blood DNA variant allele fraction (VAF), concordance for heterogenous genotype calls was 78% and increased to 97.5% using a 30% VAF threshold. Somatic mutations were observed in all 21 drug response genes, and 44% of patients had at least one somatic mutation in these genes. In tumor DNA, eight patients had a frameshift mutation in CYP2C8, which metabolizes taxanes. Overall, somatic copy number losses were more frequent than gains, including for CYP2C19 and CYP2D6 which had the most frequent copy number losses. However, copy number gains in TPMT were more than four times as common as losses. Seven % of patients had copy number gains in ABCB1, a multidrug resistance transporter of anti-cancer agents. These results demonstrate tumor-only DNA sequencing might not be reliable to call germline genotypes of drug response variants.

Identification of Reduced ERAP2 Expression and a Novel HLA Allele as Components of a Risk Score for Susceptibility to Liver Injury Due to Amoxicillin-Clavulanate.

BACKGROUND & AIMS: Drug-induced liver injury (DILI) due to amoxicillin-clavulanate (AC) has been associated with HLA-A∗02:01, HLA-DRB1∗15:01, and rs2476601, a missense variant in PTPN22. The aim of this study was to identify novel risk factors for AC-DILI and to construct a genetic risk score (GRS). METHODS: Transcriptome-wide association study and genome-wide association study analyses were performed on 444 AC-DILI cases and 10,397 population-based controls of European descent. Associations were confirmed in a validation cohort (n = 133 cases and 17,836 population-based controls). Discovery and validation AC-DILI cases were also compared with 1358 and 403 non-AC-DILI cases. RESULTS: Transcriptome-wide association study revealed a significant association of AC-DILI risk with reduced liver expression of ERAP2 (P = 3.7 × 10-7), coding for an aminopeptidase involved in antigen presentation. The lead eQTL single nucleotide polymorphism, rs1363907 (G), was associated with AC-DILI risk in the discovery (odds ratio [OR], 1.68; 95% CI, 1.23-1.66; P = 1.7 × 10-7) and validation cohorts (OR, 1.2; 95% CI, 1.04-2.05; P = .03), following a recessive model. We also identified HLA-B∗15:18 as a novel AC-DILI risk factor in both discovery (OR, 4.19; 95% CI, 2.09-8.36; P = 4.9 × 10-5) and validation (OR, 7.78; 95% CI, 2.75-21.99; P = .0001) cohorts. GRS, incorporating rs1363907, rs2476601, HLA-B∗15:18, HLA-A∗02:01, and HLA-DRB1∗15:01, was highly predictive of AC-DILI risk when cases were analyzed against both general population and non-AC-DILI control cohorts. GRS was the most significant predictor in a regression model containing known AC-DILI clinical risk characteristics and significantly improved the predictive model. CONCLUSIONS: We identified novel associations of AC-DILI risk with ERAP2 low expression and with HLA-B∗15:18. GRS based on the 5 risk variants may assist AC-DILI causality assessment and risk management.

A resource for integrated genomic analysis of the human liver.

In this study, we generated whole-transcriptome RNA-Seq from n = 192 genotyped liver samples and used these data with existing data from the GTEx Project (RNA-Seq) and previous liver eQTL (microarray) studies to create an enhanced transcriptomic sequence resource in the human liver. Analyses of genotype-expression associations show pronounced enrichment of associations with genes of drug response. The associations are primarily consistent across the two RNA-Seq datasets, with some modest variation, indicating the importance of obtaining multiple datasets to produce a robust resource. We further used an empirical Bayesian model to compare eQTL patterns in liver and an additional 20 GTEx tissues, finding that MHC genes, and especially class II genes, are enriched for liver-specific eQTL patterns. To illustrate the utility of the resource to augment GWAS analysis with small sample sizes, we developed a novel meta-analysis technique to combine several liver eQTL data sources. We also illustrate its application using a transcriptome-enhanced re-analysis of a study of neutropenia in pancreatic cancer patients. The associations of genotype with liver expression, including splice variation and its genetic associations, are made available in a searchable genome browser.

Identification and single-base gene-editing functional validation of a cis-EPO variant as a genetic predictor for EPO-increasing therapies.

Hypoxia-inducible factor prolyl hydroxylase inhibitors (HIF-PHIs) are currently under clinical development for treating anemia in chronic kidney disease (CKD), but it is important to monitor their cardiovascular safety. Genetic variants can be used as predictors to help inform the potential risk of adverse effects associated with drug treatments. We therefore aimed to use human genetics to help assess the risk of adverse cardiovascular events associated with therapeutically altered EPO levels to help inform clinical trials studying the safety of HIF-PHIs. By performing a genome-wide association meta-analysis of EPO (n = 6,127), we identified a cis-EPO variant (rs1617640) lying in the EPO promoter region. We validated this variant as most likely causal in controlling EPO levels by using genetic and functional approaches, including single-base gene editing. Using this variant as a partial predictor for therapeutic modulation of EPO and large genome-wide association data in Mendelian randomization tests, we found no evidence (at p < 0.05) that genetically predicted long-term rises in endogenous EPO, equivalent to a 2.2-unit increase, increased risk of coronary artery disease (CAD, OR [95% CI] = 1.01 [0.93, 1.07]), myocardial infarction (MI, OR [95% CI] = 0.99 [0.87, 1.15]), or stroke (OR [95% CI] = 0.97 [0.87, 1.07]). We could exclude increased odds of 1.15 for cardiovascular disease for a 2.2-unit EPO increase. A combination of genetic and functional studies provides a powerful approach to investigate the potential therapeutic profile of EPO-increasing therapies for treating anemia in CKD.

Polygenic Risk Scores for Blood Pressure to Assess the Risk of Severe Bevacizumab-Induced Hypertension in Cancer Patients (Alliance).

Hypertension is a common bevacizumab-induced toxicity. No markers are available to predict patients at risk of developing hypertension. We hypothesized that genetic risk of essential hypertension, as measured by a blood pressure polygenic risk score (PRS), would be associated with risk of severe bevacizumab-induced hypertension. PRSs were calculated for 1,027 bevacizumab-treated patients of European descent with cancer from four clinical trials (Alliance for Clinical Trials in Oncology (Alliance) / Cancer and Leukemia Group B (CALGB) 80303, 40503, 90401, 40502) using summary systolic blood pressure (SBP) and diastolic blood pressure (DBP) genome-wide association results obtained from 757,601 individuals of European descent. The association between PRS and grade 3 bevacizumab-induced hypertension (Common Toxicity Criteria for Adverse Events version 3) in each trial was performed by multivariable logistic regression. Fixed-effect meta-analyses odds ratios (ORs) per standard deviation (SD) of the association of PRS (quantitative) and hypertension across trials were estimated by inverse-variance weighting. PRSs were additionally stratified into quintiles, with the bottom quintile as the referent group. The OR of the association between hypertension and each quintile vs. the referent group was determined by logistic regression. The most significant PRS (quantitative)-hypertension association included up to 67 single-nucleotide variants (SNPs) associated with SBP (P = 0.0077, OR per SD = 1.31, 95% confidence interval (CI), 1.07-1.60), and up to 53 SNPs associated with DBP (P = 0.0209, OR per SD = 1.27, 95% CI, 1.04-1.56). Patients in the top quintile had a higher risk of developing bevacizumab-induced hypertension compared with patients in the bottom quintile using SNPs associated with SBP (P = 4.75 × 10-4 , OR = 3.72, 95% CI, 1.84-8.16) and DBP (P = 0.076, OR = 1.83, 95% CI, 0.95-3.64). Genetic variants associated with essential hypertension, mainly SBP, increase the risk of severe bevacizumab-induced hypertension.

KDR genetic predictor of toxicities induced by sorafenib and regorafenib.

No biomarkers are available to predict toxicities induced by VEGFR TKIs. This study aimed to identify markers of toxicities induced by these drugs using a discovery-validation approach. The discovery set included 140 sorafenib-treated cancer patients (TARGET study) genotyped for SNPs in 56 genes. The most significant SNPs associated with grade ≥2 hypertension, diarrhea, dermatologic toxicities, and composite toxicity (any one of the toxicities) were tested for association with grade ≥2 toxicity in a validation set of 201 sorafenib-treated patients (Alliance/CALGB 80802). The validated SNP was tested for association with grade ≥2 toxicity in 107 (LCCC 1029) and 82 (Italian cohort) regorafenib-treated patients. SNP-toxicity associations were evaluated using logistic regression, and a meta-analysis between the studies was performed by inverse variance. Variant rs4864950 in KDR increased the risk of grade ≥2 composite toxicity in TARGET, Alliance/CALGB 80802, and the Italian cohort (meta-analysis p = 6.79 × 10-4, OR = 2.01, 95% CI 1.34-3.01). We identified a predictor of toxicities induced by VEGFR TKIs. CLINICALTRIALS.GOV IDENTIFIER: NCT00073307 (TARGET), NCT01015833 (Alliance/CALGB 80802), and NCT01298570 (LCCC 1029).

Genome-wide association studies of survival in 1520 cancer patients treated with bevacizumab-containing regimens.

Germline variants might predict cancer progression. Bevacizumab improves overall survival (OS) in patients with advanced cancers. No biomarkers are available to identify patients that benefit from bevacizumab. A meta-analysis of genome-wide association studies (GWAS) was conducted in 1,520 patients from Phase III trials (CALGB 80303, 40503, 80405 and ICON7), where bevacizumab was randomized to treatment without bevacizumab. We aimed to identify genes and single nucleotide polymorphisms (SNPs) associated with survival independently of bevacizumab treatment or through interaction with bevacizumab. A cause-specific Cox model was used to test the SNP-OS association in both arms combined (prognostic), and the effect of SNPs-bevacizumab interaction on OS (predictive) in each study. The SNP effects across studies were combined using inverse variance. Findings were tested for replication in advanced colorectal and ovarian cancer patients from The Cancer Genome Atlas (TGCA). In the GWAS meta-analysis, patients with rs680949 in PRUNE2 experienced shorter OS compared to patients without it (P = 1.02 × 10-7 , hazard ratio [HR] = 1.57, 95% confidence interval [CI] 1.33-1.86), as well as in TCGA (P = .0219, HR = 1.58, 95% CI 1.07-2.35). In the GWAS meta-analysis, patients with rs16852804 in BARD1 experienced shorter OS compared to patients without it (P = 1.40 × 10-5 , HR = 1.51, 95% CI 1.25-1.82) as well as in TCGA (P = 1.39 × 10-4 , HR = 3.09, 95% CI 1.73-5.51). Patients with rs3795897 in AGAP1 experienced shorter OS in the bevacizumab arm compared to the nonbevacizumab arm (P = 1.43 × 10-5 ). The largest GWAS meta-analysis of bevacizumab treated patients identified PRUNE2 and BARD1 (tumor suppressor genes) as prognostic genes of colorectal and ovarian cancer, respectively, and AGAP1 as a potentially predictive gene that interacts with bevacizumab with respect to patient survival.

Bevacizumab-induced hypertension and proteinuria: a genome-wide study of more than 1000 patients.

BACKGROUND: Hypertension and proteinuria are common bevacizumab-induced toxicities. No validated biomarkers are available for identifying patients at risk of these toxicities. METHODS: A genome-wide association study (GWAS) meta-analysis was performed in 1039 bevacizumab-treated patients of European ancestry in four clinical trials (CALGB 40502, 40503, 80303, 90401). Grade ≥2 hypertension and proteinuria were recorded (CTCAE v.3.0). Single-nucleotide polymorphism (SNP)-toxicity associations were determined using a cause-specific Cox model adjusting for age and sex. RESULTS: The most significant SNP associated with hypertension with concordant effect in three out of the four studies (p-value <0.05 for each study) was rs6770663 (A > G) in KCNAB1, with the G allele increasing the risk of hypertension (p-value = 4.16 × 10-6). The effect of the G allele was replicated in ECOG-ACRIN E5103 in 582 patients (p-value = 0.005). The meta-analysis of all five studies for rs6770663 led to p-value = 7.73 × 10-8, close to genome-wide significance. The most significant SNP associated with proteinuria was rs339947 (C > A, between DNAH5 and TRIO), with the A allele increasing the risk of proteinuria (p-value = 1.58 × 10-7). CONCLUSIONS: The results from the largest study of bevacizumab toxicity provide new markers of drug safety for further evaluations. SNP in KCNAB1 validated in an independent dataset provides evidence toward its clinical applicability to predict bevacizumab-induced hypertension. ClinicalTrials.gov Identifier: NCT00785291 (CALGB 40502); NCT00601900 (CALGB 40503); NCT00088894 (CALGB 80303) and NCT00110214 (CALGB 90401).

PIK3R5 genetic predictors of hypertension induced by VEGF-pathway inhibitors.

No biomarkers are available to predict patients at risk of developing hypertension induced by VEGF-pathway inhibitors. This study aimed to identify predictive biomarkers of hypertension induced by these drugs using a discovery-replication approach. The discovery set included 140 sorafenib-treated patients (TARGET study) genotyped for 973 SNPs in 56 genes. The most statistically significant SNPs associated with grade ≥2 hypertension were tested for association with grade ≥2 hypertension in the replication set of a GWAS of 1039 bevacizumab-treated patients from four clinical trials (CALGB/Alliance). In the discovery set, rs444904 (G > A) in PIK3R5 was associated with an increased risk of sorafenib-induced hypertension (p = 0.006, OR = 3.88 95% CI 1.54-9.81). In the replication set, rs427554 (G > A) in PIK3R5 (in complete linkage disequilibrium with rs444904) was associated with an increased risk of bevacizumab-induced hypertension (p = 0.008, OR = 1.39, 95% CI 1.09-1.78). This study identified a predictive marker of drug-induced hypertension that should be evaluated for other VEGF-pathway inhibitors.ClinicalTrials.gov Identifier:NCT00073307 (TARGET).

Plasma levels of angiopoietin-2, VEGF-A, and VCAM-1 as markers of bevacizumab-induced hypertension: CALGB 80303 and 90401 (Alliance).

Hypertension is a common toxicity induced by bevacizumab and other antiangiogenic drugs. There are no biomarkers to predict the risk of bevacizumab-induced hypertension. This study aimed to identify plasma proteins related to the function of the vasculature to predict the risk of severe bevacizumab-induced hypertension. Using pretreated plasma samples from 398 bevacizumab-treated patients in two clinical trials (CALGB 80303 and 90401), the levels of 17 proteins were measured via ELISA. The association between proteins and grade 3 bevacizumab-induced hypertension was performed by calculating the odds ratio (OR) from logistic regression adjusting for age, sex, and clinical trial. Using the optimal cut-point of each protein, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) for hypertension were estimated. Five proteins showed no difference in levels between clinical trials and were used for analyses. Lower levels of angiopoietin-2 (p = 0.0013, OR 3.41, 95% CI 1.67-7.55), VEGF-A (p = 0.0008, OR 4.25, 95% CI 1.93-10.72), and VCAM-1 (p = 0.0067, OR 2.68, 95% CI 1.34-5.63) were associated with an increased risk of grade 3 hypertension. The multivariable model suggests independent effects of angiopoietin-2 (p = 0.0111, OR 2.71, 95% CI 1.29-6.10), VEGF-A (p = 0.0051, OR 3.66, 95% CI 1.54-9.73), and VCAM-1 (p = 0.0308, OR 2.27, 95% CI 1.10-4.92). The presence of low levels of 2-3 proteins had an OR of 10.06 (95% CI 3.92-34.18, p = 1.80 × 10-5) for the risk of hypertension, with sensitivity of 89.7%, specificity of 53.5%, PPV of 17.3%, and NPV of 97.9%. This is the first study providing evidence of plasma proteins with potential value to predict patients at risk of developing bevacizumab-induced hypertension.Clinical trial registration: ClinicalTrials.gov Identifier: NCT00088894 (CALGB 80303); and NCT00110214 (CALGB 90401).

Integration of DNA sequencing with population pharmacokinetics to improve the prediction of irinotecan exposure in cancer patients.

BACKGROUND: Irinotecan (CPT-11) is an anticancer agent widely used to treat adult solid tumours. Large interindividual variability in the clearance of irinotecan and SN-38, its active and toxic metabolite, results in highly unpredictable toxicity. METHODS: In 217 cancer patients treated with intravenous irinotecan single agent or in combination, germline DNA was used to interrogate the variation in 84 genes by next-generation sequencing. A stepwise analytical framework including a population pharmacokinetic model with SNP- and gene-based testing was used to identify demographic/clinical/genetic factors that influence the clearance of irinotecan and SN-38. RESULTS: Irinotecan clearance was influenced by rs4149057 in SLCO1B1, body surface area, and co-administration of 5-fluorouracil/leucovorin/bevacizumab. SN-38 clearance was influenced by rs887829 in UGT1A1, pre-treatment total bilirubin, and EGFR rare variant burden. Within each UGT1A1 genotype group, elevated pre-treatment total bilirubin and/or presence of at least one rare variant in EGFR resulted in significantly lower SN-38 clearance. The model reduced the interindividual variability in irinotecan clearance from 38 to 34% and SN-38 clearance from 49 to 32%. CONCLUSIONS: This new model significantly reduced the interindividual variability in the clearance of irinotecan and SN-38. New genetic factors of variability in clearance have been identified.

Genetic effects on liver chromatin accessibility identify disease regulatory variants.

Identifying the molecular mechanisms by which genome-wide association study (GWAS) loci influence traits remains challenging. Chromatin accessibility quantitative trait loci (caQTLs) help identify GWAS loci that may alter GWAS traits by modulating chromatin structure, but caQTLs have been identified in a limited set of human tissues. Here we mapped caQTLs in human liver tissue in 20 liver samples and identified 3,123 caQTLs. The caQTL variants are enriched in liver tissue promoter and enhancer states and frequently disrupt binding motifs of transcription factors expressed in liver. We predicted target genes for 861 caQTL peaks using proximity, chromatin interactions, correlation with promoter accessibility or gene expression, and colocalization with expression QTLs. Using GWAS signals for 19 liver function and/or cardiometabolic traits, we identified 110 colocalized caQTLs and GWAS signals, 56 of which contained a predicted caPeak target gene. At the LITAF LDL-cholesterol GWAS locus, we validated that a caQTL variant showed allelic differences in protein binding and transcriptional activity. These caQTLs contribute to the epigenomic characterization of human liver and help identify molecular mechanisms and genes at GWAS loci.

Genomic Analysis of Germline Variation Associated with Survival of Patients with Colorectal Cancer Treated with Chemotherapy Plus Biologics in CALGB/SWOG 80405 (Alliance).

PURPOSE: Irinotecan/5-fluorouracil (5-FU; FOLFIRI) or oxaliplatin/5-FU (FOLFOX), combined with bevacizumab or cetuximab, are approved, first-line treatments for metastatic colorectal cancer (mCRC). We aimed at identifying germline variants associated with survival in patients with mCRC treated with these regimens in Cancer and Leukemia Group B/SWOG 80405. EXPERIMENTAL DESIGN: Patients with mCRC receiving either FOLFOX or FOLFIRI were randomized to either cetuximab or bevacizumab. DNA from peripheral blood was genotyped for approximately 700,000 SNPs. The association between SNPs and overall survival (OS) was tested in 613 patients of genetically estimated European ancestry using Cox proportional hazards models. RESULTS: The four most significant SNPs associated with OS were three haplotypic SNPs between microsomal glutathione S-transferase 1 (MGST1) and LIM domain only 3 (LMO3, representative HR, 1.56; P = 1.30 × 10-6), and rs11644916 in AXIN1 (HR, 1.39, P = 4.26 × 10-6). AXIN1 is a well-established tumor suppressor gene in colorectal cancer, and rs11644916 (G>A) conferred shorter OS. Median OS for patients with the AA, AG, or GG genotypes was 18.4, 25.6, or 36.4 months, respectively. In 90 patients with stage IV colorectal cancer from The Cancer Genome Atlas (TCGA), rs11649255 in AXIN1 [in almost complete linkage disequilibrium (LD) with rs11644916], was associated with shorter OS (HR, 2.24, P = 0.0096). Using rs11648673 in AXIN1 (in very high LD with rs11644916 and with functional evidence), luciferase activity in three colorectal cancer cell lines was reduced. CONCLUSIONS: This is the first large genome-wide association study ever conducted in patients with mCRC treated with first-line standard treatment in a randomized phase III trial. A common SNP in AXIN1 conferred worse OS and the effect was replicated in TCGA. Further studies in colorectal cancer experimental models are required.

Optimal Sampling Strategies for Irinotecan (CPT-11) and its Active Metabolite (SN-38) in Cancer Patients.

Irinotecan (CPT-11) is an anticancer agent widely used in the treatment of a variety of adult solid tumors. The objective of this study was to develop an optimal sampling strategy model that accurately estimates pharmacokinetic parameters of CPT-11 and its active metabolite, SN-38. This study included 221 patients with advanced solid tumors or lymphoma receiving CPT-11 single or combination therapy with 5-fluorouracil (5-FU)/leucovorin (LV) (FOLFIRI) plus bevacizumab from 4 separate clinical trials. Population pharmacokinetic analysis of CPT-11 and SN-38 was performed by non-linear mixed effects modeling. The optimal sampling strategy model was developed using D-optimality with expected distribution approach. The pharmacokinetic profiles of CPT-11 and SN-38 were best described by a 3- and 2-compartment model, respectively, with first-order elimination. Body surface area and co-administration with 5-FU/LV plus bevacizumab were significant covariates (p < 0.01) for volumes of the central compartment of CPT-11 and SN-38, and clearance of CPT-11. Pre-treatment total bilirubin and co-administration with 5-FU/LV and bevacizumab were significant covariates (p < 0.01) for clearance of SN-38. Accurate and precise predictive performance (r2 > 0.99, -2 < bias (%ME) < 0, precision (% RMSE) < 12) of both CPT-11 and SN-38 was achieved using: (i) 6 fixed sampling times collected at 1.5, 3.5, 4, 5.75, 22, 23.5 hours post-infusion; or (ii) 1 fixed time and 2 sampling windows collected at 1.5, [3-5.75], [22-23.5] hours post-infusion. The present study demonstrates that an optimal sampling design with three blood samples achieves accurate and precise pharmacokinetic parameter estimates for both CPT-11 and SN-38.

A New Liver Expression Quantitative Trait Locus Map From 1,183 Individuals Provides Evidence for Novel Expression Quantitative Trait Loci of Drug Response, Metabolic, and Sex-Biased Phenotypes.

Expression quantitative trait locus (eQTL) studies in human liver are crucial for elucidating how genetic variation influences variability in disease risk and therapeutic outcomes and may help guide strategies to obtain maximal efficacy and safety of clinical interventions. Associations between expression microarray and genome-wide genotype data from four human liver eQTL studies (n = 1,183) were analyzed. More than 2.3 million cis-eQTLs for 15,668 genes were identified. When eQTLs were filtered against a list of 1,496 drug response genes, 187,829 cis-eQTLs for 1,191 genes were identified. Additionally, 1,683 sex-biased cis-eQTLs were identified, as well as 49 and 73 cis-eQTLs that colocalized with genome-wide association study signals for blood metabolite or lipid levels, respectively. Translational relevance of these results is evidenced by linking DPYD eQTLs to differences in safety of chemotherapy, linking the sex-biased regulation of PCSK9 expression to anti-lipid therapy, and identifying the G-protein coupled receptor GPR180 as a novel drug target for hypertriglyceridemia.