Analytical Validation of a Computational Method for Pharmacogenetic Genotyping from Clinical Whole Exome Sequencing

Return of Clinically Actionable Pharmacogenetic Results From Molecular Tumor Board DNA Sequencing Data: Workflow and Estimated Costs

Pharmacogenetic testing can prevent severe toxicities from several oncology drug therapies; it also has the potential to improve the outcomes from supportive care drugs. Paired tumor and germline sequencing is increasingly common in oncology practice; these include sequencing of pharmacogenes, but the germline pharmacogenetic variants are rarely included in the clinical reports, despite many being clinically actionable. We established an informatics workflow to evaluate the clinical sequencing results for pharmacogenetic variants. We used the Aldy computational tool, which we have previously shown to determine the variant alleles in 14 pharmacogenes in clinical sequencing data with >99% accuracy, to identify pharmacogenetic variants in the clinical whole exome sequencing from our molecular tumor board. Patients with genetic variants that are clinically actionable for their individual therapy programs, including both treatment and supportive care, are referred to a clinical pharmacogenetics testing laboratory for confirmation. Through an evaluation of our weekly informatics workflow, we determined it took approximately 3.25 hours to complete the analysis of the sequencing data from approximately 20 patients. Using a United States pharmacist's median salary, we estimated the incremental added cost of the process to be only ~$15 per patient. This adds only a minor increase to the patient's cost of testing and has the potential to improve the safety and efficacy of their treatment.

Innovation in cancer pharmacotherapy through integrative consideration of germline and tumor genomes

Precision cancer medicine is widely established, and numerous molecularly targeted drugs for various tumor entities are approved or are in development. Personalized pharmacotherapy in oncology has so far been based primarily on tumor characteristics, for example, somatic mutations. However, the response to drug treatment also depends on pharmacological processes summarized under the term ADME (absorption, distribution, metabolism, and excretion). Variations in ADME genes have been the subject of intensive research for >5 decades, considering individual patients' genetic makeup, referred to as pharmacogenomics (PGx). The combined impact of a patient's tumor and germline genome is only partially understood and often not adequately considered in cancer therapy. This may be attributed, in part, to the lack of methods for combined analysis of both data layers. Optimized personalized cancer therapies should, therefore, aim to integrate molecular information, which derives from both the tumor and the germline genome, and taking into account existing PGx guidelines for drug therapy. Moreover, such strategies should provide the opportunity to consider genetic variants of previously unknown functional significance. Bioinformatic analysis methods and corresponding algorithms for data interpretation need to be developed to integrate PGx data in cancer therapy with a special meaning for interdisciplinary molecular tumor boards, in which cancer patients are discussed to provide evidence-based recommendations for clinical management based on individual tumor profiles. SIGNIFICANCE STATEMENT: The era of personalized oncology has seen the emergence of drugs tailored to genetic variants associated with cancer biology. However, the full potential of targeted therapy remains untapped owing to the predominant focus on acquired tumor-specific alterations. Optimized cancer care must integrate tumor and patient genomes, guided by pharmacogenomic principles. An essential prerequisite for realizing truly personalized drug treatment of cancer patients is the development of bioinformatic tools for comprehensive analysis of all data layers generated in modern precision oncology programs.

Geny: a genotyping tool for allelic decomposition of killer cell immunoglobulin-like receptor genes

Introduction

Accurate genotyping of Killer cell Immunoglobulin-like Receptor (KIR) genes plays a pivotal role in enhancing our understanding of innate immune responses, disease correlations, and the advancement of personalized medicine. However, due to the high variability of the KIR region and high level of sequence similarity among different KIR genes, the generic genotyping workflows are unable to accurately infer copy numbers and complete genotypes of individual KIR genes from next-generation sequencing data. Thus, specialized genotyping tools are needed to genotype this complex region.

Methods

Here, we introduce Geny, a new computational tool for precise genotyping of KIR genes. Geny utilizes available KIR allele databases and proposes a novel combination of expectation-maximization filtering schemes and integer linear programming-based combinatorial optimization models to resolve ambiguous reads, provide accurate copy number estimation, and estimate the correct allele of each copy of genes within the KIR region.

Results & Discussion

We evaluated Geny on a large set of simulated short-read datasets covering the known validated KIR region assemblies and a set of Illumina short-read samples sequenced from 40 validated samples from the Human Pangenome Reference Consortium collection and showed that it outperforms the existing state-of-the-art KIR genotyping tools in terms of accuracy, precision, and recall. We envision Geny becoming a valuable resource for understanding immune system response and consequently advancing the field of patient-centric medicine.

