SciELO - Scientific Electronic Library Online

 
vol.44 número6Medicina de precisión, algo está cambiando en la reumatologíaDesarrollo de un cuestionario dirigido a conocer el proceso de consentimiento informado en investigación clínica desde la perspectiva del paciente índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados

Journal

Artigo

Indicadores

Links relacionados

  • Em processo de indexaçãoCitado por Google
  • Não possue artigos similaresSimilares em SciELO
  • Em processo de indexaçãoSimilares em Google

Compartilhar


Farmacia Hospitalaria

versão On-line ISSN 2171-8695versão impressa ISSN 1130-6343

Farm Hosp. vol.44 no.6 Toledo Nov./Dez. 2020  Epub 27-Dez-2021

https://dx.doi.org/10.7399/fh.11353 

ORIGINALS

Development and validation of a next-generation sequencing panel for clinical pharmacogenetics

Luis Ramudo-Cela1  2  , Jesús María López-Martí1  , Daniel Colmeiro-Echeberría1  , David de-Uña-Iglesias1  2  , José Luis Santomé-Collazo1  , Lorenzo Monserrat-Iglesias1  2  3 

1Health in Code S. L., Science Department, A Coruña. Spain.

2 Biomedical Research Institute, A Coruña. Spain.

3 Universidade da Coruña, GRINCAR (Cardiovascular Research Group), A Coruña. Spain.

Abstract

The rapid clinical implementation of next generation sequencing techniques is due to its ability to sequence a large number of genetic regions at lower costs than conventional techniques. However, its use in the field of pharmacogenetics is still very limited.

Objective:

Design, development, implementation and validation of a clinical pharmacogenetics next-generation sequencing panel.

Method:

We developed a panel of hybrid capture probes (SureSelect®) for the analysis of the genetic regions of clinical interest collected by literature search and using Illumina HiSeq 1500® sequencing platform. We developed a bioinformatic algorithm for variant annotation, haplotype inference and determination of structural variants in the genes of interest. The results obtained were validated with Coriell® reference material from the pharmacogenetic repositories.

Results:

The developed panel allows the study of a total of 12,794 regions comprised in 389 genes. Validation results showed a sensitivity greater than 99% for single nucleotide variants and small INDELs. Haplotype imputation was consistent with the consensus results in the characterized reference materials. Furthermore, the developed tool was able to correctly identify different types of CYP2D6 copy number variations as well as a wide variety of HLA-B alleles.

Conclusions:

This technology represents an appropriate alternative for its clinical use with advantages over conventional techniques in its throughput and complex gene study capabilities (CYP2D6, HLA-B).

KEYWORDS: High-throughput nucleotide sequencing; Pharmacogenetics; Precision medicine; Computational biology; CYP2D6; HLA-B; DNA copy number variation

Resumen

La rápida implantación clínica de las técnicas de secuenciación masiva en paralelo se debe a su capacidad para secuenciar un gran número de regiones genéticas con un coste menor a las técnicas convencionales. Sin embargo, su uso en el ámbito de la farmacogenética es, todavía, muy escaso.

Objetivo: Diseño, desarrollo, implementación y validación de un panel de secuenciación masiva en paralelo de farmacogenética orientado a la práctica clínica.

Método:

Se desarrolló un panel de sondas de captura híbrida (Sure-

Select®) para el análisis de las regiones genéticas de interés clínico recopiladas mediante búsqueda bibliográfica. Se empleó la plataforma de secuenciación Illumina HiSeq 1500®. Se desarrolló un algoritmo de análisis bioinformático para la anotación de variantes puntuales, inferencia de haplotipos y determinación de variantes estructurales en los genes de interés. Los resultados obtenidos se validaron con materiales de referencia Coriell® de los repositorios de farmacogenética.

Resultados:

El panel desarrollado permite el estudio de un total de 12.794 regiones comprendidas en 389 genes. Los resultados de validación mostraron una sensibilidad superior al 99% para variantes puntuales e inserciones y deleciones pequeñas. La imputación de haplotipos fue coherente con los resultados consenso de los materiales de referencia caracte rizados. Además, la herramienta desarrollada pudo identificar correctamente diferentes tipos de variaciones de número de copias de CYP2D6, así como una gran variedad de alelos de HLA-B.

