Pharmacogenomics¶
Pharmacogenomics (PGx) examines how genetic variation affects an individual's response to medications. AIVA's PGx analysis automatically assigns star alleles to 88 pharmacogenes, predicts metabolizer phenotypes, and generates clinical drug recommendations based on CPIC guidelines.
PGx results help clinicians identify patients who may need dose adjustments, alternative medications, or enhanced monitoring based on their genetic profile.
How it works¶
PGx analysis runs automatically as part of the secondary analysis pipeline when you upload FASTQ files. The pipeline proceeds through the following stages:
- Variant calling: Variants are called from the BAM file across all 88 pharmacogenes.
- Star allele assignment: Each gene receives a diplotype (a pair of star alleles, e.g.,
*1/*4) based on detected SNVs. - Phenotype prediction: Metabolizer phenotypes are predicted from the diplotype (e.g., Poor Metabolizer, Normal Metabolizer).
- Activity scoring: An activity score is calculated by summing the functional scores of each allele.
- Drug recommendations: CPIC guidelines are consulted to generate prescribing recommendations for each gene's phenotype.
- Multi-gene enrichment: Drugs that require phenotypes from two genes (e.g., fluvastatin requires both CYP2C9 and SLCO1B1) are identified and enriched with combined recommendations.
Input requirements
PGx analysis requires a BAM file generated by the Parabricks pipeline. It supports both short-read (Illumina) and long-read (PacBio, ONT) sequencing data aligned to GRCh37 or GRCh38.
Supported genes¶
AIVA analyzes 88 pharmacogenes. These genes are organized into the following categories:
CYP450 enzymes¶
Cytochrome P450 enzymes are responsible for metabolizing the majority of clinically used drugs.
| Gene | Common drug associations |
|---|---|
| CYP1A1, CYP1A2, CYP1B1 | Caffeine, theophylline, clozapine |
| CYP2A6, CYP2A13 | Nicotine, tegafur |
| CYP2B6 | Efavirenz, methadone, bupropion |
| CYP2C8 | Paclitaxel, repaglinide |
| CYP2C9 | Warfarin, phenytoin, NSAIDs |
| CYP2C19 | Clopidogrel, proton pump inhibitors, voriconazole |
| CYP2D6 | Codeine, tamoxifen, antidepressants, antipsychotics |
| CYP2E1 | Acetaminophen, isoniazid |
| CYP2F1, CYP2J2, CYP2R1, CYP2S1, CYP2W1 | Various substrates |
| CYP3A4, CYP3A5, CYP3A7, CYP3A43 | Tacrolimus, cyclosporine, statins |
| CYP4A11, CYP4A22, CYP4B1, CYP4F2 | Warfarin (CYP4F2), arachidonic acid metabolism |
| CYP17A1, CYP19A1, CYP26A1 | Steroid hormone metabolism |
Phase II enzymes¶
These enzymes catalyze conjugation reactions that facilitate drug elimination.
| Gene | Common drug associations |
|---|---|
| GSTM1, GSTP1, GSTT1 | Busulfan, platinum agents |
| NAT1, NAT2 | Isoniazid, hydralazine, sulfonamides |
| SULT1A1 | Tamoxifen, minoxidil |
| TPMT | Azathioprine, mercaptopurine, thioguanine |
| UGT1A1, UGT1A4, UGT1A6 | Irinotecan (UGT1A1), lamotrigine, acetaminophen |
| UGT2B7, UGT2B15, UGT2B17 | Morphine, lorazepam, testosterone |
Transporters¶
Drug transporters affect absorption, distribution, and elimination of medications.
| Gene | Common drug associations |
|---|---|
| ABCB1 (MDR1), ABCG2 (BCRP) | Digoxin, statins, antiretrovirals |
| SLC6A4 | SSRIs (serotonin transporter) |
| SLC15A2, SLC22A2, SLC28A3, SLC47A2 | Various drug substrates |
| SLCO1B1, SLCO1B3, SLCO2B1 | Statins, methotrexate |
Other pharmacogenes¶
| Gene | Relevance |
|---|---|
| DPYD | Fluoropyrimidines (5-FU, capecitabine) |
| NUDT15 | Thiopurines (azathioprine, mercaptopurine) |
| VKORC1 | Warfarin sensitivity |
| G6PD | Rasburicase, primaquine, dapsone |
| F2, F5 | Thrombosis risk, anticoagulant therapy |
| CACNA1S, RYR1 | Malignant hyperthermia susceptibility |
| IFNL3 | Interferon-based hepatitis C therapy |
| CFTR | Cystic fibrosis modulators |
| MTHFR | Methotrexate, folate metabolism |
| OPRM1, OPRK1 | Opioid response |
| COMT | Catecholamine metabolism, pain sensitivity |
| APOE | Statin response, Alzheimer's risk |
| MT-RNR1 | Aminoglycoside-induced hearing loss |
| DRD2, ANKK1, HTR1A, HTR2A | Psychiatric medications |
| BDNF, GRIK1, GRIK4, GRIN2B | Neuropsychiatric drug response |
| ADRA2A, ADRB2, DBH | Adrenergic drug response |
| ACYP2 | Cisplatin-induced hearing loss |
| ATM | DNA damage response |
| BCHE | Succinylcholine, mivacurium metabolism |
| ITGB3, ITPA | Platelet function, ribavirin response |
| POR | CYP enzyme electron transfer |
| PTGIS, TBXAS1 | Prostaglandin/thromboxane pathways |
| RARG | Anthracycline-induced cardiotoxicity |
| XPC | DNA repair, drug sensitivity |
Understanding results¶
Each gene in your PGx results includes the following fields:
Diplotype¶
The diplotype represents the two star alleles identified for a gene. Star alleles are standardized nomenclature for common haplotypes. For example, CYP2D6 *1/*4 indicates one normal-function allele (*1) and one no-function allele (*4).
Phenotype¶
The predicted metabolizer status based on the diplotype:
| Phenotype | Meaning |
|---|---|
| Ultrarapid Metabolizer | Increased enzyme activity; may need higher doses or alternative drugs |
| Rapid Metabolizer | Above-normal enzyme activity |
| Normal Metabolizer | Typical enzyme activity; standard dosing expected |
| Intermediate Metabolizer | Reduced enzyme activity; may need lower doses |
| Poor Metabolizer | Little to no enzyme activity; significant dose reduction or alternative drug needed |
Note
Not all genes have metabolizer phenotype predictions. Some genes have other phenotype classifications (e.g., warfarin sensitivity for VKORC1, malignant hyperthermia susceptibility for RYR1).
Activity score¶
A numeric value representing total enzyme activity, calculated by summing the functional scores of each allele. Higher scores indicate greater enzyme activity. Activity scores provide finer granularity than phenotype categories alone.
Allele functions¶
Each allele receives a functional classification:
- Normal Function: Standard activity
- Increased Function: Above-normal activity
- Decreased Function: Reduced activity
- No Function: Absent activity
- Uncertain Function: Insufficient evidence to classify
CNV detection¶
For genes with known structural variation (e.g., CYP2D6, GSTM1, GSTT1, UGT2B17), AIVA performs a separate copy number variation (CNV) analysis. CNV results are informational only and do not override the SNV-based diplotype. Possible CNV calls include Normal, WholeDel1, WholeDup1, and Tandem duplications.
Review flag
When the SNV-based diplotype is indeterminate but a CNV is detected, the result is flagged for clinical review. This may indicate a genuine structural variant requiring manual interpretation.
Drug recommendations¶
AIVA generates prescribing recommendations based on CPIC (Clinical Pharmacogenetics Implementation Consortium) guidelines.
Single-gene recommendations¶
Most drug recommendations are based on a single gene's phenotype. For example:
- CYP2D6 Poor Metabolizer + codeine: "Use alternative analgesic; increased risk of reduced efficacy"
- CYP2C19 Ultrarapid Metabolizer + clopidogrel: "Standard dosing; expected normal response"
Multi-gene recommendations¶
Some drugs require phenotype information from two genes for accurate guidance. AIVA automatically detects these cases and provides combined recommendations. For example:
- Fluvastatin: Requires both CYP2C9 and SLCO1B1 phenotypes
- Warfarin: Benefits from combined CYP2C9 and VKORC1 information
Multi-gene recommendations appear in the results for both contributing genes and are labeled accordingly.
Tip
Drug recommendations are guidelines based on population-level evidence. Clinical decisions should consider the full patient context, including other medications, organ function, and clinical presentation.