Understanding the Biology and Testing Techniques for Pharmacogenomics in Oncology: A Practical Guide for the Clinician

Pharmacogenomic (PGx) testing is a growing area of personalized medicine with demonstrated clinical utility in improving patient outcomes in oncology. PGx testing of pharmacogenes affecting drug pharmacokinetics, pharmacodynamics, and response can help inform drug selection and dosing of several anticancer therapies and supportive care medications. Several PGx testing techniques exist including polymerase chain reaction (PCR), MassARRAY, microarray, and sequencing. This review article provides a clinician-friendly guide of these techniques. Understanding the advantages, limitations, ideal use, and potential clinical applications of each platform can help clinicians choose the appropriate PGx testing platform for specific use cases.

Development of a Multifaceted Program for Pharmacogenetics Adoption at an Academic Medical Center: Practical Considerations and Lessons Learned

In 2019, Indiana University launched the Precision Health Initiative to enhance the institutional adoption of precision medicine, including pharmacogenetics (PGx) implementation, at university-affiliated practice sites across Indiana. The overarching goal of this PGx implementation program was to facilitate the sustainable adoption of genotype-guided prescribing into routine clinical care. To accomplish this goal, we pursued the following specific objectives: (i) to integrate PGx testing into existing healthcare system processes; (ii) to implement drug-gene pairs with high-level evidence and educate providers and pharmacists on established clinical management recommendations; (iii) to engage key stakeholders, including patients to optimize the return of results for PGx testing; (iv) to reduce health disparities through the targeted inclusion of underrepresented populations; (v) and to track third-party reimbursement. This tutorial details our multifaceted PGx implementation program, including descriptions of our interventions, the critical challenges faced, and the major program successes. By describing our experience, we aim to assist other clinical teams in achieving sustainable PGx implementation in their health systems.

The Potential Impact of Preemptive Pharmacogenetic Genotyping in the Neonatal Intensive Care Unit

OBJECTIVE

To evaluate the use of drugs with pharmacogenomic (PGx) guidelines from the Clinical Pharmacogenetics Implementation Consortium in early childhood.

STUDY DESIGN

A retrospective observational study of patients admitted to the neonatal intensive care (NICU) between 2005 and 2018 with at least 1 subsequent hospitalization at or after 5 years of age was performed to determine PGx drug exposure. Data regarding hospitalizations, drug exposures, gestational age, birth weight, and congenital anomalies and/or a primary genetic diagnosis were collected. Incidence of PGx drug and drug class exposures was determined and patient specific factors predictive of exposure were investigated.

RESULTS

During the study, 19 195 patients received NICU care and 4196 (22%) met study inclusion; 67% received 1-2, 28% 3-4, and 5% 5 or more PGx-drugs in early childhood. Preterm gestation, low birth weight (<2500 g), and the presence of any congenital anomalies and/or a primary genetic diagnosis were statistically significant predictors of Clinical Pharmacogenetics Implementation Consortium drug exposures (P < .01, P < .01, P < .01, respectively).

CONCLUSIONS

Preemptive PGx testing in patients in the NICU could have a significant impact on medical management during the NICU stay and throughout early childhood.

Implementation of Clinical Cytochrome P450 3A Genotyping for Tacrolimus Dosing in a Large Kidney Transplant Program

Tacrolimus is a calcineurin inhibitor with a narrow therapeutic range and is metabolized by cytochrome P450 (CYP) isoenzymes CYP3A4 and CYP3A5. The Clinical Pharmacogenetic Implementation Consortium published evidence-based guidelines for CYP3A5 normal/intermediate metabolizers prescribed tacrolimus, yet few transplant centers have implemented routine testing. The objective of this study was to implement preemptive CYP3A genotyping into clinical practice in a large kidney transplant program and to evaluate workflow feasibility, potential clinical benefit, and reimbursement to identify barriers and determine sustainability. Preemptive pharmacogenetic testing for CYP3A5 and CYP3A4 was implemented in all patients listed for a kidney transplant as part of standard clinical care. Genotyping was performed at the listing appointment, results were reported as discrete data in the electronic medical record, and education and clinical decision support alerts were developed to provide pharmacogenetic-recommended tacrolimus dosing. During this initial phase, all patients were administered standard tacrolimus dosing, and clinical and reimbursement outcomes were collected. Greater than 99.5% of genotyping claims were reimbursed by third-party payers. CYP3A5 normal/intermediate metabolizers had significantly fewer tacrolimus trough concentrations within the target range and a significantly longer time to their first therapeutic trough compared to poor metabolizers. The challenge of tacrolimus dosing is magnified in the African American population. The US Food and Drug Administration drug label recommends increased starting doses in African ancestry, yet only ≈66% of African Americans in our cohort were normal/intermediate metabolizers who required higher doses. Routine CYP3A5 genotyping may overcome this issue by using genotype over race as a more accurate predictor of drug response.

