Stargazer: a software tool for calling star alleles from next-generation sequencing data using CYP2D6 as a model.

PURPOSE: Genotyping CYP2D6 is important for precision drug therapy because the enzyme it encodes metabolizes approximately 25% of drugs, and its activity varies considerably among individuals. Genotype analysis of CYP2D6 is challenging due to its highly polymorphic nature. Over 100 haplotypes (star alleles) have been defined for CYP2D6, some involving a gene conversion with its nearby nonfunctional but highly homologous paralog CYP2D7. We present Stargazer, a new bioinformatics tool that uses next-generation sequencing (NGS) data to call star alleles for CYP2D6 ( https://stargazer.gs.washington.edu/stargazerweb/ ). Stargazer is currently being extended for other pharmacogenes. METHODS: Stargazer identifies star alleles from NGS data by detecting single nucleotide variants, insertion-deletion variants, and structural variants. Stargazer detects structural variation, including gene deletions, duplications, and conversions, by calculating paralog-specific copy numbers from read depths. RESULTS: We applied Stargazer to the NGS data of 32 ethnically diverse HapMap trios that were genotyped by TaqMan assays, long-range polymerase chain reaction, quantitative multiplex polymerase chain reaction, high-resolution melting analysis, and/or Sanger sequencing. CYP2D6 genotyping by Stargazer was 99.0% concordant with the data obtained by these methods, and showed that 28.1% of the samples had structural variation including CYP2D6/CYP2D7 hybrids. CONCLUSION: Accurate genotyping of pharmacogenes with NGS and subsequent allele calling with Stargazer will aid the implementation of precision drug therapy.

Pharmacogenetic analysis of structural variation in the 1000 genomes project using whole genome sequences.

While significant strides have been made in understanding pharmacogenetics (PGx) and gene-drug interactions, there remains limited characterization of population-level PGx variation. This study aims to comprehensively profile global star alleles (haplotype patterns) and phenotype frequencies in 58 pharmacogenes associated with drug absorption, distribution, metabolism, and excretion. PyPGx, a star-allele calling tool, was employed to identify star alleles within high-coverage whole genome sequencing (WGS) data from the 1000 Genomes Project (N = 2504; 26 global populations). This process involved detecting structural variants (SVs), such as gene deletions, duplications, hybrids, as well as single nucleotide variants and insertion-deletion variants. The majority of our PyPGx calls for star alleles and phenotype frequencies aligned with the Pharmacogenomics Knowledge Base, although notable population-specific frequencies differed at least twofold. Validation efforts confirmed known SVs while uncovering several novel SVs currently undefined as star alleles. Additionally, we identified 210 small nucleotide variants associated with severe functional consequences that are not defined as star alleles. The study serves as a valuable resource, providing updated population-level star allele and phenotype frequencies while incorporating SVs. It also highlights the burgeoning potential of cost-effective WGS for PGx genotyping, offering invaluable insights to improve tailored drug therapies across diverse populations.

Characterizing the combined effects of cytochrome P450 missense variation within star allele definitions.

Background: Cytochrome P450 (CYP) genetic variation largely impacts drug response. However, many CYP star alleles (haplotypes) lack functional annotation, impeding our understanding of drug metabolism mechanisms. We aimed to investigate the impact of missense variant combinations on CYP protein structures. Methods: Normal mode analysis was conducted on 261 missense variants within 91 CYP haplotypes. CYP2D6*2 and CYP2D6*17 were prioritized for molecular dynamics simulation. Results: Normal mode analysis and molecular dynamics highlight the effects of known CYP missense variants on protein stability and conformational dynamics. Missense variants within haplotypes may have intermodulating effects on protein structure and function. Conclusion: This study highlights the utility of multiscale modeling in interpreting CYP missense variants and particularly their combinations within various star alleles.

Characterization of CYP2C19 pharmacogenetic variation in African populations and comparison with other global populations.

Background: CYP2C19 is important in the metabolism of clopidogrel and several antidepressants. This study aimed to characterize the distribution of CYP2C19 star alleles (haplotypes) across diverse African populations compared with global populations. Methods: CYP2C19 star alleles and diplotypes were called from high coverage genomes using the StellarPGx pipeline. Results: CYP2C19*1 (51%), *2 (17%) and *17 (22%) were the most common star alleles across African populations in this study. It was observed that 3% of African participants had potentially novel CYP2C19 haplotypes. Conclusion: This study supports the necessity for CYP2C19 pharmacogenetic testing in African and global clinical settings, as well as the importance of comprehensive star allele characterization in the African context.

Characterization of POR haplotype distribution in African populations and comparison with other global populations.

Background & aim: POR is an enzyme that mediates electron transfer to enable the drug-metabolizing activity of CYP450 proteins. However, POR has been understudied in pharmacogenomics despite this vital role. This study aimed to characterize the genetic variation in POR across African populations and to compare the star allele (haplotype) distribution with that in other global populations. Materials & methods: POR star alleles were called from whole-genome sequencing data using the StellarPGx pipeline. Results: In addition to the common POR*1 and *28 (defined by rs1057868), five novel rare haplotypes were computationally inferred. No significant frequency differences were observed among the majority of African populations. However, POR*28 was observed at a higher frequency in individuals of non-African ancestry. Conclusion: This study highlights the distribution of POR alleles in Africa and across global populations with a view toward informing future precision medicine implementation.

PharmaKU: A Web-Based Tool Aimed at Improving Outreach and Clinical Utility of Pharmacogenomics.

With the tremendous advancements in genome sequencing technology in the field of pharmacogenomics, data have to be made accessible to be more efficiently utilized by broader clinical disciplines. Physicians who require the drug-genome interactome information, have been challenged by the complicated pharmacogenomic star-based classification system. We present here an end-to-end web-based pharmacogenomics tool, PharmaKU, which has a comprehensive easy-to-use interface. PharmaKU can help to overcome several hurdles posed by previous pharmacogenomics tools, including input in hg38 format only, while hg19/GRCh37 is now the most popular reference genome assembly among clinicians and geneticists, as well as the lack of clinical recommendations and other pertinent dosage-related information. This tool extracts genetic variants from nine well-annotated pharmacogenes (for which diplotype to phenotype information is available) from whole genome variant files and uses Stargazer software to assign diplotypes and apply prescribing recommendations from pharmacogenomic resources. The tool is wrapped with a user-friendly web interface, which allows for choosing hg19 or hg38 as the reference genome version and reports results as a comprehensive PDF document. PharmaKU is anticipated to enable bench to bedside implementation of pharmacogenomics knowledge by bringing precision medicine closer to a clinical reality.