Computational pharmacogenotype extraction from clinical next-generation sequencing

Background

Next-generation sequencing (NGS), including whole genome sequencing (WGS) and whole exome sequencing (WES), is increasingly being used for clinic care. While NGS data have the potential to be repurposed to support clinical pharmacogenomics (PGx), current computational approaches have not been widely validated using clinical data. In this study, we assessed the accuracy of the Aldy computational method to extract PGx genotypes from WGS and WES data for 14 and 13 major pharmacogenes, respectively.

Methods

Germline DNA was isolated from whole blood samples collected for 264 patients seen at our institutional molecular solid tumor board. DNA was used for panel-based genotyping within our institutional Clinical Laboratory Improvement Amendments- (CLIA-) certified PGx laboratory. DNA was also sent to other CLIA-certified commercial laboratories for clinical WGS or WES. Aldy v3.3 and v4.4 were used to extract PGx genotypes from these NGS data, and results were compared to the panel-based genotyping reference standard that contained 45 star allele-defining variants within CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2D6, CYP3A4, CYP3A5, CYP4F2, DPYD, G6PD, NUDT15, SLCO1B1, TPMT, and VKORC1.

Results

Mean WGS read depth was >30x for all variant regions except for G6PD (average read depth was 29 reads), and mean WES read depth was >30x for all variant regions. For 94 patients with WGS, Aldy v3.3 diplotype calls were concordant with those from the genotyping reference standard in 99.5% of cases when excluding diplotypes with additional major star alleles not tested by targeted genotyping, ambiguous phasing, and CYP2D6 hybrid alleles. Aldy v3.3 identified 15 additional clinically actionable star alleles not covered by genotyping within CYP2B6, CYP2C19, DPYD, SLCO1B1, and NUDT15. Within the WGS cohort, Aldy v4.4 diplotype calls were concordant with those from genotyping in 99.7% of cases. When excluding patients with CYP2D6 copy number variation, all Aldy v4.4 diplotype calls except for one CYP3A4 diplotype call were concordant with genotyping for 161 patients in the WES cohort.

Conclusion

Aldy v3.3 and v4.4 called diplotypes for major pharmacogenes from clinical WES and WGS data with >99% accuracy. These findings support the use of Aldy to repurpose clinical NGS data to inform clinical PGx.

CYP2C9*61, a rare missense variant identified in a Puerto Rican patient with low warfarin dose requirements

Warfarin continues to be the mainstay therapy for preventing thrombus formation. Although pharmacogenetic algorithms have shown higher predictability of the optimal warfarin dose and lower occurrence of bleeding episodes, they often do not include ethno-specific genetic variants relevant to non-Europeans. This case report describes a rare missense variant at exon 9 of CYP2C9 (rs202201137; c.1370A>G transition; p.Asn457Ser) found in a Puerto Rican patient with low warfarin dose requirements (3 mg/day). The haplotype characterized by two amino acid changes, Asn457Ser and Arg144Cys (rs1799853; c.430C>T), has been designated CYP2C9*61 by the Pharmacogene Variation Consortium. According to prediction scores assessed with the Combined Annotation Dependent Depletion tool, CYP2C9*61 (p.Asn457Ser) was classified as nondeleterious, therefore its impact on CYP2C9 enzymatic activity cannot be postulated.