Conclusiones:

Esta tecnología representa una alternativa adecuada para su empleo asistencial con ventajas frente a las técnicas convencionales en su rendimiento de producción y sus capacidades de estudio de genes complejos (CYP2D6, HLA-B).

PALABRAS CLAVE: Secuenciación de alto rendimiento; Farmacogenética; Medicina de precisión; Biología computacional; CYP2D6; HLA-B; Variaciones del número de copias de ADN

Introduction

Pharmacogenetics is the study of the influence of the genotype on the response to drugs. Clinical pharmacogenetics proposes individualized strategies for drug management based on each patient’s genotype. Its final aim is to improve the healthcare outcomes of pharmacological treatments, improve efficacy, reduce adverse effects, and improve cost-benefit ratios: that is, it is a strategy for the rational use of drugs1.

The incorporation of genetic data into healthcare processes requires high-quality sequencing technology, as well as standardized analysis and interpretation processes. In recent years, there has been a disruptive technological leap due to the development of massively parallel sequencing techniques or next-generation sequencing (NGS). Before this leap, only dozens of variants could be studied, but now tens of thousands can be studied. NGS is increasingly being incorporated into medical diagnostic processes. However, within the field of pharmacogenetics, conventional technologies such as quantitative PCR or primer extension methods such as TaqMan® (Thermofisher®) or MassARRAY® (Sequenom®) continue to be used2,3.

These conventional technologies do not allow the in-depth study of all genetic regions of interest and can only examine subsets of them. Several reviews and comparative studies have addressed pharmacogenetic techniques, finding differences between different laboratories and technologies in the analysis and interpretation of results. It has been suggested that this result is due to the fact that different laboratories analyse different genetic regions without including the entire spectrum of variants of interest4-6. This aspect reduces the clinical value of the studies, creates mistrust in the results, and delays the incorporation of high-quality pharmacogenetics into clinical care practice.

This article presents the design, development, implementation, and validation of a pharmacogenetics NGS sequencing panel oriented towards clinical practice.

Methods

Definition of genetic regions of clinical interest

A literature search was conducted to select regions of interest. A region was considered clinically relevant if it provided information that could modify the therapeutic strategy relative to a given drug (e.g. treatment selection, dosage, patient follow-up). The main literature sources used were the clinical practice guidelines of scientific pharmacogenetic societies, technical files from drug regulatory agencies, and databases related to genetics (see Table 1). The main conventional pharmacogenetic platforms reviewed were Affymetrix DMET® and Agena Bioscience iPLEX®. Some genetic regions, such as CYP2D6 or HLA-B, required a specific design. In order to study copy number variations (CNVs) and structural variants of CYP2D6, probes were designed to capture all coding regions and regions of high homology (CYP2D7 and CYP2D8). For HLA-B, probes were designed against the reference sequences of the IMGT/HLA database7. The approach used was similar to that of Wittig et al.8.

Table 1.  Information sources 

Capture probe design and sequencing process

Hybrid capture probes for these genetic regions were designed using SureSelect® Design Tool (Agilent) software. An automated genomic DNA extraction process (Qiasymphony SP®, Qiagen®) was conducted on the samples. Genomic DNA was prepared using the SureSelect XT® (Agilent Technologies®) protocol in combination with a SureSelect Custom Target Enrichment® probe panel, which selectively captures genomic regions of interest. Sequencing was done using Illumina HiSeq 1500® systems. All procedures were conducted following the manufacturer’s specifications9-11.

Data analysis and interpretation

We used an in-house bioinformatic analysis algorithm, which can annotate single nucleotide polymorphisms (SNPs) variants, Insertion-Deletion (INDELs), and structural variants (i.e. CNVs). CNVs were annotated by comparative analysis of read depth12-14. The variants identified were cross-referenced using an in-house knowledge management system. This system associates specific quality and information values to each variant using structural and SNP data from sources such as the Single Nucleotide Polymorphism Database (dbSNP), 1,000 genomes project15, and the Exome Variant Server (EVS), as well as frequency data of the alleles in different populations, or values of bioinformatic predictors such as SIFT, PMut, and CADD (see Table 1)12-15. Automated haplotype inference based on data from SNPs and small INDELs was performed using tables created in the knowledge management system based on the information sources (see Table 1), especially PharmVar and PharmGKB. CYP2D6 haplotype inference was performed using a module that filters the alignments of CYP2D7 and CYP2D8. The aim of this process is to avoid artefacts in the coverage graphs of the reads in homologous regions16. This module also normalizes the read depth data using a control sample with a known and sequenced CYP2D6 copy number on the same assay plate. The aim of this normalization procedure is to avoid artefacts due to the presence of CNVs in other samples sequenced on the same plate17 HLA-B allele were inferred using a specific module based on the work of Wittig et al.8.

