Pharmacogenomic Next-Generation DNA Sequencing: Lessons from the Identification and Functional Characterization of Variants of Unknown Significance in CYP2C9 and CYP2C19.

CYP2C9 and CYP2C19 are highly polymorphic pharmacogenes; however, clinically actionable genetic variability in drug metabolism due to these genes has been limited to a few common alleles. The identification and functional characterization of less-common open reading frame sequence variation might help to individualize therapy with drugs that are substrates for the enzymes encoded by these genes. The present study identified seven uncharacterized variants each in CYP2C9 and CYP2C19 using next-generation sequence data for 1013 subjects, and functionally characterized the encoded proteins. Constructs were created and transiently expressed in COS-1 cells for the assay of protein concentration and enzyme activities using fluorometric substrates and liquid chromatography- tandem mass spectrometry with tolbutamide (CYP2C9) and (S)-mephenytoin (CYP2C19) as prototypic substrates. The results were compared with the SIFT, Polyphen, and Provean functional prediction software programs. Cytochrome P450 oxidoreductase (CPR) activities were also determined. Positive correlations were observed between protein content and fluorometric enzyme activity for variants of CYP2C9 (P < 0.05) and CYP2C19 (P < 0.0005). However, CYP2C9 709G>C and CYP2C19 65A>G activities were much lower than predicted based on protein content. Substrate intrinsic clearance values for CYP2C9 218C>T, 343A>C, and CYP2C19 337G>A, 518C>T, 556C>T, and 557G>A were less than 25% of wild-type allozymes. CPR activity levels were similar for all variants. In summary, sequencing of CYP2C9 and CYP2C19 in 1013 subjects identified low-frequency variants that had not previously been functionally characterized. In silico predictions were not always consistent with functional assay results. These observations emphasize the need for high-throughput methods for pharmacogene variant mutagenesis and functional characterization.

A phase 1 study of eribulin mesylate (E7389), a novel microtubule-targeting chemotherapeutic agent, in children with refractory or recurrent solid tumors: A Children's Oncology Group Phase 1 Consortium study (ADVL1314).

BACKGROUND: Eribulin mesylate is a novel anticancer agent that inhibits microtubule growth, without effects on shortening, and promotes nonproductive tubulin aggregate formation. We performed a phase 1 trial to determine the dose-limiting toxicities (DLTs), maximum tolerated or recommended phase 2 dose (MTD/RP2D), and pharmacokinetics (PK) of eribulin in children with refractory or recurrent solid (excluding central nervous system) tumors. METHODS: Eribulin was administered intravenously on days 1 and 8 in 21-day cycles. Three dose levels (1.1, 1.4, and 1.8 mg/m2 /dose) were evaluated using the rolling six design with additional patients enrolled into a PK expansion cohort at the MTD. PK samples were obtained following the day 1, cycle 1 dose. RESULTS: Twenty-three patients, ages 3-17 (median 14) years were enrolled; 20 were evaluable for toxicity. DLTs occurred in 0/6 and 1/6 subjects at the 1.1 and 1.4 mg/m2 /dose, respectively. One subject at the 1.4 mg/m2 /dose had grade 4 neutropenia and grade 3 fatigue. At the 1.8 mg/m2 /dose, 2/5 subjects experienced dose-limiting (grade 4) neutropenia. Grade 3/4 non-DLTs included lymphopenia and hypokalemia, while low-grade toxicities included anorexia and nausea. No episodes of grade > 2 corrected QT interval prolongation or peripheral neuropathy were reported. Eribulin pharmacokinetic parameters were highly variable; the median elimination half-life was 39.6 (range 24.2-96.4) hr. A partial response was observed in one patient (Ewing sarcoma). CONCLUSIONS: Eribulin was well tolerated in children with refractory or recurrent solid tumors with neutropenia identified as the primary DLT. The RP2D of eribulin is 1.4 mg/m2 /dose on days 1 and 8 of a 21-day cycle.

CYP2C9 and CYP2C19: Deep Mutational Scanning and Functional Characterization of Genomic Missense Variants.

Single nucleotide variants in the open reading frames (ORFs) of pharmacogenes are important causes of interindividual variability in drug response. The functional characterization of variants of unknown significance within ORFs remains a major challenge for pharmacogenomics. Deep mutational scanning (DMS) is a high-throughput technique that makes it possible to analyze the functional effect of hundreds of variants in a parallel and scalable fashion. We adapted a "landing pad" DMS system to study the function of missense variants in the ORFs of cytochrome P450 family 2 subfamily C member 9 (CYP2C9) and cytochrome P450 family 2 subfamily C member 19 (CYP2C19). We studied 230 observed missense variants in the CYP2C9 and CYP2C19 ORFs and found that 19 of 109 CYP2C9 and 36 of 121 CYP2C19 variants displayed less than ~ 25% of the wild-type protein expression, a level that may have clinical relevance. Our results support DMS as an efficient method for the identification of damaging ORF variants that might have potential clinical pharmacogenomic application.

Anastrozole Aromatase Inhibitor Plasma Drug Concentration Genome-Wide Association Study: Functional Epistatic Interaction Between SLC38A7 and ALPPL2.

Anastrozole is a widely prescribed aromatase inhibitor for the therapy of estrogen receptor positive (ER+) breast cancer. We performed a genome-wide association study (GWAS) for plasma anastrozole concentrations in 687 postmenopausal women with ER+ breast cancer. The top single-nucleotide polymorphism (SNP) signal mapped across SLC38A7 (rs11648166, P = 2.3E-08), which we showed to encode an anastrozole influx transporter. The second most significant signal (rs28845026, P = 5.4E-08) mapped near ALPPL2 and displayed epistasis with the SLC38A7 signal. Both of these SNPs were cis expression quantitative trait loci (eQTL)s for these genes, and patients homozygous for variant genotypes for both SNPs had the highest drug concentrations, the highest SLC38A7 expression, and the lowest ALPPL2 expression. In summary, our GWAS identified a novel gene encoding an anastrozole transporter, SLC38A7, as well as epistatic interaction between SNPs in that gene and SNPs near ALPPL2 that influenced both the expression of the transporter and anastrozole plasma concentrations.