Whole-exome sequencing (WES) has been widely used for analysis of human genetic diseases, but its value for the pharmacogenomic profiling of individuals is not well studied. Initially, we performed an in-depth evaluation of the accuracy of WES variant calling in the pharmacogenes CYP2D6 and CYP2C19 by comparison with MiSeq(®) amplicon sequencing data (n = 36). This analysis revealed that the concordance rate between WES and MiSeq(®) was high, achieving 99.60% for variants that were called without exceeding the truth-sensitivity threshold (99%), defined during variant quality score recalibration (VQSR). Beyond this threshold, the proportion of discordant calls increased markedly. Subsequently, we expanded our findings beyond CYP2D6 and CYP2C19 to include more genes genotyped by the iPLEX(®) ADME PGx Panel in the subset of twelve samples. WES performed well, agreeing with the genotyping panel in approximately 99% of the selected pass-filter variant calls. Overall, our results have demonstrated WES to be a promising approach for pharmacogenomic profiling, with an estimated error rate of lower than 1%. Quality filters, particularly VQSR, are important for reducing the number of false variants. Future studies may benefit from examining the role of WES in the clinical setting for guiding drug therapy.
Next generation DNA sequencing and analysis of amplicons spanning the pharmacogene CYP2D6 suggested that the Nextera transposase used for fragmenting and providing sequencing priming sites displayed a targeting bias. This manifested as dramatically lower sequencing coverage at sites in the amplicon that appeared likely to form G-quadruplex structures. Since secondary DNA structures such as G-quadruplexes are abundant in the human genome, and are known to interact with many other proteins, we further investigated these sites of low coverage. Our investigation revealed that G-quadruplex structures are formed in vitro within the CYP2D6 pharmacogene at these sites, and G-quadruplexes can interact with the hyperactive Tn5 transposase (EZ-Tn5) with high affinity. These findings indicate that secondary DNA structures such as G-quadruplexes may represent preferential transposon integration sites and provide additional evidence for the role of G-quadruplex structures in transposition or viral integration processes.