Analytical validation of results using reference materials

The determination of SNPs and small INDELs (< 20 bp) was validated using the Coriell® NA12878 sample. The sample was processed in triplicate and the results were compared to the reference data. This data was the result of integrating several datasets from entire-genome sequencing on different next-generation sequencing platforms18. We calculated the median values for analytical sensitivity (Se), specificity (Sp), and positive (PPV) and negative (NPV) predictive values. This analysis was conducted with three different quality filters using two NGS quality parameters: Quality Score (Qual) and read depth per sample (DP). The sequencing quality filters were defined as high when Qual was more than 49 and DP more than 29, intermediate when Qual was more than 49 and DP more than 14, and low when Qual was more than 0 and DP more than 019.

The determination of pharmacogenetic haplotypes based on SNPs and small INDELs was validated using the Coriell® sample NA12878. The results were compared to the reference data, which had been constructed by integrating various data sets from different pharmacogenetic platforms (GeT-RM project)20,21. These data include the following genes: CYP1A1, CYP1A2, CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2E1, CYP3A4, CYP3A5, CYP4F2, DPYD, GSTM1, GSTP1, GSTT1, NAT1, NAT2, SLC15A2, SLC22A2, SLCO1B1, SLCO2B1, TPMT, UGT1A1, UGT2B7, UGT2B15, UGT2B17, VKORC1. Due to the complexity of HLA-B and CYP2D6, more probands were included: Coriell® samples NA12878, NA02016, NA17254, and NA17281 include different types of structural variants of CYP2D6 and HLA-B alleles of clinical interest (HLA-B*58:01)(8.20.21). This analysis established concordance between the assayed haplotypes and the reference haplotypes.

Results

Genetic regions of interest included in the panel

The panel included a total of 12,794 base pairs distributed among 389 genes. The Supplementary Table 5, Table 6, Table 7, Table 8 shows the size of the region under study for each gene and the average coverage of the region. These regions were classified into three groups: primary, secondary, and candidate genes according to the following internal criteria:

Primary genes (17): CYP2C19, CYP2C9, CYP2D6, CYP3A5, CYP4F2, DPYD, F5, G6PD, HLA-B, IFNL3, RARG, SLC28A3, SLCO1B1, TPMT, UGT1A1, UGT1A6, VKORC1. These genes are described in the main clinical practice guidelines drawn up by the main clinical pharmacogenetic consortia (see Table 1). They have the highest levels of evidence, and customized dosing strategies based on genotype have been published for all of them.

Secondary genes (50): ABCB1, ABCC2, ABCG2, ACE, ADH1A, ADH1B, ADH1C, ADRB1, ADRB2, AHR, ALDH1A1, ALOX5, ASS1, COMT, CPS1, CYP1A1, CYP1A2, CYP2A6, CYP2C8, CYP2E1, CYP2J2, CYP3A4, DRD2, GSTM1, GSTP1, GSTT1, HMGCR, KCNJ11, MTHFR, NAT1, NAT2, NQO1, NR1I2, P2RY1, P2RY12, POLG, PTGIS, SCN1A, SLC15A2, SLC19A1, SLC22A1, SLC22A6, SLC6A4, SLCO1B3, SULT1A1,

TYMS, UGT2B15, UGT2B17, UGT2B7, VDR. Although secondary genes are not included in clinical practice guidelines, we consider that they may be of clinical interest according to the sources used (Table 1). This group includes the FDA lists of pharmacogenetic biomarkers, PharmGKB genes with levels of evidence 2 or higher, Clinical Pharmacogenetics Implementation Consortium (CPIC) priority categories A and B, or “core” genes from the pharmaADME list. Although no strategies for individualised treatment based on genotype have yet been published, we consider that the presence of variants in these genes may be useful for assessing the response to treatment in individual patients following a clinical finding.

