Next Article in Journal
A COL7A1 Variant in a Litter of Neonatal Basset Hounds with Dystrophic Epidermolysis Bullosa
Next Article in Special Issue
Distinct Effects of Inflammation on Cytochrome P450 Regulation and Drug Metabolism: Lessons from Experimental Models and a Potential Role for Pharmacogenetics
Previous Article in Journal
Development and Cell Cycle Activity of the Root Apical Meristem in the Fern Ceratopteris richardii
Previous Article in Special Issue
Using Personal Genomic Data within Primary Care: A Bioinformatics Approach to Pharmacogenomics
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 12
10.3390/genes11121456
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Technologies for Pharmacogenomics: A Review

by
Maaike van der Lee
1,2,
Marjolein Kriek
3,
Henk-Jan Guchelaar
1,2 and
Jesse J. Swen
1,2,*
1
Department of Clinical Pharmacy and Toxicology, Leiden University Medical Center, 2333ZA Leiden, The Netherlands
2
Leiden Network of Personalized Therapeutics, 2333ZA Leiden, The Netherlands
3
Department of Clinical Genetics, Leiden University Medical Center, 2333ZA Leiden, The Netherlands
*
Author to whom correspondence should be addressed.
Genes 2020, 11(12), 1456; https://doi.org/10.3390/genes11121456
Submission received: 13 November 2020 / Revised: 30 November 2020 / Accepted: 2 December 2020 / Published: 4 December 2020
(This article belongs to the Special Issue Pharmacogenomic Determinants of Interindividual Drug Response Variability: From Discovery to Implementation)
Browse Figure
Versions Notes

Abstract

:
The continuous development of new genotyping technologies requires awareness of their potential advantages and limitations concerning utility for pharmacogenomics (PGx). In this review, we provide an overview of technologies that can be applied in PGx research and clinical practice. Most commonly used are single nucleotide variant (SNV) panels which contain a pre-selected panel of genetic variants. SNV panels offer a short turnaround time and straightforward interpretation, making them suitable for clinical practice. However, they are limited in their ability to assess rare and structural variants. Next-generation sequencing (NGS) and long-read sequencing are promising technologies for the field of PGx research. Both NGS and long-read sequencing often provide more data and more options with regard to deciphering structural and rare variants compared to SNV panels—in particular, in regard to the number of variants that can be identified, as well as the option for haplotype phasing. Nonetheless, while useful for research, not all sequencing data can be applied to clinical practice yet. Ultimately, selecting the right technology is not a matter of fact but a matter of choosing the right technique for the right problem.

1. Introduction

The field of pharmacogenomics (PGx) is developing rapidly. The first PGx dose recommendations for antidepressant and psychiatric drugs were published in 2001, even before the first human genome was sequenced [1]. An increase in available evidence and the ambition to implement PGx in clinical practice has led to the need for more comprehensive dosing guidelines and genotyping strategies. In 2005, the Dutch Pharmacogenetics Working Group (DPWG) was formed to develop evidence-based PGx guidelines [2]. In 2011, the Clinical Pharmacogenomics Implementation Consortium (CPIC) was founded [3]. Currently, CPIC and the DPWG combined have issued PGx dose recommendations covering more than 50 drugs and 21 genes (Table 1) [4,5].
In the last 20 years, there has not only been much progress in development of PGx guidelines, there has also been a significant technological advancement and rise of new technologies for assessing genetic variants. At the time of development of the first PGx guideline, only Sanger-based sequencing techniques and SNV (single nucleotide variant) arrays were available as methods for variant identification. To date, SNV panel testing remains the most commonly used technology in clinical practice. However, while it is efficient and comes at low costs, SNV panels cannot detect all important genetic variation such as rare and structural variants. Currently, multiple high throughput whole genome sequencing techniques are available, yielding an abundance of genetic information at a fraction of the costs of 20 years ago [6,7,8,9,10]. Nevertheless, these approaches are not yet routinely used in other clinical fields, despite their potential [11,12].
In this paper, we review the application of genotyping technologies for PGx. We first discuss the use of SNV panels, which are the most commonly used approach for clinical PGx. Next, we discuss the potential of next generation sequencing (NGS) and long-read sequencing and their current use in PGx. We review the challenges in PGx, both for clinical as well as research purposes, and the way PGx technologies can help in solving them. We mainly focus on germline variants and their role in PGx. Nonetheless, the outlined principles hold true for somatic mutations [13,14,15]. Additionally, implementation and adoption of PGx in clinical practice is outside the scope of this review and has been extensively discussed in previous publications [16,17,18].

2. SNV Panels: Current Clinical Practice

SNV panel testing is the most commonly used technology in PGx practice, either through commercially available micro-array platforms or with custom arrays. The arrays typically contain a preselected set of SNVs, which, depending on the array and platform, can range from a few variants in a single gene to thousands of variants genome wide. Commercially available PGx arrays typically contain variants that are linked to drug response in PGx guidelines or on PharmGKB [19]. The evidence underlying the selected variants can vary, from small arrays containing only the most strongly associated variants, to very large arrays containing all variants potentially or theoretically associated with drug response—for example, including all known drug-related genes. Almost all available arrays use PCR, sequencing by synthesis and nanospheres or beads, combined with a form of fluorescence or chemiluminescence detection to identify which variant is present at the site of interest [20,21,22]. Another technology is the use of mass spectrometry, relying on differences in mass between wildtype and mutant nucleotides [23]. Detailed descriptions of these techniques have been described previously [3,24,25,26,27,28]. The pre-selection of variants and the relatively low amount of data to process allow for a quick result at low costs (Table 2).

2.1. Commercial Arrays