Cytochrome P450 Oxidoreductase (POR) Associated with Severe Paclitaxel-Induced Peripheral Neuropathy in Patients of European Ancestry from ECOG-ACRIN E5103

PURPOSE

Paclitaxel is a widely used anticancer therapeutic. Peripheral neuropathy is the dose-limiting toxicity and negatively impacts quality of life. Rare germline gene markers were evaluated for predicting severe taxane-induced peripheral neuropathy (TIPN) in the patients of European ancestry. In addition, the impact of Cytochrome P450 (CYP) 2C8, CYP3A4, and CYP3A5 metabolizer status on likelihood of severe TIPN was also assessed.

EXPERIMENTAL DESIGN

Whole-exome sequencing analyses were performed in 340 patients of European ancestry who received a standard dose and schedule of paclitaxel in the adjuvant, randomized phase III breast cancer trial, E5103. Patients who experienced grade 3-4 (n = 168) TIPN were compared to controls (n = 172) who did not experience TIPN. For the analyses, rare variants with a minor allele frequency ≤ 3% and predicted to be deleterious by protein prediction programs were retained. A gene-based, case-control analysis using SKAT was performed to identify genes that harbored an imbalance of deleterious variants associated with increased risk of severe TIPN. CYP star alleles for CYP2C8, CYP3A4, and CYP3A5 were called. An additive logistic regression model was performed to test the association of CYP2C8, CYP3A4, and CYP3A5 metabolizer status with severe TIPN.

RESULTS

Cytochrome P450 oxidoreductase (POR) was significantly associated with severe TIPN (P value = 1.8 ×10-6). Six variants were predicted to be deleterious in POR. There were no associations between CYP2C8, CYP3A4, or CYP3A5 metabolizer status with severe TIPN.

CONCLUSIONS

Rare variants in POR predict an increased risk of severe TIPN in patients of European ancestry who receive paclitaxel.

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.

An efficient genotyper and star-allele caller for pharmacogenomics

High-throughput sequencing provides sufficient means for determining genotypes of clinically important pharmacogenes that can be used to tailor medical decisions to individual patients. However, pharmacogene genotyping, also known as star-allele calling, is a challenging problem that requires accurate copy number calling, structural variation identification, variant calling, and phasing within each pharmacogene copy present in the sample. Here we introduce Aldy 4, a fast and efficient tool for genotyping pharmacogenes that uses combinatorial optimization for accurate star-allele calling across different sequencing technologies. Aldy 4 adds support for long reads and uses a novel phasing model and improved copy number and variant calling models. We compare Aldy 4 against the current state-of-the-art star-allele callers on a large and diverse set of samples and genes sequenced by various sequencing technologies, such as whole-genome and targeted Illumina sequencing, barcoded 10x Genomics, and Pacific Biosciences (PacBio) HiFi. We show that Aldy 4 is the most accurate star-allele caller with near-perfect accuracy in all evaluated contexts, and hope that Aldy remains an invaluable tool in the clinical toolbox even with the advent of long-read sequencing technologies.

Development of an extensive workflow for comprehensive clinical pharmacogenomic profiling: lessons from a pilot study on 100 whole exome sequencing data

This pilot study is aimed at implementing an approach for comprehensive clinical pharmacogenomics (PGx) profiling. Fifty patients with cardiovascular diseases and 50 healthy individuals underwent whole-exome sequencing. Data on 1800 PGx genes were extracted and analyzed through deep filtration separately. Theoretical drug induced phenoconversion was assessed for the patients, using sequence2script. In total, 4539 rare variants (including 115 damaging non-synonymous) were identified. Four publicly available PGx bioinformatics algorithms to assign PGx haplotypes were applied to nine selected very important pharmacogenes (VIP) and revealed a 45-70% concordance rate. To ensure availability of the results at point-of-care, actionable variants were stored in a web-hosted database and PGx-cards were developed for quick access and handed to the study subjects. While a comprehensive clinical PGx profile could be successfully extracted from WES data, available tools to interpret these data demonstrated inconsistencies that complicate clinical application.