Candidate genes: these correspond to genomic regions with lower levels of evidence, but which are still included in one of the main pharmacogenomic platforms commonly used in clinical practice(22.23), or to regions described in the information sources (see Table 1), such as PharmGKB level 3 or 4, or CPIC categories C/D. These studies were set within an investigational framework. The supplementary material shows the full list of genes and regions.

Validation of the determination of SNPs and small INDELs

Table 2 shows the results of the determination of genetic variants of SNPs and INDELs of 20 bp or less. The analytical sensitivity and specificity of the assay were more than 99%. PPV and NPV were also more than 99% for all quality values.

Table 2.  Validation of SNPs and INDELs < 20 bp 

Note: Sample NA12878 was processed in triplicate and the results were compared with reference data generated on different mass sequencing platforms18. The median values for the three runs are shown. The sequencing quality filters were defined as high (Qual > 49 and DP > 29), intermediate (Qual > 49 and DP > 14), and low (Qual > 0 and DP > 0

Validation of the determination of pharmacogenetic alleles (haplotypes)

Table 3 shows the results of the determination of pharmacogenetic haplotypes using the Coriell® NA12878 sample. Reference data is available for 28 genes each with 2 haplotypes (i.e. 56 haplotypes)(20.21). The results were different for four genes: CYP2D6, NAT2, SLCO1B1, and UGT2B17. In the case of genes CYP2D6, NAT2, and SLCO1B1, the differences were due to the fact that the NGS platform examined more genetic regions, and thus improved the detection capabilities of the platforms used by the GeT-RM project(20.21). Regarding CYP2D6, our NGS platform detected the presence of a tandem D6/D7 hybrid type rearrangement, which has been previously described(24.25). In relation to NAT2, the technique was able to detect allele *12. This allele is only included in the Agena Bioscience iPLEX ADME PGx Pro® platform. In the case of SLCO1B1, our technique identified two posible combinations of haplotypes *1A/*15 and *1B/*5. This result is similar to that obtained using the Affimetrix DMET® platform and is due to the inclusion of the suballeles *1A and *1B. Our technique did not detect the presence of haplotype *2 in UGT2B17, which indicates a gene deletion. This requirement was not included in the design specifications. Overall, there was more than 98% concordance in the determination of alleles. Concordance would rise to 100% if the analysis is restricted to genes included in the main clinical practice guidelines (Table 1).

Table 3.  Validation of the determination of pharmacogenetic alleles (haplotypes) 

(1) GeT-RM results for sample NA12878 and CYP2E1 and UGT2B15 genes indicate "no consensus". This is due to the use of two different techniques: Affimetrix DMET® and Agena Bioscience iPLEX®. iPLEX® is unable to detect the CYP2E1*4 or UGT2B15*5 alleles. Therefore, the table includes the results from Affimetrix DMET®. (2) The consensus result for GSTP1 is ambiguous: two possible allele combinations are established for a given combination of variants.

Validation of the results of the identification of CYP2D6 structural variants

Our technique was based on the comparative study of read depth, and was able to correctly identify the structural variants present in the four selected control samples. Figure 1 shows the results of the coverage graph analysis for CYP2D6. Sample NA17254 with the CYP2D6*4/*41 genotype presents two alleles of CYP2D6. This is the control sample against which the read depth of CYP2D6 is normalised. NA17281 with the CYP2D6*5/*9 genotype presents a deletion of one of the CYP2D6 alleles, and the normalized read depth decreased by approximately 50%. NA02016 with the CYP2D6*2x2/*17 genotype presents a duplication of one of the alleles, and had a 50% increase in read depth. Sample NA12878 with genotype*3/*4+68 presents a tandem D6/D7 hybrid, and had a 50% increase in read depth in part of the CYP2D6 region and part of the CYP2D7 region.

Figure 1.  Read depth and normalised read depth graphs. 

Validation of the results of the identification of HLA variants

The imputation technique developed to infer alleles in the HLA-B region was able to identify the alleles present in the three control samples. Table 4 shows the established consensus alleles and the results obtained using our technique.

Table 4.  Results of HLA-B allele identification 

Discussion

This article describes the development and analytical validation of a next generation sequencing platform for clinical pharmacogenetics, and presents the results obtained with the bioinformatics tools used to interpret the data.

Our panel includes a total of 12,794 base pairs distributed among 389 genes. The panel is much larger than that of conventional platforms: Affymetrix DMET® (1936 markers), Agena Bioscience® formerly (Sequenom®), iPLEX ADME PGx Pro® (270 markers), Illumina VeraCode ADME® (184 markers), and Roche Amplichip® (36 CYP2D6 and CYP2C19 alleles)22,23.

The results from the control samples were consistent with the reference results of the GeT-RM project. Differences in the results were due to the higher capacity of our platform, except in the case of UGT2B7. For this gene, deletion-type structural variants have been described that would have required a specific design on our platform. However, this was not done because it has not been included in any of the main pharmacogenetic clinical practice guidelines.

It is difficult to study certain genomic regions of great pharmacogenetic interest, such as CYP2D6 and HLA-B. CYP2D6 metabolises approximately 25% of commercialized drugs26 and its genotype is key to individualizing antidepressant treatments such as paroxetine27, antipsychotics such as aripiprazole28, tamoxifen29, antiemetics such as ondansetron30, or analgesics such as codeine31, among others. It is difficult to study because of the existence of two regions with very high homology (i.e. CYP2D7 and CYP2D8), and the existence of CNV-like variants and pseudogene rearrangement (conversion-type and tandem hybrids with functional alleles). Other tools for the bioinformatic analysis of CYP2D6 have obtained good results in the analysis of SNP-type variants or small INDELs, but are of limited use in the determination of structural variants, especially hybrids16,17,24,32,33. This issue was addressed as described in the materials and methods section: our NGS technique filters the reads from the homologous regions and normalizes the read depth against a known control sample. The NGS plat form appropriately identified different types of structural variants and CNVs analysed in the reference materials and provided graphs with the read depth showing the behaviour of coverage throughout the entire genomic region (see Figure 1).

The Human Leukocyte Antigen (HLA) region is the most densely polymorphic region of the genome. Its relevance in pharmacogenetics stems from its association with hypersensibility to immunomediated drugs34. Conventional typing techniques can only analyse a few haplotypes, while general procedures using NGS technologies do not yield good results due to variability and sequence homology between genes and pseudogenes35. The current NGS version incorporates a specific probe design and previously validated analysis software that enables the comprehensive study of HLA-B.

Compared to conventional platforms, the main limitations of the NGS platform are assay time, the need for qualified technicians, development and commissioning costs, and higher costs of data processing and storage on servers36. It takes approximately 10 days to obtain the results after a sample has been received (i.e. including sample preparation, sequencing, and bioinformatic analysis). These processes require technical staff with specific training and bioinformatics personnel to fine-tune the technique and analyse the sequencing data.

Although our technique and bioinformatics data process identified rare variants and other previously-undescribed variants, these results cannot be extrapolated to clinical care settings since there are no data in the literature for their interpretation. These rare variants account for over 80% of variations in the genes of interest in pharmacogenetics: however, conventional platforms are not oriented to their analysis37,38. These variants may be incorporated into clinical care practice when scientific articles on genotype-phenotype correlation data are published, consensus on interpretation is built, and protocols become available39.

Finally, the evaluation and validation of this platform in prospective clinical studies would demonstrate the clinical value of this platform.

Supplement table 1.  Genomic regions, size, and average coverage by gene 

Supplement table 1 (cont.).  Genomic regions, size, and average coverage by gene 

Supplement table 1 (cont.).  Genomic regions, size, and average coverage by gene 

Supplement table 1 (cont.).  Genomic regions, size, and average coverage by gene 

Contribution to the scientific literature

This article describes the development, implementation, and validation of a next-generation sequencing panel for clinical pharmacogenetics.

Acknowledgments

We would like to thank the staff of Health in Code S. L

Bibliography

Roden DM, McLeod HL, Relling MV, Williams MS, Mensah GA, Peterson JF, et al. Pharmacogenomics. Lancet (Internet). 2019;6736(19):1-12. Available at: https://linkinghub.elsevier.com/retrieve/pii/S0140673619312760Links ]

Han SM, Park J, Lee JH, Lee SS, Kim H, Han H, et al. Targeted Next-Generation Sequencing for Comprehensive Genetic Profiling of Pharmacogenes. Clin Pharmacol Ther. 2017;101(3):396-405. [ Links ]

Ji Y, Si Y, McMillin GA, Lyon E. Clinical pharmacogenomics testing in the era of next generation sequencing: challenges and opportunities for precision medicine. Expert Rev Mol Diagn (Internet). 2018;18(5):411-21. Available at: https://doi.org/10.1080/14737159.2018.1461561Links ]

Bank P, Caudle K, Swen J, Gammal R, Whirl-Carrillo M, Klein T, et al. Comparison of the Guidelines of the Clinical Pharmacogenetics Implementation Consortium and the Dutch Pharmacogenetics Working Group. Clin Pharmacol Ther (Internet). 2017;4457. Available at: http://onlinelibrary.wiley.com/doi/10.1002/cpt.762/fullLinks ]

Bousman CA, Jaksa P, Pantelis C. Systematic evaluation of commercial pharmacogenetic testing in psychiatry. Pharmacogenet Genomics (Internet). 2017;27(11):387-93. Available at: https://www.ncbi.nlm.nih.gov/pubmed/28777243 [ Links ]

Kalman LV, Agúndez JAG, Appell ML, Black JL, Bell GC, Boukouvala S, et al. Pharmacogenetic Allele Nomenclature: International Workgroup Recommendations for Test Result Reporting. Clin Pharmacol Ther (Internet). 2016;99(2):172-85. Available at: https://www.ncbi.nlm.nih.gov/pubmed/26479518 [ Links ]

Robinson J, Barker DJ, Georgiou X, Cooper MA, Flicek P, Marsh SGE. IPD-IMGT/ HLA Database. Nucleic Acids Res (Internet). 2019 Oct 31. Available at: https://academic.oup.com/nar/advance-article/doi/10.1093/nar/gkz950/5610347Links ]

Wittig M, Anmarkrud JA, Kassens JC, Koch S, Forster M, Ellinghaus E, et al. Development of a high-resolution NGS-based HLA-typing and analysis pipeline. Nucleic Acids Res (Internet). 2015;43(11):e70-e70. Available at: https://academic.oup.com/nar/article-lookup/doi/10.1093/nar/gkv184Links ]

Agilent Technologies. SureSelectXT Target Enrichment System for Illumina Paired-End Multiplexed Sequencing Library (Internet). 2019. Available at: https://www.agilent.com/cs/library/usermanuals/Public/G7530-90000.pdfLinks ]

Qiagen. Manual de instrucciones de uso de QIAsymphony® DSP DNA (Internet). 2015. Available at: https://www.qiagen.com/fr/resources/download.aspx?id=8bc88dad-4140-467e-a1f0-e390fd193865&lang=es-ESLinks ]

Illumina. HiSeq 1500 System Guide (15035788 v03) (Internet). 2019. Available at: https://emea.support.illumina.com/downloads/hiseq_1500_user_guide_15035788.htmlLinks ]

Ochoa JP, Sabater-Molina M, García-Pinilla JM, Mogensen J, Restrepo-Córdoba A ,Palomino-Doza J, et al. Formin Homology 2 Domain Containing 3 (FHOD3) Is a Genetic Basis for Hypertrophic Cardiomyopathy. J Am Coll Cardiol (Internet). 2018;72(20): 2457-67. Available at: http://www.ncbi.nlm.nih.gov/pubmed/30442288 [ Links ]

Ortiz-Genga MF, Cuenca S, Dal Ferro M, Zorio E, Salgado-Aranda R, Climent V, et al. Truncating FLNC Mutations Are Associated With High-Risk Dilated and Arrhythmogenic Cardiomyopathies. J Am Coll Cardiol (Internet). 2016;68(22):2440-51. Available at: http://www.ncbi.nlm.nih.gov/pubmed/27908349 [ Links ]

Trujillo-Quintero JP, Gutiérrez-Agulló M, Ochoa JP, Martínez-Martínez JG, de Uña D, García-Fernández A. Familial Brugada Syndrome Associated With a Complete Deletion of the SCN5A and SCN10A Genes. Rev Esp Cardiol (Engl Ed) (Internet). 2019;72(2):176-8. Available at: http://www.ncbi.nlm.nih.gov/pubmed/29650450 [ Links ]

The 1000 Genomes Project Consortium. An integrated map of genetic variation from 1,092 human genomes. Nature (Internet). 2012;491(7422):56-65. Available at: http://www.nature.com/articles/nature11632Links ]

Numanagić I, Malikić S, Pratt VM, Skaar TC, Flockhart DA, Sahinalp SC. Cypiripi: exact genotyping of CYP2D6 using high-throughput sequencing data. Bioinformatics (Internet). 2015;31(12):i27-34. Available at: https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/btv232Links ]

Lee S, Wheeler MM, Patterson K, McGee S, Dalton R, Woodahl EL, et al. Stargazer: a software tool for calling star alleles from next-generation sequencing data using CYP2D6 as a model. Genet Med (Internet). 2019;21(2):361-72. Available at: http://dx.doi.org/10.1038/s41436-018-0054-0 [ Links ]

Zook JM, Chapman B, Wang J, Mittelman D, Hofmann O, Hide W, et al. Integrating human sequence data sets provides a resource of benchmark SNP and indel genotype calls. Nat Biotechnol (Internet). 2014;32(3):246-51. Available at: http://www.ncbi.nlm.nih.gov/pubmed/24531798 [ Links ]

BCFtools - Variant Calling (Internet) (accessed 06/08/2020). Available at: https://samtools.github.io/bcftools/howtos/variant-calling.htmlLinks ]

Pratt VM, Everts RE, Aggarwal P, Beyer BN, Broeckel U, Epstein-Baak R, et al. Characterization of 137 Genomic DNA Reference Materials for 28 Pharmacogenetic Genes: A GeT-RM Collaborative Project. J Mol Diagnostics (Internet). 2016;18(1):109-23. Available at: http://dx.doi.org/10.1016/j.jmoldx.2015.08.005 [ Links ]

Pratt VM, Zehnbauer B, Wilson JA, Baak R, Babic N, Bettinotti M, et al. Characterization of 107 Genomic DNA Reference Materials for CYP2D6, CYP2C19, CYP2C9, VKORC1, and UGT1A1. J Mol Diagnostics (Internet). 2010;12(6):835-46. Available at: https://linkinghub.elsevier.com/retrieve/pii/S1525157810601341Links ]

Arbitrio M, Di Martino MT, Scionti F, Agapito G, Hiram Guzzi P, Cannataro M, et al. DMETTM (Drug Metabolism Enzymes and Transporters): a Pharmacogenomic platform for precision medicine. Oncotarget (Internet). 2016;5(0). Available at: http://www.oncotarget.com/abstract/9927Links ]

Agena Bioscience: iPLEX® PGx Pro Panel Flyer (Internet) (accessed 06/04/2020). Available at: https://agenabio.com/resources/product-literature/Links ]

Twist GP, Gaedigk A, Miller NA, Farrow EG, Willig LK, Dinwiddie DL, et al. Constellation: a tool for rapid, automated phenotype assignment of a highly polymorphic pharmacogene, CYP2D6, from whole-genome sequences. npj Genomic Med (Internet). 2016;1(August 2015):15007. Available at: http://www.nature.com/articles/npjgenmed20157Links ]

Gaedigk A, Turner A, Everts RE, Scott SA, Aggarwal P, Broeckel U, et al. Characterization of Reference Materials for Genetic Testing of CYP2D6 Alleles: A GeT-RM Collaborative Project. J Mol Diagnostics (Internet). 2019;21(6):1034-52. Available at: https://doi.org/10.1016/j.jmoldx.2019.06.007 [ Links ]

Owen RP, Sangkuhl K, Klein TE, Altman RB. Cytochrome P450 2D6. Pharmacogenet Genomics (Internet). 2009;19(7):559-62. Available at: http://www.ncbi.nlm.nih.gov/pubmed/19512959 [ Links ]

Hicks JK, Bishop JR, Sangkuhl K, Muller DJ, Ji Y, Leckband SG, et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) guideline for CYP2D6 and CYP2C19 genotypes and dosing of selective serotonin reuptake inhibitors. Clin Pharmacol Ther. 2015;98(2):127-34. [ Links ]

Swen JJ, Nijenhuis M, de Boer A, Grandia L, Maitland-van der Zee AH, Mulder H, et al. Pharmacogenetics: from bench to byte--an update of guidelines. Clin Pharmacol Ther (Internet). 2011;89(5):662-73. Available at: http://www.ncbi.nlm.nih.gov/pubmed/21412232 [ Links ]

Drögemöller BI, Wright GEB, Shih J, Monzon JG, Gelmon KA, Ross CJD, et al. CYP2D6 as a treatment decision aid for ER-positive non-metastatic breast cancer patients: a systematic review with accompanying clinical practice guidelines. Breast Cancer Res Treat (Internet). 2019;173(3):521-32. Available at: http://dx.doi.org/10.1007/s10549-018-5027-0 [ Links ]

Bell GC, Caudle KE, Whirl-Carrillo M, Gordon RJ, Hikino K, Prows CA, et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) guideline for CYP2D6 genotype and use of ondansetron and tropisetron. Clin Pharmacol Ther (Internet). 2017;102(2):213-8. Available at: http://www.pharmgkb.org/page/cpic/Links ]

Madadi P, Amstutz U, Rieder M, Ito S, Fung V, Hwang S, et al. Clinical practice guideline: CYP2D6 genotyping for safe and efficacious codeine therapy. J Popul Ther Clin Pharmacol. 2013;20(3):369-96. [ Links ]

Klein TE, Ritchie MD. PharmCAT: A Pharmacogenomics Clinical Annotation Tool. Clin Pharmacol Ther (Internet). 2018;104(1):19-22. Available at: http://doi.wiley.com/10.1002/cpt.928 [ Links ]

Numanagić I, Malikić S, Ford M, Qin X, Toji L, Radovich M, et al. Allelic decomposition and exact genotyping of highly polymorphic and structurally variant genes. Nat Commun (Internet). 2018;9(1):828. Available at: http://www.nature.com/articles/s41467-018-03273-1Links ]

Barbarino JM, Kroetz DL, Klein TE, Altman RB. PharmGKB summary: very important pharmacogene information for human leukocyte antigen B. Pharmacogenet Genomics (Internet). 2015;25(4):205-21. Available at: http://www.ncbi.nlm.nih.gov/pubmed/25647431 [ Links ]

Hosomichi K, Shiina T, Tajima A, Inoue I. The impact of next-generation sequencing technologies on HLA research. J Hum Genet (Internet). 2015;60(11):665-73. Available at: http://www.ncbi.nlm.nih.gov/pubmed/26311539 http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=PMC4660052 [ Links ]

Roy S, Coldren C, Karunamurthy A, Kip NS, Klee EW, Lincoln SE, et al. Standards and Guidelines for Validating Next-Generation Sequencing Bioinformatics Pipelines. J Mol Diagnostics (Internet). 2018;20(1):4-27. Available at: https://linkinghub.elsevier.com/retrieve/pii/S1525157817303732Links ]

Hovelson DH, Xue Z, Zawistowski M, Ehm MG, Harris EC, Stocker SL, et al. Characterization of ADME gene variation in 21 populations by exome sequencing. Pharmacogenet Genomics (Internet). 2017;27(3):89-100. Available at: http://content.wkhealth.com/linkback/openurl?sid=WKPTLP:landingpage&an=01213011-900000000-99273Links ]

Kozyra M, Ingelman-Sundberg M, Lauschke VM. Rare genetic variants in celular transporters, metabolic enzymes, and nuclear receptors can be important determinants of interindividual differences in drug response. Genet Med. 2017;19(1):20-9. [ Links ]

Lauschke VM, Ingelman-Sundberg M. How to Consider Rare Genetic Variants in Personalized Drug Therapy. Clin Pharmacol Ther (Internet). 2018;103(5):745-8. Available at: http://dx.doi.org/10.1002/cpt.976 [ Links ]

Received: April 12, 2020; Accepted: June 13, 2020

Author of correspondence. Email: luis.ramudo@healthincode.com (Luis Ramudo Cela).

Conflicts of interest No conflict of interests

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License