There are many arrays available that can be used for PGx; a full overview of these arrays is beyond the scope of this review. Two of the smaller commercial arrays are the VeraCode ADME core panel (Illumina Inc. San Diego, CA, USA) and the VeriDose core panel (Agena BioScience, San Diego, CA, USA). The VeraCode consists of 184 variants in 34 pharmacogenes [20] and the Veridose contains 68 variants in 20 genes and 5 copy number targets for CYP2D6 [30]. The ADME core variant list is based on an expert gene panel and contains the most biologically relevant variants within these genes [20,31]. For the VeriDose, genes with a known clinical impact and their common clinically actionable variants are selected. Additionally, it is possible to expand the panel if so desired [30]. Both panels provide sufficient coverage for clinical PGx by covering the most common variants in actionable pharmacogenes. More extensive panels are the pharmacofocus with (2000 variants in 150 genes including CNVs (Copy Number Variants) (ThermoFisher Scientific, Waltham, MA, USA) [32] and the pharmacoscan with 4627 variants in 1191 genes (ThermoFisher Scientific) [22]. The latter contain nearly all variants from the DMET and Illumina ADME core panel, in addition to all genes and variants with clinical annotations in CPIC and PharmGKB (Pharmacogenomics knowledge base), HLA genes and sample ID and tracking markers. These types of arrays are widely used in clinical PGx implementation studies. For example, the DMET array is used in the PG4KDS study from St. Jude’s children’s research hospital [3] and the VeraCode ADME core panel is used in the PREDICT study [33]. Both these studies use only a subset of the variants available on the panel for the clinical implementation part of the studies [3,33]. Only the variants in the genes of interest and the variants with sufficient data are reported for clinical practice. The remainder of the genetic data is, with informed consent permission, stored and can be used for research purposes or for later clinical use. For a full overview of studies using PGx panel approaches, we refer to previous publications providing such an overview [11,34].

2.2. Custom Arrays

Several of the commercial arrays contain a high number of variants, making a fast turnaround time and interpretation challenging. Additionally, these arrays will include variants which may not be of direct interest in a clinical setting as panels often include all known PGx variants regardless of the level of evidence supporting their clinical utility. This has driven many institutions to develop their own custom clinical array with a more focused set of clinically actionable variants. Typically, these custom arrays are performed using single gene-based testing covering only the variants with direct clinical applicability yet limiting broad applicability. To be able to test a broad set of variants while maintaining a rapid turnaround time, companies have developed customizable arrays. One example frequently used for clinical PGx is the OpenArray (Thermofisher scientific). This array can detect between 12 and 240 variants using standard TaqMan technology [35]. Selected TaqMan assays are spotted on a chip based on the clients need. Each chip contains wells for the patient samples; in each of these wells, through-holes are present which contain the assays desired. This format allows for a thorough characterization of a few genes as well as for a broader approach focused on a panel of common variants in multiple genes. This array has been used in clinical implementation studies as well. For example, in the INGENIOUS study it was used to interrogate a panel of 43 variants in 14 genes [36]. A similar approach was used by the Ubiquitous pharmacogenomics (U-PGx) consortium. The U-PGx consortium’s initiated the PREPARE study aimed at collecting evidence of the clinical utility of a pre-emptive PGx panel consisting of 58 SNVs in 14 pharmacogenes [34]. The panel covers the most common variants in all actionable genes included in the DPWG guidelines [37] and is analyzed with KASP technology using the SNPline (LGC) [38].

2.3. Array Developments

The above-mentioned commercial and custom arrays are developed specifically for PGx. However, there are also multiple arrays with genome-wide coverage available. An added benefit of genome wide arrays is that they also allow for GWAS (Genome Wide Association Study) analysis in addition to providing PGx information. Examples of such arrays are the Illumina GSA (Global Screening Array) which contains over 600,000 SNVs genome-wide, including 17,750 PGx markers [39] and the Axiom arrays (Thermofisher scientific), which contain genome-wide coverage specifically for a certain population [40]. Nonetheless, these arrays often miss dense coverage of PGx regions and not all critical SNVs are available on the array. For example, in the case of the GSA v3.0, the SNVs for CYP2D6*4 and CYP2C19*9 and target CNV testing are not included. This is particularly concerning for the CYP2D6 gene. The CYP2D6 enzyme is responsible for 25–30% of commonly prescribed drugs, making sufficient coverage on variants on the CYP2D6 gene clinical important [41]. The CYP2D6*4 allele is the most frequent null-allele in Caucasians, with the key SNV (rs3892097; NC_000022.11:g.42128945C>T) occurring in 19% of the European (non-Finnish) population [42,43]. Additionally, the GSA v3.0 does not contain probes for direct CNV detection.

3. Next Generation Sequencing

3.1. Next Generation Sequencing Technologies

NGS technologies are not yet routinely applied in clinical PGx. However, they are often used in PGx research and disease genetics. While SNV panels only cover a limited set of selected variants, sequencing data cover the full exome or genome. Technical details of NGS technologies have been extensively reviewed elsewhere [44,45,46]. In short, NGS technologies are capable of sequencing reads of 100–200 bp in a high throughput manner, allowing for the sequencing of a full genome in a matter of hours. These reads are aligned to the reference genome and variants are identified based on deviations from the reference.
NGS applications can be roughly categorized into three approaches. First, whole exome sequencing (WES) focusing on sequencing the coding regions of the genome, covering approximately 1–2% of the entire genome. Secondly, whole genome sequencing (WGS) which is aimed at sequencing the entire genome, both coding and non-coding regions. Lastly, targeted sequencing of a region or panel of genes of interest [44,45,46]. While NGS can be performed at relatively low costs, the large amount of data makes processing more challenging (Table 2).

3.2. Use of NGS for Pharmacogenomics

While NGS has become the standard for clinical diagnostics and in research, it is yet to be widely adopted for clinical PGx. Nonetheless, in a research setting NGS has been used for PGx for several years. Multiple studies have been conducted investigating the accuracy of NGS technologies in PGx as well as the applicability of an NGS sequencing panel or of the repurposing of clinical NGS data for PGx [9,10,47,48,49,50,51]. Yang et al. performed a three-way analysis with the DMET, WES and WGS, to investigate the concordance between PGx genotyping calls based on these different technologies. They showed a 94% concordance between the DMET and WES, and a 96% concordance between the DMET and WGS [47]. Similar results were reported by other groups, all of which report the superior results obtained from sequencing compared to orthogonal testing [49,50,51,52,53]. The difference in concordance between WES and the DMET array (94%), and WGS and the DMET array (96%) can be explained by the genomic coverage of each approach. WES only covers the exons and can therefore, by definition, not cover all relevant variants if they are located in the intronic or intergenic regions. WGS, on the other hand, also covers intronic regions leading to an expanded coverage. Nonetheless, intronic variants are of clinical importance in PGx. For example, one of the key CYP2C19*17 variants (rs12248560; NC_000010.11:g.94761900C>T) is located upstream of the CYP2C19 locus. Other examples are CYP3A5*3 and *5 as well as CYP2D6*4 and *41 [54]. A targeted sequencing approach can combine the lower costs of WES with the advances of WGS data. This type of panel only captures genes of interest, both intronic and exonic regions. This results in lower costs while maintaining the accuracy and abundance of data of WGS. One such approach is the PGRNseq panel [52]. This panel is based on full-gene sequencing of a panel of 84 pharmacogenes using NGS, it has also been used in clinical implementation studies showing promising results [49,55]. Another approach is the PGxSeq panel described by Gulilat et al., which covers 100 pharmacogenes [56].

3.3. Repurposing of Clinical Genetics Data

The abovementioned approaches are aimed at generating novel sequencing data with the goal of providing PGx results, whether it be panel, WES or WGS-based. However, in clinical diagnostics, the use of NGS is already standard practice in many centers leading to vast amounts of sequencing data. Several groups have investigated the feasibility of repurposing these data for PGx, by extracting a panel of evidence-based PGx variants from the data and translating this into a clinically applicable result [9,10]. The same type of analysis is performed for large populations studies, such as the Estonian biobank [8], the University of Colorado [57], SWEDEGENE [7] and the AllofUs initiative [58,59]. Unfortunately, the utility of the repurposing of data is dependent on the capture panel used in the original sequencing. This is particularly a problem in the use of WES data, as mentioned above, several important PGx variants are located in intronic regions which are not included in WES capture kits [54]. Nonetheless, even with these limitations, high percentages (>85%) of individuals with actionable phenotypes are identified [9,10,60], but one should bear in mind that important variants are missing.

4. Long-Read Sequencing

4.1. Long-Read Sequencing

Long-read sequencing technologies have emerged in the playing field and are slowly gaining ground over the short-read approaches in the field of research [46,61]. Both Pacific Bioscience (PacBio) technology as well as Oxford Nanopore Technologies (ONT) are becoming an integrated part of genetic approaches [12]. PacBio uses SMRT (single molecule real-time)-sequencing to be able to sequence reads up to 45 kB. SMRT cells make use of microwells, each of which contains one single strand of DNA which is then sequenced by assembly and recorded in real time [62]. Oxford Nanopore Technologies uses nanopores through which the DNA strand is pulled, the disruption in the current is specific to a codon, allowing for the full assembly of the DNA sequence [63]. By correcting for the randomly distributed errors in single cell sequencing, the consensus reads can obtain very high accuracy [12,62,64]. While an abundance of data can be generated by long-read sequencing, the processing is significantly more intensive compared to SNV panels and NGS (Table 2).

4.2. Long-Read Sequencing for PGx

Long-read sequencing has been shown to be capable of solving complex loci genome wide [12,64]. Current disease diagnostics already use long-read sequencing for complex genetic diseases such as the ATXN10 repeats in Parkinson’s disease [65] and tandem repeats in the FMR1 gene associated with Fragile X syndrome [66]. Nonetheless, only a few long-read sequencing studies for pharmacogenes have been conducted. The most thoroughly investigated complex locus in PGx is the CYP2D6 gene, which contains both SNVs and SVs. (Structural Variants) It has been shown that with long-read sequencing, the CYP2D6 locus (~6.6 kbp) can be sequenced in one full read and be fully resolved into phased haplotypes, including structural variants [67,68]. The same has been observed for the notoriously complex HLA genes [12,64]. Long-read sequencing in PGx is currently, to our knowledge, limited to single gene studies and no large-scale studies applying long-read sequencing for clinical PGx to a panel of genes have been conducted.

5. Challenges

5.1. Drug Metabolizer Phenotype Inference

Prior to application in clinical practice, SNV data are translated into predicted drug metabolizer phenotypes. Many of the arrays mentioned above are developed based on the SNVs which are present in the genotype nomenclature. For CYP enzymes, haplotypes are named with the star (*) nomenclature [42,69]. All variants making up a *-allele are described by the Pharmacogene variation Consortium (PharmVar, https://www.pharmvar.org). The combination of the two *-alleles present in an individual is subsequently translated into a predicted drug metabolizer phenotype. For most pharmacogenes, there are four metabolizer phenotypes: normal metabolizers (NMs) displaying full protein function, intermediate metabolizers (IMs) associated with decreased protein function, poor metabolizers (PMs) indicating absence of protein function and ultra-rapid metabolizers (UMs), which are associated with increased enzyme function.
Depending on the number of variants and the presence of the variants in translation guidelines, the interpretation is relatively straightforward. However, if there are many variants leading to *-haplotypes of unknown function present on the array, the interpretation is challenging. Furthermore, PharmVar describes the *-haplotypes extensively, defining large numbers haplotypes and sub haplotypes which are constantly increasing. In theory, all variants defined by PharmVar at a certain point in time could be included on an SNV panel. Nonetheless, with over 2000 variants known in CYP2D6 alone, this would quickly grow to a very large panel which is difficult to apply and interpret in clinical practice [42]. As there is no standardization in regard to the variants which need to be tested in a clinical setting, every array contains its own set of variants. This can lead to differences in regard to assigned haplotypes when testing the same DNA in different laboratories [49]. For example, results from the pharmacoscan (4627 variants) will be much more extensive and detailed compared to the VeriDose Core panel (68 variants) [22,30]. The pharmacoscan analysis might also result in more variants and haplotypes of unknown effect while, on the other hand, the VeriDose core panel might miss variants with a known effect.
The use of sequencing data enables the inclusion of almost every known variant in the *-haplotype assignments and subsequently into predicted phenotypes. However, with this abundance of variant data comes an increased difficulty in haplotype assignments. Several tools have been developed to assign *-allele haplotypes based on sequencing data, incorporating all variants in the assignment. Yet, as Caspar et al. have shown, these tools are not performing perfectly with errors in the assigned haplotypes compared to the consensus genotype from GeT-RM [60]. An error is here defined as a haplotype assignment which differs from the consensus. The best performing tool was Aldy [70] with 2 out of 21 errors. Astrolab [71] and Stargazer [72] both performed worse with 9 and 10 errors, respectively. Similar results were reported by Twesigomwe et al. in a CYP2D6 specific study [73]. All tools use the variants and *-allele translation in PharmVar as their database. However, PharmVar is updated continuously leading to potential differences in assignments if not every tool is updated at the same time. Additionally, depending on the data on which the tools are trained and tested, they might be more sensitive for specific variants and alleles [60]. Even the most extensive tools in regard to *-haplotype calling might not be suitable for clinical practice as they contain variants of which the effect is unknown. Furthermore, for clinical implementation only, the variants of known functional effect are relevant, resulting in a need for translation tools focused on clinical implementation only. However, selecting which variants are of direct clinical relevance remains challenging and requires attention and standardization [37,74]. SNV panels are usually designed to contain known variants, often with known clinical effect. This makes them easy to implement in clinical practice with standardized variant to haplotype translations. Sequencing data, on the other hand, contain more variant information and allow for the extraction of additional variants should they become of interest. Therefore, sequencing based translations can be updated with the development of more guidelines and insights into variant effects. Nonetheless, this does come with more intense data processing.

5.2. Imputation

To expand the number of interrogated SNVs, imputation can be used for technologies that do not cover the entire genome (Table 2). Imputation is predominantly used in GWAS analysis and to expand the PGx panel in genome wide arrays [8]. With imputation, the presence of a genomic variant is inferred based on the absence or presence of a linked SNV. These predictions all come with a probability for the occurrence of the SNV of interest. Often only imputations with a high probability are included (e.g., >90%) to avoid inaccurate assignments. Reisberg et al. have shown that imputation accuracies as high as 99% can be reached for PGx variants [8]. Nevertheless, a probability of 90% also means that there is a 10% change that the imputed variant is not correct. While this is certainly acceptable for population studies, it is not sufficient for tailoring drug treatment in an individual patient. Furthermore, to reach high imputation accuracy, an imputation dataset specific to each patient’s ethnical background is needed as the level of linkage disequilibrium (LD) between two SNVs can differ between different populations [8,75]. One clear example of the differences in LD between populations is the HLA tagging SNVs; to identify the HLA-A*3101 allele, associated with carbamazepine toxicity, a linked SNV is used. In Caucasians, the rs1061235 (NC_000006.12:g.29945521A>T) variant is in full LD with the *3101 haplotype, therefore the presence of the HLA-A*3101 allele can be inferred based on the presence of the rs1061235 variant [76]. However, in the Asian population, this variant is not in LD with the *3101 allele. For individuals of Asian descent, the rs1633021 (NC_000006.12:g.29779092T>C) variant can be used as a linked SNV as this variant is in LD with HLA-A*3101 in this population [77]. Using the Caucasian-linked SNV in the wrong population can lead to errors in the inferred haplotype, phenotype and ultimately lead to treatment errors. Therefore, the application of imputation should be limited to research purposes until the reliability for an individual patient has been proven.

5.3. Haplotype Phasing

In addition to variant expansion, imputation can also be used for haplotype phasing (Table 2). With haplotype phasing it can be determined if variants are located on the same allele or if they are on different alleles, potentially leading to differences in phenotype assignment [10,75,78]. The problem of phasing only exists when >1 heterozygous variant is present. However, given the polymorphic nature of many pharmacogenes, the likelihood of identifying multiple heterozygous variants within the gene locus of interest is highly likely [79]. The *-nomenclature is designed to describe only the variants in one allele, assuming that the variants have been phased into two separate alleles. Nonetheless, there are no clear guidelines available describing in what manner variants should be assigned to either one of the alleles. Both CPIC and the DPWG report which diplotypes translate into which phenotypes and occasionally which variants are needed to assign a specific haplotype but not on the phasing of variants. Haplotype phasing can, however, make the difference between a Poor Metabolizer (two loss-of-function variants on different alleles) and an Intermediate Metabolizer (two loss-of-function variants on the same allele and no variants on the other allele). Especially since many pharmacogenomic haplotypes are characterized by multiple variants, in theory all of them can have an impact on protein function. Including phasing in pharmacogenomics can, in some patients, improve haplotype assignments and therefore phenotype prediction. One example is the CYP2B6 gene for which phasing has been shown to be relevant [10]. When the rs3745274 (NC_000019.10:g.41006936G>T) and the rs2279343 (NC_000019.10:g.41009358A>G) variant are both detected, conventional methods assume they are located on the same allele based on linkage disequilibrium and assign a *6 haplotype (Figure 1). Additionally, CYP2B6*6 is the more common haplotype in most populations. The CYP2B6*6 haplotype occurs around 10% in Asians and up to 40% in the African population. CYP2B6*4 and *9 occur between 0 and 5% in all populations [80]. Clinical data show that in 1.5% of the individuals who carry both these variants, they are located on different alleles, resulting in a *4/*9 haplotype. In this case, both a CYP2B6*1/*6 and a CYP2B6*4/*9 call result in the same phenotype making the clinical impact limited [5]. However, for other variants, this might not be the case. Imputation could be used to infer haplotype phasing, by using the linkage between two observed SNVs to predict if they are located on the same allele or on different alleles. Nevertheless, the same limitations to imputation as described above apply [10].
With NGS, it is possible to phase read to their allele of origin without the need for pedigree information or computational phasing. For NGS, linked reads can be used for this purpose. Linked-read sequencing is based on partitioning the DNA with barcoded gel beads resulting in barcoded short fragments which can then be sequenced with conventional short-read methods. Due to the barcodes, every read can be linked back to the original position and artificial long input DNA can be reconstructed [81]. For long-read sequencing, the length of the reads in itself can be utilized for haplotype phasing. By overlapping the long-reads, large haploblocks can be formed. Nonetheless, large regions which are homozygous can cause an inability to phase large haploblocks [64].

5.4. Structural Variants

It has been shown that the majority of pharmacogenes is largely characterized by complex regions, such as CNVs, structural rearrangements and repetitive sections (Table 1) [82]. Full gene deletions or duplications of CYP2D6 occur in 5–10% of the population [83]. Nonetheless, not all arrays contain probes that can directly detect CNVs. To still be able to obtain CNV data from microarrays, several tools are available. Examples of these types of tools are PennCNV [84], QuantiSNP [85], GenoCN [86] and Nexus [87]. These tools make use of the B allele frequency and the log R ratio which are extracted from the array data. SNP array make use of common SNPs indicated by an A or B (allele) variant. In a normal situation of two copies, there is either an AA, AB or BB at a specific locus. The presence of a deletion or a duplication can be derived from an aberrant number of either an A or a B allele frequency, which is reflected in the B allele frequency parameter. The log R ratio is the log of the ratio between observed and expected intensity values at each variant. It reflects the intensity of the signal at each variant site, a deviation from the expected intensity signal. [88,89]. While results seem promising [8], full validation of these tools for PGx is still needed.
NGS data can resolve all SNVs in the sequenced region, nevertheless, it is difficult to assess CNVs based on sequencing data alone. To aid CNV calling with NGS, several tools have been developed; XHMM [90], CoNIFER [91], Varseq [92] and CNVnator [93], all of which use sequencing depth as an indication of a gene deletion or duplication. These tools do require large datasets and a sufficient range in depth to identify CNVs; details of the use of sequencing depth for CNV calling have been reviewed previously [90,94,95]. Yao et al. tested these three tools on their performance on CNVs of different sizes. Unfortunately, the agreement between the methods was low and there was a bias towards the smaller CNVs as opposed to large CNVs, potentially caused by the limitations of read length [94]. Nonetheless, the advancements in this field evolve rapidly, leading to several laboratories which are now able to reliably identify CNV based on short-read sequencing data. Efforts of using sequencing based CNV calling for CYP2D6 have shown mixed results [51,56]. Cohn et al. were able to accurately determine CNV status for 87 out of 98 patients. For nine, the CNV calls were inconclusive, and for a further two, there was a discrepant call between sequencing-based CNV calling and a targeted panel [51]. Gulilat et al., including 235 subjects, were able to confirm all sequencing-based CNV variants with panel-based testing [56].
The distinction between a pharmacogene and a pseudogene can be even more challenging. For example, CYP2D6 and CYP2D7 share >98% of their sequence, making it difficult to determine from which gene a sequencing read originates [67,96,97]. Due to the relatively short reads (100–200 bp), these complex regions cannot always be well characterized by NGS as reads are not long enough to distinguish between different locations in the complex region [98]. Long-read sequencing, on the other hand, allows for unambiguous mapping of a sequence to the gene of origin without interference of pseudogenes. Additionally, complex regions can be solved in one long read. Indels (deletions or duplications of 1–1000 base pairs) cannot always be accurately determined with short reads as the indel length might surpass the maximum read length. With a read length around 10 kbp in long-read sequencing any structural variants within this maximum length can be covered in one read [64,67,68].
The detection of CNVs in CYP2D6 is routine in clinical practice. However, the full characterization of the complexity of pharmacogenes is still in the research phase, relying on the further development of long-read sequencing and bioinformatic tools.

5.5. Variants of Unknown Effect

A clear benefit of sequencing over SNV panels is the increase in the number of variants that can be identified. While SNV-panel approaches remain limited to the pre-selected variants, sequencing data can help identify variants in the entire sequenced region, including rare variants (Table 2). As mentioned previously, over 90% of the identified variants in pharmacogenes are classified as a rare variant [99,100,101] (Table 1). Additionally, rare variants are expected to be more deleterious than common variants resulting in a potential higher impact on protein function [82]. To collect the most data on rare and novel variants, a WGS or targeted whole gene sequencing approach would be most suitable. Nevertheless, due to a lack of knowledge regarding the impact of these variants, they cannot yet be applied in clinical practice [83,98,101,102,103]. As they are by definition not commonly observed, it is difficult to assign a functional effect. Several strategies have been proposed to detect the impact of rare variants, of which the most common options are the use of cell-line models, in silico predictions or by studying patients displaying the most extreme phenotypes [104]. For clinical application, in vivo studies are most suitable. However, due to the low frequency of these variants, it is nearly impossible to have an appropriate sample size to study these variants [99,105]. In vitro analyses are more easily accessible and generally show a good indication of the effect of a particular variant. Nonetheless, in vitro findings can deviate from the in vivo situation and can still be too laborious for high throughput analyses. Therefore, for high throughput variant predictions, the in silico approach is most desirable. In silico models are based on sequence conservation, the physiochemical and crystal structure of the protein, or on evolutionary scores [102,106]. One, or a combination, of these factors is used to predict the impact the variant will have on enzyme function. To assess the applicability of in silico tools in pharmacogenetics, Han et al. conducted a study to test the accuracy of these tools. They showed that for 10 selected SNPs, the best models accurately predicted the functionality of 80% of the SNPs [107]. In addition, Hao et al. showed that 68% of non-synonymous SNPs in phase II enzyme genes were correctly predicted to be damaging [108]. These results indicate that the applicability of these assays is still limited. Ultimately, collecting more genetic, accompanying clinical data and better prediction models can help us understand the role of these rare variants to be able to use them in clinical practice.

5.6. Pharmacogenomics and Disease Genes

Moreover, variants that are disease predictors can be encountered. Complex examples of this are genes which are both pharmacogenes as well as disease-causing genes. RYR1 is linked to an increased risk of malignant hyperthermia which could classify it as a disease gene and, as such, it is included in the ACMG guidelines [109]. However, one of the factors that could cause the MH in susceptible patients is volatile anesthetics which could classify RYR1 as a pharmacogene, as well as it interacting with drugs [110].
In summary, the number of available genotyping technologies for PGx has evolved rapidly in recent years and continues to expand. Ultimately, selecting the right technology is not a matter of fact but a matter of choosing the right technique for the right problem.

Author Contributions

Writing—original draft preparation, M.v.d.L. and J.J.S.; writing—review and editing, H.-J.G. and M.K. All authors have read and agreed to the published version of the manuscript.

Funding

The research leading to these results has received funding from the European Community’s Horizon 2020 Program under grant agreement no. 668353 (U-PGx).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kirchheiner, J.; Brøsen, K.; Dahl, M.L.; Gram, L.F.; Kasper, S.; Roots, I.; Sjöqvist, F.; Spina, E.; Brockmöller, J. CYP2D6 and CYP2C19 genotype-based dose recommendations for antidepressants: A first step towards subpopulation-specific dosages. Acta Psychiatr. Scand. 2001, 104, 173–192. [Google Scholar] [CrossRef] [PubMed]
  2. Swen, J.J.; Wilting, I.; de Goede, A.L.; Grandia, L.; Mulder, H.; Touw, D.J.; de Boer, A.; Conemans, J.M.; Egberts, T.C.; Klungel, O.H.; et al. Pharmacogenetics: From bench to byte. Clin. Pharmacol. Ther. 2008, 83, 781–787. [Google Scholar] [CrossRef] [PubMed]
  3. Hoffman, J.M.; Haidar, C.E.; Wilkinson, M.R.; Crews, K.R.; Baker, D.K.; Kornegay, N.M.; Yang, W.; Pui, C.H.; Reiss, U.M.; Gaur, A.H.; et al. PG4KDS: A model for the clinical implementation of pre-emptive pharmacogenetics. Am. J. Med. Genet. C Semin. Med. Genet. 2014, 166c, 45–55. [Google Scholar] [CrossRef] [Green Version]
  4. Dutch Pharmacogenetics Working group. Pharmacogenetics Guidelines; Royal Dutch Pharmacists Association (KNMP Kennisbank): The Hague, The Netherlands, 2020. [Google Scholar]
  5. Clinical Pharmacogenetics Implementation Consortium. CPIC-guidelines. Available online: https://cpicpgx.org/ (accessed on 16 October 2020).
  6. National Human Genome Research Institute. DNA Sequencing Costs: Data. Available online: https://www.genome.gov/about-genomics/fact-sheets/DNA-Sequencing-Costs-Data (accessed on 16 October 2020).
  7. Hallberg, P.; Yue, Q.Y.; Eliasson, E.; Melhus, H.; Ås, J.; Wadelius, M. SWEDEGENE-a Swedish nation-wide DNA sample collection for pharmacogenomic studies of serious adverse drug reactions. Pharm. J. 2020, 20, 579–585. [Google Scholar] [CrossRef] [PubMed]
  8. Reisberg, S.; Krebs, K.; Kals, M.; Magi, R.; Metsalu, K.; Lauschke, V.M.; Vilo, J.; Milani, L. Translating genotype data of 44,000 biobank participants into clinical pharmacogenetic recommendations: Challenges and solutions. Genet. Med. 2019, 21, 1345–1354. [Google Scholar] [CrossRef] [Green Version]
  9. Cousin, M.A.; Matey, E.T.; Blackburn, P.R.; Boczek, N.J.; McAllister, T.M.; Kruisselbrink, T.M.; Babovic-Vuksanovic, D.; Lazaridis, K.N.; Klee, E.W. Pharmacogenomic findings from clinical whole exome sequencing of diagnostic odyssey patients. Mol. Genet. Genom. Med. 2017, 5, 269–279. [Google Scholar] [CrossRef] [Green Version]
  10. Van der Lee, M.; Allard, W.G.; Bollen, S.; Santen, G.W.E.; Ruivenkamp, C.A.L.; Hoffer, M.J.V.; Kriek, M.; Guchelaar, H.J.; Anvar, S.Y.; Swen, J.J. Repurposing of Diagnostic Whole Exome Sequencing Data of 1,583 Individuals for Clinical Pharmacogenetics. Clin. Pharmacol. Ther. 2020, 107, 617–627. [Google Scholar] [CrossRef] [Green Version]
  11. Krebs, K.; Milani, L. Translating pharmacogenomics into clinical decisions: Do not let the perfect be the enemy of the good. Hum. Genom. 2019, 13, 39. [Google Scholar] [CrossRef] [Green Version]
  12. Ameur, A.; Kloosterman, W.P.; Hestand, M.S. Single-molecule sequencing: Towards clinical applications. Trends Biotechnol. 2019, 37, 72–85. [Google Scholar] [CrossRef]
  13. Sakamoto, Y.; Sereewattanawoot, S.; Suzuki, A. A new era of long-read sequencing for cancer genomics. J Hum. Genet. 2020, 65, 3–10. [Google Scholar] [CrossRef] [Green Version]
  14. Aganezov, S.; Goodwin, S.; Sherman, R.M.; Sedlazeck, F.J.; Arun, G.; Bhatia, S.; Lee, I.; Kirsche, M.; Wappel, R.; Kramer, M.; et al. Comprehensive analysis of structural variants in breast cancer genomes using single-molecule sequencing. Genome Res. 2020, 30, 1258–1273. [Google Scholar] [CrossRef] [PubMed]
  15. Sakamoto, Y.; Xu, L.; Seki, M.; Yokoyama, T.T.; Kasahara, M.; Kashima, Y.; Ohashi, A.; Shimada, Y.; Motoi, N.; Tsuchihara, K.; et al. Long-read sequencing for non-small-cell lung cancer genomes. Genome Res. 2020, 30, 1243–1257. [Google Scholar] [CrossRef] [PubMed]
  16. Just, K.S.; Steffens, M.; Swen, J.J.; Patrinos, G.P.; Guchelaar, H.J.; Stingl, J.C. Medical education in pharmacogenomics-results from a survey on pharmacogenetic knowledge in healthcare professionals within the European pharmacogenomics clinical implementation project Ubiquitous Pharmacogenomics (U-PGx). Eur. J. Clin. Pharmacol. 2017, 73, 1247–1252. [Google Scholar] [CrossRef] [PubMed]
  17. Rollinson, V.; Turner, R.; Pirmohamed, M. Pharmacogenomics for Primary Care: An Overview. Genes 2020, 11, 1337. [Google Scholar] [CrossRef]
  18. Bank, P.C.D.; Swen, J.J.; Guchelaar, H.J. Implementation of Pharmacogenomics in Everyday Clinical Settings. Stud. Surf. Sci. Catal. 2018, 83, 219–246. [Google Scholar] [CrossRef]
  19. PharmGKB. DPWG: Dutch Pharmacogenetics Working Group. Available online: https://www.pharmgkb.org/page/dpwg (accessed on 23 June 2020).
  20. Illumina Inc. VeraCode ADME Core Panel. Available online: https://www.illumina.com/documents/products/datasheets/datasheet_veracode_adme_core_panel.pdf (accessed on 21 July 2020).
  21. Arbitrio, M.; Di Martino, M.T.; Scionti, F.; Agapito, G.; Guzzi, P.H.; Cannataro, M.; Tassone, P.; Tagliaferri, P. DMET™ (Drug Metabolism Enzymes and Transporters): A pharmacogenomic platform for precision medicine. Oncotarget 2016, 7, 54028–54050. [Google Scholar] [CrossRef]
  22. ThermoFisher Scientific. Pharmacoscan Assay. Available online: https://www.thermofisher.com/order/catalog/product/903010TS (accessed on 16 October 2020).
  23. Gabriel, S.; Ziaugra, L.; Tabbaa, D. SNP genotyping using the Sequenom MassARRAY iPLEX platform. Curr. Protoc. Hum. Genet. 2009. [Google Scholar] [CrossRef]
  24. Spierings, G.; Dunbar, S.A. Pharmacogenetics using Luminex(R) xMAP(R) technology: A method for developing a custom multiplex single nucleotide polymorphism mutation assay. Methods Mol. Biol. 2013, 1015, 115–126. [Google Scholar] [CrossRef]
  25. Chen, C.; Li, S.; Lu, X.; Tan, B.; Huang, C.; Qin, L. High resolution melting method to detect single nucleotide polymorphism of VKORC1 and CYP2C9. Int. J. Clin. Exp. Pathol. 2014, 7, 2558–2564. [Google Scholar]
  26. Jannetto, P.J.; Laleli-Sahin, E.; Wong, S.H. Pharmacogenomic genotyping methodologies. Clin. Chem. Lab. Med. 2004, 42, 1256–1264. [Google Scholar] [CrossRef]
  27. Ghasemi, Z.; Hashemi, M.; Ejabati, M.; Ebrahimi, S.M.; Kheiri Manjili, H.; Sharafi, A.; Ramazani, A. Development of a High-Resolution Melting Analysis Method for CYP2C19*17 Genotyping in Healthy Volunteers. Avicenna J. Med. Biotechnol. 2016, 8, 193–199. [Google Scholar] [PubMed]
  28. Mukerjee, G.; Huston, A.; Kabakchiev, B.; Piquette-Miller, M.; van Schaik, R.; Dorfman, R. User considerations in assessing pharmacogenomic tests and their clinical support tools. NPJ Genom. Med. 2018, 3, 26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Vilella, A. Next-Generation-Sequencing, v1.5.7. Available online: https://docs.google.com/spreadsheets/d/1GMMfhyLK0-q8XkIo3YxlWaZA5vVMuhU1kg41g4xLkXc/htmlview?hl=en_GB (accessed on 27 October 2020).
  30. Agenda Bioscience. VeriDose Core Panel. Available online: https://agenabio.com/products/panel/veridose-core-panel/ (accessed on 2 November 2020).
  31. Huang, S.M.; Goodsaid, F.; Rahman, A.; Frueh, F.; Lesko, L.J. Application of pharmacogenomics in clinical pharmacology. Toxicol. Mech. Methods. 2006, 16, 89–99. [Google Scholar] [CrossRef] [PubMed]
  32. ThermoFisher Scientifuic. Axiom Pharmacofocus. Available online: https://www.thermofisher.com/order/catalog/product/952425?SID=srch-hj-952425#/952425?SID=srch-hj-952425 (accessed on 3 November 2020).
  33. Pulley, J.M.; Denny, J.C.; Peterson, J.F.; Bernard, G.R.; Vnencak-Jones, C.L.; Ramirez, A.H.; Delaney, J.T.; Bowton, E.; Brothers, K.; Johnson, K.; et al. Operational implementation of prospective genotyping for personalized medicine: The design of the Vanderbilt PREDICT project. Clin. Pharmacol. Ther. 2012, 92, 87–95. [Google Scholar] [CrossRef] [Green Version]
  34. Van der Wouden, C.H.; Cambon-Thomsen, A.; Cecchin, E.; Cheung, K.C.; Davila-Fajardo, C.L.; Deneer, V.H.; Dolzan, V.; Ingelman-Sundberg, M.; Jonsson, S.; Karlsson, M.O.; et al. Implementing Pharmacogenomics in Europe: Design and Implementation Strategy of the Ubiquitous Pharmacogenomics Consortium. Clin. Pharmacol. Ther. 2017, 101, 341–358. [Google Scholar] [CrossRef]
  35. ThermoFisher Scientific. Open Array Technology Overview. Available online: https://www.thermofisher.com/nl/en/home/life-science/pcr/real-time-pcr/real-time-openarray/open-array-technology.html (accessed on 3 November 2020).
  36. Eadon, M.T.; Desta, Z.; Levy, K.D.; Decker, B.S.; Pierson, R.C.; Pratt, V.M.; Callaghan, J.T.; Rosenman, M.B.; Carpenter, J.S.; Holmes, A.M.; et al. Implementation of a pharmacogenomics consult service to support the INGENIOUS trial. Clin. Pharmacol. Ther. 2016, 100, 63–66. [Google Scholar] [CrossRef]
  37. Van der Wouden, C.H.; van Rhenen, M.H.; Jama, W.O.M.; Ingelman-Sundberg, M.; Lauschke, V.M.; Konta, L.; Schwab, M.; Swen, J.J.; Guchelaar, H.J. Development of the PGx-Passport: A Panel of Actionable Germline Genetic Variants for Pre-emptive Pharmacogenetic Testing. Clin. Pharmacol. Ther. 2019, 106, 866–873. [Google Scholar] [CrossRef]
  38. Biosearch Technologies. SNPline Genotyping Automation. Available online: https://www.biosearchtech.com/products/instruments-and-consumables/genotyping-instruments/snpline-genotyping-automation (accessed on 3 November 2020).
  39. Illumina Inc. Illumina Global Screening Array. Available online: https://emea.illumina.com/products/by-type/microarray-kits/infinium-global-screening.html (accessed on 23 June 2020).
  40. ThermoFisher Scientific. Axiom Genotyping Solutions. Available online: https://assets.thermofisher.com/TFS-Assets/LSG/brochures/axiom_solution_brochure.pdf (accessed on 3 November 2020).
  41. Ingelman-Sundberg, M. Pharmacogenetics of cytochrome P450 and its applications in drug therapy: The past, present and future. Trends Pharmacol. Sci. 2004, 25, 193–200. [Google Scholar] [CrossRef]
  42. Gaedigk, A.; Ingelman-Sundberg, M.; Miller, N.A.; Leeder, J.S.; Whirl-Carrillo, M.; Klein, T.E. The Pharmacogene Variation (PharmVar) Consortium: Incorporation of the Human Cytochrome P450 (CYP) Allele Nomenclature Database. Clin. Pharmacol. Ther. 2018, 103, 399–401. [Google Scholar] [CrossRef] [Green Version]
  43. Broadinstitute. GnomAD. Available online: https://gnomad.broadinstitute.org/ (accessed on 26 October 2020).
  44. Mardis, E.R. Next-generation sequencing platforms. Annu. Rev. Anal. Chem. 2013, 6, 287–303. [Google Scholar] [CrossRef] [Green Version]
  45. Slatko, B.E.; Gardner, A.F.; Ausubel, F.M. Overview of Next-Generation Sequencing Technologies. Curr. Protoc. Mol. Biol. 2018, 122, e59. [Google Scholar] [CrossRef] [PubMed]
  46. Levy, S.E.; Myers, R.M. Advancements in Next-Generation Sequencing. Annu. Rev. Genom. Hum. Genet. 2016, 17, 95–115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Yang, W.; Wu, G.; Broeckel, U.; Smith, C.A.; Turner, V.; Haidar, C.E.; Wang, S.; Carter, R.; Karol, S.E.; Neale, G.; et al. Comparison of genome sequencing and clinical genotyping for pharmacogenes. Clin. Pharmacol. Ther. 2016, 100, 380–388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Londin, E.R.; Clark, P.; Sponziello, M.; Kricka, L.J.; Fortina, P.; Park, J.Y. Performance of exome sequencing for pharmacogenomics. Pers. Med. 2014, 12, 109–115. [Google Scholar] [CrossRef] [Green Version]
  49. Rasmussen-Torvik, L.J.; Almoguera, B.; Doheny, K.F.; Freimuth, R.R.; Gordon, A.S.; Hakonarson, H.; Hawkins, J.B.; Husami, A.; Ivacic, L.C.; Kullo, I.J.; et al. Concordance between Research Sequencing and Clinical Pharmacogenetic Genotyping in the eMERGE-PGx Study. J. Mol. Diagn. 2017, 19, 561–566. [Google Scholar] [CrossRef] [Green Version]
  50. Ng, D.; Hong, C.S.; Singh, L.N.; Johnston, J.J.; Mullikin, J.C.; Biesecker, L.G. Assessing the capability of massively parallel sequencing for opportunistic pharmacogenetic screening. Genet. Med. 2017, 19, 357–361. [Google Scholar] [CrossRef] [Green Version]
  51. Cohn, I.; Paton, T.A.; Marshall, C.R.; Basran, R.; Stavropoulos, D.J.; Ray, P.N.; Monfared, N.; Hayeems, R.Z.; Meyn, M.S.; Bowdin, S.; et al. Genome sequencing as a platform for pharmacogenetic genotyping: A pediatric cohort study. NPJ Genom. Med. 2017, 2, 19. [Google Scholar] [CrossRef]
  52. Gordon, A.S.; Fulton, R.S.; Qin, X.; Mardis, E.R.; Nickerson, D.A.; Scherer, S. PGRNseq: A targeted capture sequencing panel for pharmacogenetic research and implementation. Pharm. Genom. 2016, 26, 161–168. [Google Scholar] [CrossRef]
  53. Chua, E.W.; Cree, S.L.; Ton, K.N.; Lehnert, K.; Shepherd, P.; Helsby, N.; Kennedy, M.A. Cross-Comparison of Exome Analysis, Next-Generation Sequencing of Amplicons, and the iPLEX((R)) ADME PGx Panel for Pharmacogenomic Profiling. Front. Pharmacol. 2016, 7, 1. [Google Scholar] [CrossRef] [Green Version]
  54. Ingelman-Sundberg, M.; Sim, S.C. Intronic polymorphisms of cytochromes P450. Hum. Genom. 2010, 4, 402–405. [Google Scholar] [CrossRef] [Green Version]
  55. Bush, W.S.; Crosslin, D.R.; Owusu-Obeng, A.; Wallace, J.; Almoguera, B.; Basford, M.A.; Bielinski, S.J.; Carrell, D.S.; Connolly, J.J.; Crawford, D.; et al. Genetic variation among 82 pharmacogenes: The PGRNseq data from the eMERGE network. Pharmacol. Ther. 2016, 100, 160–169. [Google Scholar] [CrossRef] [PubMed]
  56. Gulilat, M.; Lamb, T.; Teft, W.A.; Wang, J.; Dron, J.S.; Robinson, J.F.; Tirona, R.G.; Hegele, R.A.; Kim, R.B.; Schwarz, U.I. Targeted next generation sequencing as a tool for precision medicine. BMC Med. Genom. 2019, 12, 81. [Google Scholar] [CrossRef] [PubMed]
  57. Aquilante, C.L.; Kao, D.P.; Trinkley, K.E.; Lin, C.T.; Crooks, K.R.; Hearst, E.C.; Hess, S.J.; Kudron, E.L.; Lee, Y.M.; Liko, I.; et al. Clinical implementation of pharmacogenomics via a health system-wide research biobank: The University of Colorado experience. Pharmacogenomics 2020, 21, 375–386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. National Institute of Health. AllofUs Research Program. Available online: https://allofus.nih.gov/ (accessed on 23 October 2020).
  59. Precision Medicine Initiative Work Group. The Precision Medicine Initiative Cohort Program—Building a Research Foundation for 21st Century Medicine; National Institutes of Health. 2015. Available online: https://www.nih.gov/sites/default/files/research-training/initiatives/pmi/pmi-working-group-report-20150917-2.pdf (accessed on 23 October 2020).
  60. Caspar, S.M.; Schneider, T.; Meienberg, J.; Matyas, G. Added Value of Clinical Sequencing: WGS-Based Profiling of Pharmacogenes. Int. J. Mol. Sci. 2020, 21, 2308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Mantere, T.; Kersten, S.; Hoischen, A. Long-Read Sequencing Emerging in Medical Genetics. Front. Genet. 2019, 10, 426. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Rhoads, A.; Au, K.F. PacBio Sequencing and Its Applications. Genom. Proteom. Bioinform. 2015, 13, 278–289. [Google Scholar] [CrossRef] [Green Version]
  63. Bowden, R.; Davies, R.W.; Heger, A.; Pagnamenta, A.T.; de Cesare, M.; Oikkonen, L.E.; Parkes, D.; Freeman, C.; Dhalla, F.; Patel, S.Y.; et al. Sequencing of human genomes with nanopore technology. Nat. Commun. 2019, 10, 1869. [Google Scholar] [CrossRef]
  64. Wenger, A.M.; Peluso, P.; Rowell, W.J.; Chang, P.C.; Hall, R.J.; Concepcion, G.T.; Ebler, J.; Fungtammasan, A.; Kolesnikov, A.; Olson, N.D.; et al. Accurate circular consensus long-read sequencing improves variant detection and assembly of a human genome. Nat. Biotechnol. 2019, 37, 1155–1162. [Google Scholar] [CrossRef]
  65. Schüle, B.; McFarland, K.N.; Lee, K.; Tsai, Y.C.; Nguyen, K.D.; Sun, C.; Liu, M.; Byrne, C.; Gopi, R.; Huang, N.; et al. Parkinson’s disease associated with pure ATXN10 repeat expansion. NPJ Park. Dis. 2017, 3, 27. [Google Scholar] [CrossRef]
  66. Ardui, S.; Race, V.; Zablotskaya, A.; Hestand, M.S.; Van Esch, H.; Devriendt, K.; Matthijs, G.; Vermeesch, J.R. Detecting AGG Interruptions in Male and Female FMR1 Premutation Carriers by Single-Molecule Sequencing. Hum. Mutat. 2017, 38, 324–331. [Google Scholar] [CrossRef] [Green Version]
  67. Qiao, W.; Yang, Y.; Sebra, R.; Mendiratta, G.; Gaedigk, A.; Desnick, R.J.; Scott, S.A. Long-Read Single Molecule Real-Time Full Gene Sequencing of Cytochrome P450-2D6. Hum. Mutat. 2016, 37, 315–323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Buermans, H.P.; Vossen, R.H.; Anvar, S.Y.; Allard, W.G.; Guchelaar, H.J.; White, S.J.; den Dunnen, J.T.; Swen, J.J.; van der Straaten, T. Flexible and Scalable Full-Length CYP2D6 Long Amplicon PacBio Sequencing. Hum. Mutat. 2017, 38, 310–316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Robarge, J.D.; Li, L.; Desta, Z.; Nguyen, A.; Flockhart, D.A. The star-allele nomenclature: Retooling for translational genomics. Clin. Pharmacol. Ther. 2007, 82, 244–248. [Google Scholar] [CrossRef] [PubMed]
  70. Numanagic, I.; Malikic, S.; Ford, M.; Qin, X.; Toji, L.; Radovich, M.; Skaar, T.C.; Pratt, V.M.; Berger, B.; Scherer, S.; et al. Allelic decomposition and exact genotyping of highly polymorphic and structurally variant genes. Nat. Commun. 2018, 9, 828. [Google Scholar] [CrossRef] [PubMed]
  71. Twist, G.P.; Gaedigk, A.; Miller, N.A.; Farrow, E.G.; Willig, L.K.; Dinwiddie, D.L.; Petrikin, J.E.; Soden, S.E.; Herd, S.; Gibson, M.; et al. Constellation: A tool for rapid, automated phenotype assignment of a highly polymorphic pharmacogene, CYP2D6, from whole-genome sequences. NPJ Genom. Med. 2016, 1, 15007. [Google Scholar] [CrossRef] [PubMed]
  72. Lee, S.B.; Wheeler, M.M.; Thummel, K.E.; Nickerson, D.A. Calling Star Alleles with Stargazer in 28 Pharmacogenes With Whole Genome Sequences. Clin. Pharmacol. Ther. 2019, 106, 1328–1337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Twesigomwe, D.; Wright, G.E.B.; Drögemöller, B.I.; da Rocha, J.; Lombard, Z.; Hazelhurst, S. A systematic comparison of pharmacogene star allele calling bioinformatics algorithms: A focus on CYP2D6 genotyping. NPJ Genom. Med. 2020, 5, 30. [Google Scholar] [CrossRef]
  74. Pratt, V.M.; Everts, R.E.; Aggarwal, P.; Beyer, B.N.; Broeckel, U.; Epstein-Baak, R.; Hujsak, P.; Kornreich, R.; Liao, J.; Lorier, R.; et al. Characterization of 137 Genomic DNA Reference Materials for 28 Pharmacogenetic Genes: A GeT-RM Collaborative Project. J. Mol. Diagn. 2016, 18, 109–123. [Google Scholar] [CrossRef] [Green Version]
  75. Browning, S.R.; Browning, B.L. Haplotype phasing: Existing methods and new developments. Nat. Rev. Genet. 2011, 12, 703–714. [Google Scholar] [CrossRef] [Green Version]
  76. McCormack, M.; Alfirevic, A.; Bourgeois, S.; Farrell, J.J.; Kasperavičiūtė, D.; Carrington, M.; Sills, G.J.; Marson, T.; Jia, X.; de Bakker, P.I.; et al. HLA-A*3101 and carbamazepine-induced hypersensitivity reactions in Europeans. N. Engl. J. Med. 2011, 364, 1134–1143. [Google Scholar] [CrossRef] [Green Version]
  77. Ozeki, T.; Mushiroda, T.; Yowang, A.; Takahashi, A.; Kubo, M.; Shirakata, Y.; Ikezawa, Z.; Iijima, M.; Shiohara, T.; Hashimoto, K.; et al. Genome-wide association study identifies HLA-A*3101 allele as a genetic risk factor for carbamazepine-induced cutaneous adverse drug reactions in Japanese population. Hum. Mol. Genet. 2011, 20, 1034–1041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Browning, S.R.; Browning, B.L. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am. J. Hum. Genet. 2007, 81, 1084–1097. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Jin, Y.; Wang, J.; Bachtiar, M.; Chong, S.S.; Lee, C.G.L. Architecture of polymorphisms in the human genome reveals functionally important and positively selected variants in immune response and drug transporter genes. Hum. Genom. 2018, 12, 43. [Google Scholar] [CrossRef] [PubMed]
  80. Zanger, U.M.; Klein, K.; Saussele, T.; Blievernicht, J.; Hofmann, M.H.; Schwab, M. Polymorphic CYP2B6: Molecular mechanisms and emerging clinical significance. Pharmacogenomics 2007, 8, 743–759. [Google Scholar] [CrossRef] [PubMed]
  81. Zhou, X.; Batzoglou, S.; Sidow, A.; Zhang, L. HAPDeNovo: A haplotype-based approach for filtering and phasing de novo mutations in linked read sequencing data. BMC Genom. 2018, 19, 467. [Google Scholar] [CrossRef] [PubMed]
  82. Ingelman-Sundberg, M.; Sim, S.C. Pharmacogenetic biomarkers as tools for improved drug therapy; emphasis on the cytochrome P450 system. Biochem. Biophys. Res. Commun. 2010, 396, 90–94. [Google Scholar] [CrossRef] [PubMed]
  83. Zhou, Y.; Ingelman-Sundberg, M.; Lauschke, V.M. Worldwide distribution of cytochrome P450 alleles: A meta-analysis of population-scale sequencing projects. Clin. Pharmacol. Ther. 2017, 102, 688–700. [Google Scholar] [CrossRef] [Green Version]
  84. Wang, K.; Li, M.; Hadley, D.; Liu, R.; Glessner, J.; Grant, S.F.; Hakonarson, H.; Bucan, M. PennCNV: An integrated hidden Markov model designed for high-resolution copy number variation detection in whole-genome SNP genotyping data. Genome Res. 2007, 17, 1665–1674. [Google Scholar] [CrossRef] [Green Version]
  85. Colella, S.; Yau, C.; Taylor, J.M.; Mirza, G.; Butler, H.; Clouston, P.; Bassett, A.S.; Seller, A.; Holmes, C.C.; Ragoussis, J. QuantiSNP: An Objective Bayes Hidden-Markov Model to detect and accurately map copy number variation using SNP genotyping data. Nucleic Acids Res. 2007, 35, 2013–2025. [Google Scholar] [CrossRef] [Green Version]
  86. Sun, W.; Wright, F.A.; Tang, Z.; Nordgard, S.H.; Van Loo, P.; Yu, T.; Kristensen, V.N.; Perou, C.M. Integrated study of copy number states and genotype calls using high-density SNP arrays. Nucleic Acids Res. 2009, 37, 5365–5377. [Google Scholar] [CrossRef]
  87. Darvishi, K. Application of Nexus copy number software for CNV detection and analysis. Curr. Protoc. Hum. Genet. 2010. [Google Scholar] [CrossRef] [PubMed]
  88. Seiser, E.L.; Innocenti, F. Hidden Markov Model-Based CNV Detection Algorithms for Illumina Genotyping Microarrays. Cancer Inform. 2014, 13, 77–83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Dellinger, A.E.; Saw, S.M.; Goh, L.K.; Seielstad, M.; Young, T.L.; Li, Y.J. Comparative analyses of seven algorithms for copy number variant identification from single nucleotide polymorphism arrays. Nucleic Acids Res. 2010, 38, e105. [Google Scholar] [CrossRef]
  90. Fromer, M.; Moran, J.L.; Chambert, K.; Banks, E.; Bergen, S.E.; Ruderfer, D.M.; Handsaker, R.E.; McCarroll, S.A.; O’Donovan, M.C.; Owen, M.J.; et al. Discovery and statistical genotyping of copy-number variation from whole-exome sequencing depth. Am. J. Hum. Genet. 2012, 91, 597–607. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Krumm, N.; Sudmant, P.H.; Ko, A.; O’Roak, B.J.; Malig, M.; Coe, B.P.; Quinlan, A.R.; Nickerson, D.A.; Eichler, E.E. Copy number variation detection and genotyping from exome sequence data. Genome Res. 2012, 22, 1525–1532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. GoldenHelix VarSeq. Available online: https://www.goldenhelix.com/products/VarSeq/ (accessed on 25 November 2020).
  93. Abyzov, A.; Urban, A.E.; Snyder, M.; Gerstein, M. CNVnator: An approach to discover, genotype, and characterize typical and atypical CNVs from family and population genome sequencing. Genome Res. 2011, 21, 974–984. [Google Scholar] [CrossRef] [Green Version]
  94. Yao, R.; Zhang, C.; Yu, T.; Li, N.; Hu, X.; Wang, X.; Wang, J.; Shen, Y. Evaluation of three read-depth based CNV detection tools using whole-exome sequencing data. Mol. Cytogenet. 2017, 10, 30. [Google Scholar] [CrossRef] [Green Version]
  95. Tremmel, R.; Klein, K.; Battke, F.; Fehr, S.; Winter, S.; Scheurenbrand, T.; Schaeffeler, E.; Biskup, S.; Schwab, M.; Zanger, U.M. Copy number variation profiling in pharmacogenes using panel-based exome resequencing and correlation to human liver expression. Hum. Genet. 2019, 139, 137–149. [Google Scholar] [CrossRef]
  96. Gaedigk, A. Complexities of CYP2D6 gene analysis and interpretation. Int. Rev. Psychiatry 2013, 25, 534–553. [Google Scholar] [CrossRef]
  97. Gaedigk, A.; Jaime, L.K.; Bertino, J.S., Jr.; Berard, A.; Pratt, V.M.; Bradfordand, L.D.; Leeder, J.S. Identification of Novel CYP2D7-2D6 Hybrids: Non-Functional and Functional Variants. Front. Pharmacol. 2010, 1, 121. [Google Scholar] [CrossRef] [Green Version]
  98. Lauschke, V.M.; Milani, L.; Ingelman-Sundberg, M. Pharmacogenomic Biomarkers for Improved Drug Therapy-Recent Progress and Future Developments. AAPS J. 2017, 20, 4. [Google Scholar] [CrossRef] [PubMed]
  99. Kozyra, M.; Ingelman-Sundberg, M.; Lauschke, V.M. Rare genetic variants in cellular transporters, metabolic enzymes, and nuclear receptors can be important determinants of interindividual differences in drug response. Genet. Med. 2017, 19, 20–29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  100. Gordon, A.S.; Tabor, H.K.; Johnson, A.D.; Snively, B.M.; Assimes, T.L.; Auer, P.L.; Ioannidis, J.P.; Peters, U.; Robinson, J.G.; Sucheston, L.E.; et al. Quantifying rare, deleterious variation in 12 human cytochrome P450 drug-metabolism genes in a large-scale exome dataset. Hum. Mol. Genet. 2014, 23, 1957–1963. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Fujikura, K.; Ingelman-Sundberg, M.; Lauschke, V.M. Genetic variation in the human cytochrome P450 supergene family. Pharm. Genom. 2015, 25, 584–594. [Google Scholar] [CrossRef]
  102. Drogemoller, B.I.; Wright, G.E.; Warnich, L. Considerations for rare variants in drug metabolism genes and the clinical implications. Expert Opin. Drug Metab. Toxicol. 2014, 10, 873–884. [Google Scholar] [CrossRef]
  103. Ingelman-Sundberg, M.; Mkrtchian, S.; Zhou, Y.; Lauschke, V.M. Integrating rare genetic variants into pharmacogenetic drug response predictions. Hum. Genom. 2018, 12, 26. [Google Scholar] [CrossRef]
  104. Lauschke, V.M.; Ingelman-Sundberg, M. Requirements for comprehensive pharmacogenetic genotyping platforms. Pharmacogenomics 2016, 17, 917–924. [Google Scholar] [CrossRef]
  105. Lauschke, V.M.; Ingelman-Sundberg, M. How to Consider Rare Genetic Variants in Personalized Drug Therapy. Clin. Pharmacol. Ther. 2018, 103, 745–748. [Google Scholar] [CrossRef]
  106. Li, B.; Seligman, C.; Thusberg, J.; Miller, J.L.; Auer, J.; Whirl-Carrillo, M.; Capriotti, E.; Klein, T.E.; Mooney, S.D. In silico comparative characterization of pharmacogenomic missense variants. BMC Genom. 2014, 15 (Suppl. S4), S4. [Google Scholar] [CrossRef] [Green Version]
  107. Han, S.M.; Park, J.; Lee, J.H.; Lee, S.S.; Kim, H.; Han, H.; Kim, Y.; Yi, S.; Cho, J.Y.; Jang, I.J.; et al. Targeted Next-Generation Sequencing for Comprehensive Genetic Profiling of Pharmacogenes. Clin. Pharmacol. Ther. 2017, 101, 396–405. [Google Scholar] [CrossRef]
  108. Hao, D.; Xiao, P.; Chen, S. Phenotype prediction of nonsynonymous single nucleotide polymorphisms in human phase II drug/xenobiotic metabolizing enzymes: Perspectives on molecular evolution. Sci. China Life Sci. 2010, 53, 1252–1262. [Google Scholar] [CrossRef] [PubMed]
  109. Kalia, S.S.; Adelman, K.; Bale, S.J.; Chung, W.K.; Eng, C.; Evans, J.P.; Herman, G.E.; Hufnagel, S.B.; Klein, T.E.; Korf, B.R.; et al. Recommendations for reporting of secondary findings in clinical exome and genome sequencing, 2016 update (ACMG SF v2.0): A policy statement of the American College of Medical Genetics and Genomics. Genet. Med. 2017, 19, 249–255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  110. Gonsalves, S.G.; Dirksen, R.T.; Sangkuhl, K.; Pulk, R.; Alvarellos, M.; Vo, T.; Hikino, K.; Roden, D.; Klein, T.E.; Poler, S.M.; et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for the Use of Potent Volatile Anesthetic Agents and Succinylcholine in the Context of RYR1 or CACNA1S Genotypes. Clin. Pharmacol. Ther. 2019, 105, 1338–1344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Genes 11 01456 g001 550
Figure 1. Haplotype phasing in CYP2B6. Inability to phase the rs3745274 (NC_000019.10:g.41006936G>T) and the rs2279343 (NC_000019.10:g.41009358A>G) to the correct allele can result in differences in *-haplotype assignment. Panel (A) shows the most common situation in individuals who are heterozygous for both variants, the CYP2B6*1/*6 diplotype. Panel (B) shows the alternative conformation where variants are located on opposing alleles leading to a CYP2B6*4/*9 diplotype. Both situations result in the same predicted CYP2B6 metabolizer phenotype (intermediate metabolizer). Conventional methods using linkage disequilibrium assume the variants to be located on the same allele, resulting in a CYP2B6*6 assignment. wt: Wild type.
Figure 1. Haplotype phasing in CYP2B6. Inability to phase the rs3745274 (NC_000019.10:g.41006936G>T) and the rs2279343 (NC_000019.10:g.41009358A>G) to the correct allele can result in differences in *-haplotype assignment. Panel (A) shows the most common situation in individuals who are heterozygous for both variants, the CYP2B6*1/*6 diplotype. Panel (B) shows the alternative conformation where variants are located on opposing alleles leading to a CYP2B6*4/*9 diplotype. Both situations result in the same predicted CYP2B6 metabolizer phenotype (intermediate metabolizer). Conventional methods using linkage disequilibrium assume the variants to be located on the same allele, resulting in a CYP2B6*6 assignment. wt: Wild type.
Genes 11 01456 g001
Table
Table 1. Characteristics of pharmacogenes in CPIC and DPWG guidelines. Related drugs are all drugs with dose recommendations. A drug can be related to multiple genes and therefore counted more than once. The locus size and the number of known rare variants (the number of variants with a minor allele frequency < 1% in an aggregated population, including singletons) are extracted from Gnomad (GRCh38). Part of the locus defined as complex is the percentage of the locus defined as a repeat or segmental duplication extracted from UCSC browser (https://genome.ucsc.edu). CPIC: Clinical Pharmacogenetic Implementation Consortium, DPWG: Dutch Pharmacogenetics Working Group.
Table 1. Characteristics of pharmacogenes in CPIC and DPWG guidelines. Related drugs are all drugs with dose recommendations. A drug can be related to multiple genes and therefore counted more than once. The locus size and the number of known rare variants (the number of variants with a minor allele frequency < 1% in an aggregated population, including singletons) are extracted from Gnomad (GRCh38). Part of the locus defined as complex is the percentage of the locus defined as a repeat or segmental duplication extracted from UCSC browser (https://genome.ucsc.edu). CPIC: Clinical Pharmacogenetic Implementation Consortium, DPWG: Dutch Pharmacogenetics Working Group.
ProteinGeneRelated DrugsLocus Size (bp)Rare Variants, n (% of Known Variants)Part of Locus Defined as Complex, %(bp)
CPICDPWG
CACNA1SCACNA1S7-73,0552520 (98%)33.3
CFTRCFTR1-250,1871684 (99%)42.2
CYP2B6CYP2B61127,149761 (98%)100.0
CYP2C9CYP2C910250,734632 (98%)72.0
CYP2C19CYP2C19151090,525712 (99%)83.6
CYP2D6CYP2D614214408992 (97%)100.0
CYP3A5CYP3A51131,833643 (98%)49.4
CYP4F2CYP4F21-20,098766 (97%)51.4
DPDDPYD24917,2581211 (98%)40.0
FACT. V LEIDENFACT. V LEIDEN-1 *72,4231679 (97%)41.9
G6PDG6PD1-16,183465 (98%)36.4
HLA-AHLA-A214625423 (71%)100.0
HLA-BHLA-B6787,698308 (78%)62.1
IFNL3IFNL32-1577317 (95%)100.0
IFNL4IFNL42-3543404 (97%)100.0
NUDT15NUDT15339656244 (99%)64.7
RYR-1RYR17-153,8666584 (98%)51.4
SLCO1B1SLCO1B112108,045951 (96%)69.6
TPMTTPMT3326,764346 (97%)52.3
UGT1A1UGT1A11113,052470 (99%)40.3
VKORC1VKORC1135139370 (98%)41.8
* This interaction is aimed at the entire group of drugs classified as oral contraceptives with estrogen.
Table
Table 2. Performance and applicability of available genotyping methods used for PGx. PGx: pharmacogenomics, WES: whole exome sequencing, WGS: whole genome sequencing. NA: not applicable, in this case due to a whole gene/genome coverage and therefore no need for imputation in this region. This table aims to serve as a guide to help select the best technology for the problem at hand. ++ indicates the best score on the parameter, + indicates a good score, - indicates a bad score, -- indicates the worst score on the parameters. Depending on the specific purpose, the weight of the parameters in the selection of the appropriate technology may vary.
Table 2. Performance and applicability of available genotyping methods used for PGx. PGx: pharmacogenomics, WES: whole exome sequencing, WGS: whole genome sequencing. NA: not applicable, in this case due to a whole gene/genome coverage and therefore no need for imputation in this region. This table aims to serve as a guide to help select the best technology for the problem at hand. ++ indicates the best score on the parameter, + indicates a good score, - indicates a bad score, -- indicates the worst score on the parameters. Depending on the specific purpose, the weight of the parameters in the selection of the appropriate technology may vary.
SNV PanelShort-Read SeqLong-Read Seq
PGx PanelWhole Genome PanelPGx PanelWESWGSPGx PanelWGS
Turnaround Time Wetlab * ++++++/----
Haplotype phasingComputational-+/-++/-+++++
Direct-----++++
Imputation -+/-+/-+/-NANANA
Coverage of PGx variation ++/-+++/-++++++
Detection of rare variants † +++++/-++++++
Detection of variants outside the predefined gene/variant panel -----/+-/+++++++
Detection of structural and complex variants ----++/-+++++
Turnaround time data processing * +++++++/----
Costs ‡ [29]Investment+++-----
Running costs per sample+/-+++/-+/--+/--
† For SNV panels, it is assumed that the variants are present in the selected variant panel. * A short turnaround time is indicated by the +. ‡ Lower costs are indicated by the +.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

van der Lee, M.; Kriek, M.; Guchelaar, H.-J.; Swen, J.J. Technologies for Pharmacogenomics: A Review. Genes 2020, 11, 1456. https://doi.org/10.3390/genes11121456

AMA Style

van der Lee M, Kriek M, Guchelaar H-J, Swen JJ. Technologies for Pharmacogenomics: A Review. Genes. 2020; 11(12):1456. https://doi.org/10.3390/genes11121456

Chicago/Turabian Style

van der Lee, Maaike, Marjolein Kriek, Henk-Jan Guchelaar, and Jesse J. Swen. 2020. "Technologies for Pharmacogenomics: A Review" Genes 11, no. 12: 1456. https://doi.org/10.3390/genes11121456

APA Style

van der Lee, M., Kriek, M., Guchelaar, H. -J., & Swen, J. J. (2020). Technologies for Pharmacogenomics: A Review. Genes, 11(12), 1456. https://doi.org/10.3390/genes11121456

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop