Multi-ethnic cytochrome-P450 copy number profiling: novel pharmacogenetic alleles and mechanism of copy number variation formation

Unveiling novel genetic variants in 370 challenging medically relevant genes using the long read sequencing data of 41 samples from 19 global populations

BACKGROUND

A large number of challenging medically relevant genes (CMRGs) are situated in complex or highly repetitive regions of the human genome, hindering comprehensive characterization of genetic variants using next-generation sequencing technologies. In this study, we employed long-read sequencing technology, extensively utilized in studying complex genomic regions, to characterize genetic alterations, including short variants (single nucleotide variants and short insertions and deletions) and copy number variations, in 370 CMRGs across 41 individuals from 19 global populations.

RESULTS

Our analysis revealed high levels of genetic variants in CMRGs, with 68.73% exhibiting copy number variations and 65.20% containing short variants that may disrupt protein function across individuals. Such variants can influence pharmacogenomics, genetic disease susceptibility, and other clinical outcomes. We observed significant differences in CMRG variation across populations, with individuals of African ancestry harboring the highest number of copy number variants and short variants compared to samples from other continents. Notably, 15.79% to 33.96% of short variants were exclusively detectable through long-read sequencing. While the T2T-CHM13 reference genome significantly improved the assembly of CMRG regions, thereby facilitating variant detection in these regions, some regions still lacked resolution.

CONCLUSION

Our results provide an important reference for future clinical and pharmacogenetic studies, highlighting the need for a comprehensive representation of global genetic diversity in the reference genome and improved variant calling techniques to fully resolve medically relevant genes.

Characterization of CYP2B6 and CYP2A6 Pharmacogenetic Variation in Sub-Saharan African Populations

Genetic variation in CYP2B6 and CYP2A6 is known to impact interindividual response to antiretrovirals, nicotine, and bupropion, among other drugs. However, the full catalogue of clinically relevant pharmacogenetic variants in these genes is yet to be established, especially across African populations. This study therefore aimed to characterize the star allele (haplotype) distribution in CYP2B6 and CYP2A6 across diverse and understudied sub-Saharan African (SSA) populations. We called star alleles from 961 high-depth full genomes using StellarPGx, Aldy, and PyPGx. In addition, we performed CYP2B6 and CYP2A6 star allele frequency comparisons between SSA and other global biogeographical groups represented in the new 1000 Genomes Project high-coverage dataset (n = 2,000). This study presents frequency information for star alleles in CYP2B6 (e.g., *6 and *18; frequency of 21-47% and 2-19%, respectively) and CYP2A6 (e.g., *4, *9, and *17; frequency of 0-6%, 3-10%, and 6-20%, respectively), and predicted phenotypes (for CYP2B6), across various African populations. In addition, 50 potentially novel African-ancestry star alleles were computationally predicted by StellarPGx in CYP2B6 and CYP2A6 combined. For each of these genes, over 4% of the study participants had predicted novel star alleles. Three novel star alleles in CYP2A6 (*54, *55, and *56) and CYP2B6 apiece, and several suballeles were further validated via targeted Single-Molecule Real-Time resequencing. Our findings are important for informing the design of comprehensive pharmacogenetic testing platforms, and are highly relevant for personalized medicine strategies, especially relating to antiretroviral medication and smoking cessation treatment in Africa and the African diaspora. More broadly, this study highlights the importance of sampling diverse African ethnolinguistic groups for accurate characterization of the pharmacogene variation landscape across the continent.

Copy number variations of cytochrome P450 genes in Kinh Vietnamese

Background

The cytochrome P450 (CYP450) family is well known as a major group of drug metabolizing enzymes. The polymorphism of CYP450 genes is the main factor having an impact on the interindividual difference in drug response, including drug efficacy and drug safety. The single nucleotide polymorphism (SNPs) of Vietnamese Kinh has been widely studied, but information about the copy number variations (CNVs) of other CYP450 genes is still unknown.

Objective

To identify the CNV variability of CYP450 in 154 healthy unrelated Kinh Vietnamese, except eCYP2D6, which was previously reported.

Methods

Multiplex Ligation-Dependent Probe Amplification (MLPA) was applied for determination of copy number of 10 CYP450 genes. Later, PCR or quantitative PCR (qPCR) was used to confirm the detected CNVs in randomly chosen subjects.

Results

Of the 154 subjects, along with CYP2D6, 4 other CYP450 genes showed CNVs including duplications (CYP1B1), deletions (CYP2A6 and CYP2C9), and both duplications and deletions (CYP2E1). Among these, CYP2A6 exhibited the greatest frequency of CNVs compared with other CYP450, in which CYP2A6Del accounted for 11%. Meanwhile, allele CYP2E1Del showed the lowest frequency with only 0.3%.

Conclusions

The present study provides new insight into CYP450 CNVs in the Kinh Vietnamese cohort. Our data have contributed to genetic profiling of CYP450 CNVs in Vietnam, which would be helpful for facilitating implementation of pharmacogenetics in drug dosing adjustment in Vietnam.

Clinical pharmacogenetic analysis in 5,001 individuals with diagnostic Exome Sequencing data

Exome sequencing is utilized in routine clinical genetic diagnosis. The technical robustness of repurposing large-scale next-generation sequencing data for pharmacogenetics has been demonstrated, supporting the implementation of preemptive pharmacogenetic strategies based on adding clinical pharmacogenetic interpretation to exomes. However, a comprehensive study analyzing all actionable pharmacogenetic alleles contained in international guidelines and applied to diagnostic exome data has not been performed. Here, we carried out a systematic analysis based on 5001 Spanish or Latin American individuals with diagnostic exome data, either Whole Exome Sequencing (80%), or the so-called Clinical Exome Sequencing (20%) (60 Mb and 17 Mb, respectively), to provide with global and gene-specific clinical pharmacogenetic utility data. 788 pharmacogenetic alleles, distributed through 19 genes included in Clinical Pharmacogenetics Implementation Consortium guidelines were analyzed. We established that Whole Exome and Clinical Exome Sequencing performed similarly, and 280 alleles in 11 genes (CACNA1S, CYP2B6, CYP2C9, CYP4F2, DPYD, G6PD, NUDT15, RYR1, SLCO1B1, TPMT, and UGT1A1) could be used to inform of pharmacogenetic phenotypes that change drug prescription. Each individual carried in average 2.2 alleles and overall 95% (n = 4646) of the cohort could be informed of at least one actionable pharmacogenetic phenotype. Differences in variant allele frequency were observed among the populations studied and the corresponding gnomAD population for 7.9% of the variants. In addition, in the 11 selected genes we uncovered 197 novel variants, among which 27 were loss-of-function. In conclusion, we provide with the landscape of actionable pharmacogenetic information contained in diagnostic exomes, that can be used preemptively in the clinics.

Genophenotypic Factors and Pharmacogenomics in Adverse Drug Reactions

Adverse drug reactions (ADRs) rank as one of the top 10 leading causes of death and illness in developed countries. ADRs show differential features depending upon genotype, age, sex, race, pathology, drug category, route of administration, and drug–drug interactions. Pharmacogenomics (PGx) provides the physician effective clues for optimizing drug efficacy and safety in major problems of health such as cardiovascular disease and associated disorders, cancer and brain disorders. Important aspects to be considered are also the impact of immunopharmacogenomics in cutaneous ADRs as well as the influence of genomic factors associated with COVID-19 and vaccination strategies. Major limitations for the routine use of PGx procedures for ADRs prevention are the lack of education and training in physicians and pharmacists, poor characterization of drug-related PGx, unspecific biomarkers of drug efficacy and toxicity, cost-effectiveness, administrative problems in health organizations, and insufficient regulation for the generalized use of PGx in the clinical setting. The implementation of PGx requires: (i) education of physicians and all other parties involved in the use and benefits of PGx; (ii) prospective studies to demonstrate the benefits of PGx genotyping; (iii) standardization of PGx procedures and development of clinical guidelines; (iv) NGS and microarrays to cover genes with high PGx potential; and (v) new regulations for PGx-related drug development and PGx drug labelling.

Applying Next-Generation Sequencing Platforms for Pharmacogenomic Testing in Clinical Practice

Pharmacogenomics (PGx) studies the use of genetic data to optimize drug therapy. Numerous clinical centers have commenced implementing pharmacogenetic tests in clinical routines. Next-generation sequencing (NGS) technologies are emerging as a more comprehensive and time- and cost-effective approach in PGx. This review presents the main considerations for applying NGS in guiding drug treatment in clinical practice. It discusses both the advantages and the challenges of implementing NGS-based tests in PGx. Moreover, the limitations of each NGS platform are revealed, and the solutions for setting up and management of these technologies in clinical practice are addressed.

DNA copy number and structural variation (CNV) contributions to adult and childhood obesity

In recent years, obesity has reached epidemic proportions globally and has become a major public health concern. The development of obesity is likely caused by several behavioral, environmental, and genetic factors. Genomic variability among individuals is largely due to copy number variations (CNVs). Recent genome-wide association studies (GWAS) have successfully identified many loci containing CNV related to obesity. These obesity-related CNVs are informative to the diagnosis and treatment of genomic diseases. A more comprehensive classification of CNVs may provide the basis for determining how genomic diversity impacts the mechanisms of expression for obesity in children and adults of a variety of genders and ethnicities. In this review, we summarize current knowledge on the relationship between obesity and the CNV of several genomic regions, with an emphasis on genes at the following loci: 11q11, 1p21.1, 10q11.22, 10q26.3, 16q12.2, 16p12.3, and 4q25.

Han Chinese specific cytochrome P450 polymorphisms and their impact on the metabolism of anti-hypertensive drugs with adrenoreceptor blocking properties

Introduction: Cytochrome P450 (CYP) is a monooxygenase superfamily mediating the elimination of anti-hypertensive drugs. Polymorphisms of CYP would lead to differential drug efficacy. Building relationships between genotype and phenotype will benefit individual medical treatment of hypertension.Areas covered: The review systematically summarizes the polymorphisms of four CYPs (CYP2C9, CYP2C19, CYP2D6, and CYP3A4) concentrated distributed in the Han Chinese population. Moreover, the activity of variants on metabolizing anti-hypertensive drugs are reviewed, especially drugs with adrenoceptor blocking properties, as well as their clinical relevancies.Expert opinion: The polymorphisms of CYP can cause stratification in drug exposure of antihypertensive drugs. Although the clinical relevance has been built partially, the translational medicine still lacks reliable data support. Furthermore, the studies have demonstrated that even the same CYP variant will exhibit different catalytic capability for different drugs, which is another obstacle to hinder its application. With the deepening of multiomics research and structural biology, nucleotide polymorphisms can be combined with transcriptome, proteome, metabolome and molecular structure analyses to study the susceptibility to hypertension and drug efficacy. A complete data chain would be further estabolished by combining studies of pharmacokinetics-pharmacodynamics, which can effectively promote the precise application of anti-hypertensive drugs.

How Can Drug Metabolism and Transporter Genetics Inform Psychotropic Prescribing?

Many genetic variants in drug metabolizing enzymes and transporters have been shown to be relevant for treating psychiatric disorders. Associations are strong enough to feature on drug labels and for prescribing guidelines based on such data. A range of commercial tests are available; however, there is variability in included genetic variants, methodology, and interpretation. We herein provide relevant background for understanding clinical associations with specific variants, other factors that are relevant to consider when interpreting such data (such as age, gender, drug-drug interactions), and summarize the data relevant to clinical utility of pharmacogenetic testing in psychiatry and the available prescribing guidelines. We also highlight areas for future research focus in this field.

Identification and population genetic analyses of copy number variations in six domestic goat breeds and Bezoar ibexes using next-generation sequencing

Background

Copy number variations (CNVs) are a major form of genetic variations and are involved in animal domestication and genetic adaptation to local environments. We investigated CNVs in the domestic goat (Capra hircus) using Illumina short-read sequencing data, by comparing our lab data for 38 goats from three Chinese breeds (Chengdu Brown, Jintang Black, and Tibetan Cashmere) to public data for 26 individuals from three other breeds (two Moroccan and one Chinese) and 21samples from Bezoar ibexes.

Results

We obtained a total of 2394 CNV regions (CNVRs) by merging 208,649 high-confidence CNVs, which spanned ~ 267 Mb of total length and accounted for 10.80% of the goat autosomal genome. Functional analyses showed that 2322 genes overlapping with the CNVRs were significantly enriched in 57 functional GO terms and KEGG pathways, most related to the nervous system, metabolic process, and reproduction system. Clustering patterns of all 85 samples generated separately from duplications and deletions were generally consistent with the results from SNPs, agreeing with the geographical origins of these goats. Based on genome-wide F_ST at each CNV locus, some genes overlapping with the highly divergent CNVs between domestic and wild goats were mainly enriched for several immunity-related pathways, whereas the genes overlapping with the highly differentiated CNVs between highland and lowland goats were mainly related to vitamin and lipid metabolism. Remarkably, a 507-bp deletion at ~ 14 kb downstream of FGF5 on chromosome 6 showed highly divergent (F_ST = 0.973) between the highland and lowland goats. Together with an enhancer activity of this sequence shown previously, the function of this duplication in regulating fiber growth deserved to be further investigated in detail.

Conclusion

We generated a comprehensive map of CNVs in goats. Many genetically differentiated CNVs among various goat populations might be associated with the population characteristics of domestic goat breeds.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07267-6.

Development and Analytical Validation of a 29 Gene Clinical Pharmacogenetic Genotyping Panel: Multi-Ethnic Allele and Copy Number Variant Detection

To develop a novel pharmacogenetic genotyping panel, a multidisciplinary team evaluated available evidence and selected 29 genes implicated in interindividual drug response variability, including 130 sequence variants and additional copy number variants (CNVs). Of the 29 genes, 11 had guidelines published by the Clinical Pharmacogenetics Implementation Consortium. Targeted genotyping and CNV interrogation were accomplished by multiplex single‐base extension using the MassARRAY platform (Agena Biosciences) and multiplex ligation‐dependent probe amplification (MRC Holland), respectively. Analytical validation of the panel was accomplished by a strategic combination of > 500 independent tests performed on 170 unique reference material DNA samples, which included sequence variant and CNV accuracy, reproducibility, and specimen (blood, saliva, and buccal swab) controls. Among the accuracy controls were 32 samples from the 1000 Genomes Project that were selected based on their enrichment of sequence variants included in the pharmacogenetic panel (VarCover.org). Coupled with publicly available samples from the Genetic Testing Reference Materials Coordination Program (GeT‐RM), accuracy validation material was available for the majority (77%) of interrogated sequence variants (100% with average allele frequencies > 0.1%), as well as additional structural alleles with unique copy number signatures (e.g., CYP2D6*5, *13, *36, *68; CYP2B6*29; and CYP2C19*36). Accuracy and reproducibility for both genotyping and copy number were > 99.9%, indicating that the optimized panel platforms were precise and robust. Importantly, multi‐ethnic allele frequencies of the interrogated variants indicate that the vast majority of the general population carries at least one of these clinically relevant pharmacogenetic variants, supporting the implementation of this panel for pharmacogenetic research and/or clinical implementation programs.

Structural variation at the CYP2C locus: Characterization of deletion and duplication alleles

The human CYP2C locus harbors the polymorphic CYP2C18, CYP2C19, CYP2C9, and CYP2C8 genes, and of these, CYP2C19 and CYP2C9 are directly involved in the metabolism of ~15% of all medications. All variant CYP2C19 and CYP2C9 star (*) allele haplotypes currently cataloged by the Pharmacogene Variation (PharmVar) Consortium are defined by sequence variants. To determine if structural variation also occurs at the CYP2C locus, the 10q23.33 region was interrogated across deidentified clinical chromosomal microarray (CMA) data from 20,642 patients tested at two academic medical centers. Fourteen copy number variants that affected the coding region of CYP2C genes were detected in the clinical CMA cohorts, which ranged in size from 39.2 to 1,043.3 kb. Selected deletions and duplications were confirmed by MLPA or ddPCR. Analysis of the clinical CMA and an additional 78,839 cases from the Database of Genomic Variants (DGV) and ClinGen (total n = 99,481) indicated that the carrier frequency of a CYP2C structural variant is ~1 in 1,000, with ~1 in 2,000 being a CYP2C19 full gene or partial-gene deletion carrier, designated by PharmVar as CYP2C19*36 and *37, respectively. Although these structural variants are rare in the general population, their detection will likely improve metabolizer phenotype prediction when interrogated for research and/or clinical testing.

Copy number variation profiling in pharmacogenes using panel-based exome resequencing and correlation to human liver expression

Structural variants including copy number variations (CNV) have gained widespread attention, especially in pharmacogenomics but for several genes functional relevance and clinical evidence are still lacking. Detection of CNVs in next-generation sequencing data is challenging but offers widespread applications. We developed a cohort-based CNV detection workflow to extract CNVs from read counts of targeted NGS of 340 genes involved in absorption, distribution, metabolism and excretion (ADME) of drugs. We applied our method to 150 human liver tissue samples and correlated identified CNVs to mRNA expression levels. In total, we identified 445 deletions (73%) and 167 duplications (27%) in 36 pharmacogenes including all well-known CNVs of CYPs, GSTs, SULTs, UGTs, numerous described rare CNVs of CYP2E1, SLC16A3 or UGT2B15 as well as novel observations, e.g., for SLC22A12, SLC22A17 and GPS2 (G Protein Pathway Suppressor 2). We were able to fine-map complex CNVs of CYP2A6 and CYP2D6 with exon resolution. Correlation analysis confirmed known expression patterns for common CNVs and suggested an influence on expression variability for some rare CNVs. Our straightforward CNV detection workflow can be easily applied to any NGS coverage data and helped to analyze CNVs in an ADME-NGS panel of 340 pharmacogenes to improve genotype-phenotype correlations.

The importance of both CYP2C19 and CYP2B6 germline variations in cyclophosphamide pharmacokinetics and clinical outcomes

Cyclophosphamide is an alkylating agent used in the treatment of solid and haematological malignancies and as an immunosuppressive agent. As a prodrug, it is dependent on bioactivation to the active phosphoramide mustard metabolite to elicit its therapeutic effect. This focused review will highlight the evidence for the role of germline pharmacogenetic variation in both plasma pharmacokinetics and clinical outcomes. There is a substantial indication from 13 pharmacokinetic and 17 therapeutic outcome studies, in contexts as diverse as haematological malignancy, breast cancer, systemic lupus erythematosus and myeloablation, that pharmacogenetic variation in both CYP2C19 and CYP2B6 influence the bioactivation of cyclophosphamide. An additional role for pharmacogenetic variation in ALDH1A1 has also been reported. Future studies should comprehensively assess these 3 pharmacogenes and undertake appropriate statistical analysis of gene-gene interactions to confirm these findings and may allow personalised treatment regimens.

Copy number variation profiling in pharmacogenetics CYP-450 and GST genes in Colombian population

BACKGROUND

Copy Number variation (CNVs) in genes related to drug absorption, distribution, metabolism and excretion (ADME) are relevant in the interindividual variability of drug response. Studies of the CNVs in ADME genes in Latin America population are lacking. The objective of the study was to identify the genetic variability of CNVs in CYP-450 and GST genes in a subgroup of individuals of Colombian origin.

METHODS

Genomic DNA was isolated from 123 healthy individuals from a Colombian population. Multiplex Ligation-Dependent Probe Amplification (MLPA) was performed for the identification of CNVs in 40 genomic regions of 11 CYP-450 and 3 GST genes. The genetic variability, allelic and genotypic frequencies were analyzed.

RESULTS

We found that 13 out of 14 genes had CNVs: 5 (35.7%) exhibited deletions and duplications, while 8 (57.1%) presented either deletions or duplications.. 33.3% of individuals carried deletions and duplications while 49.6% had a unique type of CNV (deletion or duplication). The allelic frequencies of the CYP and GST genes were 0 to 47.6% (allele null), 0 to 17.5% (duplicated alleles) and 37 to 100% (normal alleles).

CONCLUSIONS

Our results describe, for the first time, the genomic profile of CNVs in a subgroup of Colombian population in GST and CYP-450 genes. GST genes indicated greater genetic variability than CYP-450 genes. The data obtained contributes to the knowledge of genetic profiles in Latin American subgroups. Although the clinical relevance of CNVs has not been fully established, it is a valuable source of pharmacogenetic variability data with potential involvement in the response to medications.

Pharmacogenomics, CYP2D6, and Tamoxifen: A Survey of the Reasons Sustaining European Clinical Practice Paradigms

Tamoxifen is a drug that is often used in the clinical management of breast cancer. CYP2D6 is a key metabolizing enzyme that is involved in the conversion of tamoxifen to its active drug metabolites. CYP2D6 has several alleles that metabolize tamoxifen and other drugs at different rates that can alter therapeutic impact, a characteristic that renders it one of the most studied enzymes in the field of pharmacogenetics. Background and objectives: Portugal has no implemented measures based on pharmacogenomics analysis prior to therapy that might function as a cultural sample control when analyzing the individual and economic factors present in clinical practice paradigms. Therefore, we aim to investigate the impact of CYP2D6 genotyping of the tamoxifen metabolizing enzymes in the clinical management of breast cancer patients. Materials and Methods: Qualitative/quantitative studies regarding the impact of pharmacogenomics in breast cancer; personal interviews in different Portuguese laboratories within hospital setting using a survey. Analysis of data through interviews to management board and/or decision makers from major oncological centers. Results: Reasons for common adoption of pharmacogenomics practice are contradictory and based both in economic factors and cultural/clinical bias. Conclusions: This research study identifies specific cultural and/or clinical bias that act as obstacles to pharmacogenomic implementation and proposes viable courses of action that might bring about change in cultural/medical habits.

Acute Respiratory Distress Syndrome Associated With Clopidogrel in a Young Male Patient

Background: Clopidogrel is a commonly prescribed antiplatelet drug in patients with stents and histories of arterial vascular disease. It generally has a favorable side effect profile with increasing bleeding risk as the main concern as an adverse event.

Case Presentation: A 19-year-old previously healthy male presented to the neurological intensive care unit with a subarachnoid hemorrhage requiring a flow diverting stent to secure the aneurysm. The patient was stable for 2 weeks and had no changes to management or medication within 48 h of antiplatelet therapy. Within hours of first-time dosing of clopidogrel, the patient experienced a syncopal episode and dyspnea. He was difficult to arouse and using accessory muscles to breath with an otherwise benign exam. He was hypoxic with bibasilar crackles requiring bilevel positive airway pressure (BiPap). Imaging showed bilateral pulmonary edema and he was diagnosed with moderate acute respiratory distress syndrome (ARDS). Infectious, cardiogenic, and contrast-induced ARDS were ruled out. Upon cessation of clopidogrel, his pulmonary function and mental status improved.

Conclusions: This is the first reported case of a young and immunocompetent patient's severe pulmonary edema leading to acute respiratory distress syndrome in association with first- time dosing of clopidogrel.

CYP2D6 Pharmacogenetics Testing and Post-Cesarean Section Pain Scores-a Preliminary Study

OBJECTIVE

Prospective observational study to analyze CYP2D6 pharmacogenetics in 55 Portuguese adult parturients undergoing elective cesarean section and to investigate the association between CYP2D6 alleles and pain score.

METHODS

DNA was extracted from peripheral blood by standard methods. Genetic analysis included allelic discrimination (CYP2D6*1, *2, *3, *4, *5, *6, *10, *17, and *41) and copy number determination with TaqMan probes by real-time polymerase chain reaction (PCR). Allele duplications were confirmed (long PCR and PCR-restriction fragment length polymorphism). Theoretical metabolic profiles prediction was based on genetic data and activity scores. Association was investigated between genotypes and predicted phenotypes with pain scores. Statistical analysis was performed by using a χ2 test, and significance was set at P < 0.05.

RESULTS

The percentage of poor, intermediate, extensive, and ultrarapid metabolizers found were 9%, 38%, 46%, and 7%, respectively. The results reveal a positive association between alleles *4, *10, and pain.

CONCLUSIONS

A positive association was found between predicted reduced or null activity of CYP2D6 and increased pain. It can be hypothesized that if CYP2D6 activity is reduced, tyramine metabolism will decrease, resulting in reduced formation of endogenous dopamine. Consequently, activation of the signal transduction pathways that controls pain and analgesic effect may be reduced, leading to an increase in pain. Therefore, we would recommend CYP2D6 genotyping to anticipate the needs for analgesia, which will help to adjust opioid dose and maximize clinical efficacy while reducing side effects.

Integrated CYP2D6 interrogation for multiethnic copy number and tandem allele detection

AIM

To comprehensively interrogate CYP2D6 by integrating genotyping, copy number analysis and novel strategies to identify CYP2D6*36 and characterize CYP2D6 duplications.

METHODS

Genotyping of 16 CYP2D6 alleles, multiplex ligation-dependent probe amplification (MLPA) and CYP2D6*36 and duplication allele-specific genotyping were performed on 427 African-American, Asian, Caucasian, Hispanic, and Ashkenazi Jewish individuals.

RESULTS

A novel PCR strategy determined that almost half of all CYP2D6*10 (100C>T) alleles are actually *36 (isolated or in tandem with *10) and all identified duplication alleles were characterized. Integrated results from all testing platforms enabled the refinement of genotype frequencies across all studied populations.

CONCLUSION

The polymorphic CYP2D6 gene requires comprehensive interrogation to characterize allelic variation across ethnicities, which was enabled in this study by integrating multiplexed genotyping, MLPA copy number analysis, novel PCR strategies and duplication allele-specific genotyping.

Genetic Basis of Delayed Hypersensitivity Reactions to Drugs in Jewish and Arab Populations

Genetic variation can affect drug pharmacokinetics and pharmacodynamics and contribute to variability between individuals in response to medications. Specifically, differences in allele frequencies among individuals and ethnic groups have been associated with variation in their propensity to develop drug hypersensitivity reactions (HSRs). This article reviews the current knowledge on the genetic background of HSRs and its relevance to Jewish and Arab populations. The focus is on human leukocyte antigen (HLA) alleles and haplotypes as predictive markers of HSRs ("immunopharmacogenetics"), but other genes and alleles are described as well. Also discussed is the translation of the pharmacogenetic information to practice recommendations.

Comprehensive overview of the pharmacogenetic diversity in Ashkenazi Jews

BACKGROUND

Adverse drug reactions are a major concern in drug development and clinical therapy. Genetic polymorphisms in genes involved in drug metabolism and transport are major determinants of treatment efficacy and adverse reactions, and constitute important biomarkers for drug dosing, efficacy and safety. Importantly, human populations and subgroups differ substantially in their pharmacogenetic variability profiles, with important consequences for personalised medicine strategies and precision public health approaches. Despite their long migration history, Ashkenazi Jews constitute a rather isolated population with a unique genetic signature that is distinctly different from other populations.

OBJECTIVE

To provide a comprehensive overview of the pharmacogenetic profile in Ashkenazim.

METHODS

We analysed next-generation sequencing data from 5076 Ashkenazim individuals and used sequence data from 117 425 non-Jewish individuals as reference.

RESULTS

We derived frequencies of 164 alleles in 17 clinically relevant pharmacogenes and derived profiles of putative functional consequences, providing the most comprehensive data set of Jewish pharmacogenetic diversity published to date. Furthermore, we detected 127 variants with an aggregated frequency of 20.7% that were specifically found in Ashkenazim, of which 55 variants were putatively deleterious (aggregated frequency of 9.4%).

CONCLUSION

The revealed pattern of pharmacogenetic variability in Ashkenazi Jews is distinctly different from other populations and is expected to translate into unique functional consequences, especially for the metabolism of CYP2A6, CYP2C9, NAT2 and VKORC1 substrates. We anticipate that the presented data will serve as a powerful resource for the guidance of pharmacogenetic treatment decisions and the optimisation of population-specific genotyping strategies in the Ashkenazi diaspora.

Multi-ethnic SULT1A1 copy number profiling with multiplex ligation-dependent probe amplification

AIM

To develop a SULT1A1 multiplex ligation-dependent probe amplification assay and to investigate multi-ethnic copy number variant frequencies.

METHODS

A novel multiplex ligation-dependent probe amplification assay was developed and tested on 472 African-American, Asian, Caucasian, Hispanic and Ashkenazi Jewish individuals.

RESULTS

The frequencies of atypical total copy number (i.e., greater or less than two) were 38.7% for Hispanics, 38.9% for Ashkenazi Jewish, 43.2% for Caucasians, 53.6% for Asians and 64.1% for African-Americans. Heterozygous SULT1A1 deletion carriers (slow sulfators) were most common among Caucasians (8.4%), whereas African-Americans had the highest frequencies of three or more copies (rapid sulfators; 60.9%).

CONCLUSION

Different ethnic and racial populations have varying degrees of SULT1A1-mediated sulfation activity, which warrants further research and that may have utility for drug response prediction among SULT1A1-metabolized medications.

Pharmacogenetics and interstitial lung disease

PURPOSE OF REVIEW

Interstitial lung disease (ILD) is comprised of a heterogeneous group of disorders with highly variable natural histories and response to therapies. Pharmacogenetics focuses on the variability in drug response because of the presence of genetic factors that influence drug metabolism or disease activity. In this article, we review relevant drug-specific and disease-specific polymorphisms that may influence therapeutic response, and then highlight a recently identified drug-gene interaction in patients with idiopathic pulmonary fibrosis (IPF).

RECENT FINDINGS

The emergence of high-throughput genomic technology has allowed for identification of gene polymorphisms associated with susceptibility to specific disease states, including IPF and several connective tissue diseases known to cause ILD. IPF risk loci span a diverse group of genes, while most associated with connective tissue disease are critical to immune signaling. A recent pharmacogenetic analysis of patients enrolled in an IPF clinical trial identified a variant within TOLLIP to be associated with differential response to N-acetylcysteine therapy.

SUMMARY

Though few pharmacogenetic investigations have been conducted in patients with ILD to date, ample opportunities for pharmacogenetic exploration exist in this patient population. Such exploration will advance our understanding of specific ILDs and help usher in an era of personalized medicine.

Novel copy-number variations in pharmacogenes contribute to interindividual differences in drug pharmacokinetics

PurposeVariability in pharmacokinetics and drug response is shaped by single-nucleotide variants (SNVs) as well as copy-number variants (CNVs) in genes with importance for drug absorption, distribution, metabolism, and excretion (ADME). While SNVs have been extensively studied, a systematic assessment of the CNV landscape in ADME genes is lacking.MethodsWe integrated data from 2,504 whole genomes from the 1000 Genomes Project and 59,898 exomes from the Exome Aggregation Consortium to identify CNVs in 208 relevant pharmacogenes.ResultsWe describe novel exonic deletions and duplications in 201 (97%) of the pharmacogenes analyzed. The deletions are population-specific and frequencies range from singletons up to 1%, accounting for >5% of all loss-of-function alleles in up to 42% of the genes studied. We experimentally confirmed novel deletions in CYP2C19, CYP4F2, and SLCO1B3 by Sanger sequencing and validated their allelic frequencies in selected populations.ConclusionCNVs are an additional source of pharmacogenetic variability with important implications for drug response and personalized therapy. This, together with the important contribution of rare alleles to the variability of pharmacogenes, emphasizes the necessity of comprehensive next-generation sequencing-based genotype identification for an accurate prediction of the genetic variability of drug pharmacokinetics.

Characterization of ADME gene variation in 21 populations by exome sequencing

OBJECTIVE

Proteins involving absorption, distribution, metabolism, and excretion (ADME) play a critical role in drug pharmacokinetics. The type and frequency of genetic variation in the ADME genes differ among populations. The aim of this study was to systematically investigate common and rare ADME coding variation in diverse ethnic populations by exome sequencing.

MATERIALS AND METHODS

Data derived from commercial exome capture arrays and next-generation sequencing were used to characterize coding variation in 298 ADME genes in 251 Northeast Asians and 1181 individuals from the 1000 Genomes Project.

RESULTS

Approximately 75% of the ADME coding sequence was captured at high quality across the joint samples harboring more than 8000 variants, with 49% of individuals carrying at least one 'knockout' allele. ADME genes carried 50% more nonsynonymous variation than non-ADME genes (P=8.2×10) and showed significantly greater levels of population differentiation (P=7.6×10). Out of the 2135 variants identified that were predicted to be deleterious, 633 were not on commercially available ADME or general-purpose genotyping arrays. Forty deleterious variants within important ADME genes, with frequencies of at least 2% in at least one population, were identified as candidates for future pharmacogenetic studies.

CONCLUSION

Exome sequencing was effective in accurately genotyping most ADME variants important for pharmacogenetic research, in addition to identifying rare or potentially de novo coding variants that may be clinically meaningful. Furthermore, as a class, ADME genes are more variable and less sensitive to purifying selection than non-ADME genes.

Institutional profile: translational pharmacogenomics at the Icahn School of Medicine at Mount Sinai

For almost 50 years, the Icahn School of Medicine at Mount Sinai has continually invested in genetics and genomics, facilitating a healthy ecosystem that provides widespread support for the ongoing programs in translational pharmacogenomics. These programs can be broadly cataloged into discovery, education, clinical implementation and testing, which are collaboratively accomplished by multiple departments, institutes, laboratories, companies and colleagues. Focus areas have included drug response association studies and allele discovery, multiethnic pharmacogenomics, personalized genotyping and survey-based education programs, pre-emptive clinical testing implementation and novel assay development. This overview summarizes the current state of translational pharmacogenomics at Mount Sinai, including a future outlook on the forthcoming expansions in overall support, research and clinical programs, genomic technology infrastructure and the participating faculty.

Differential Regulation of CYP3A4 and CYP3A5 and its Implication in Drug Discovery

BACKGROUND

Cancer cells use several mechanisms to resist the cytotoxic effects of drugs, resulting in tumor progression and invasion. One such mechanism capitalizes on the body's natural defense against xenobiotics by increasing the rate of xenobiotic efflux and metabolic inactivation. Xenobiotic metabolism typically involves conversion of parent molecules to more soluble and easily excreted derivatives in reactions catalyzed by Phase I and Phase II drug metabolizing enzymes.

METHODS

We performed a structured search of peer-reviewed literature on P450 (CYP) 3A, with a focus on CYP3A4 and CYP3A5.

RESULTS

Recent reports indicate that components of the xenobiotic response system are upregulated in some diseases, including many cancers. Such components include the pregnane X receptor (PXR), CYP3A4 and CYP3A5 enzymes. The CYP3A enzymes are a subset of the numerous enzymes that are transcriptionally activated following the interaction of PXR and many ligands.

CONCLUSION

Intense research is ongoing to understand the functional ramifications of aberrant expression of these components in diseased states with the goal of designing novel drugs that can selectively target them.

High frequency of CYP2D6 ultrarapid metabolizers in Spain: controversy about their misclassification in worldwide population studies

A high frequency (7-10%) of CYP2D6 ultrarapid metabolizers estimated from the genotype (gUMs) has been claimed to exist among Spaniards and Southern Europeans. However, methodological aspects such as the inclusion of individuals carrying non-active multiplied alleles as gUMs may have led to an overestimation. Thus, this study aimed to analyze the gUM frequency (considering only those carrying more than two active genes) in 805 Spanish healthy volunteers studied for CYP2D6*2, *3, *4, *5, *6, *10, *17, *35, *41, and multiplications. Second, all worldwide studies reporting gUM frequencies were reviewed in order to evaluate potential misclassifications. The gUM frequency in this Spanish population was 5.34%, but increased to 8.3% if all individuals with CYP2D6 multiplications were classified as gUMs without considering the activity of the multiplied alleles. Moreover, among all reviewed worldwide studies only 55.6% precisely determined whether the multiplied alleles were active. Present results suggest that the evaluation of CYP2D6 ultrarapid metabolism should be standarized, and that the frequency of gUMs should be reconsidered in Spaniards and globally.

High frequency of CYP2D6 ultrarapid metabolizer genotypes in an Ashkenazi Jewish population from Argentina

A twofold higher frequency of CYP2D6 ultrarapid metabolizers (estimated from genotype: gUMs) was reported among Ashkenazi Jews (AJ) living in New York (USA) than in other North American Caucasians, which might be important to guide the prescription for CYP2D6 substrates in AJ communities around the world. This study was aimed to determine whether the high frequency of CYP2D6 gUMs described in AJ from USA was replicated in AJ from Argentina when compared with other multiethnic admixture Argentines (GA). The frequency of the most common allelic variants and of CYP2D6 gUMs (>2 active genes) and poor metabolizers (0 active genes, gPMs) was also compared among the studied Argentine populations. CYP2D6 genotyping was performed in 173 AJ and 246 GA DNA samples of unrelated donors from the metropolitan area of Buenos Aires. CYP2D6 alleles (*2, *3, *4, *5, *6, *10, *17, *35, *41 and multiple copies), genotypes and functional phenotype frequencies were determined. The frequencies of gUMs and gPMs in AJ from Argentina were 11.5% and 5.2%, respectively, whereas in GA, the frequencies of gUM and gPMs were 6.5% and 4.9%, respectively. Comparisons between AJ and GA showed that gUMs frequencies were twofold higher (P<0.05) in AJ than GA. CYP2D6*35 allele was more frequent in GA than AJ, whereas CYP2D6*41 and *1xN were more frequent in AJ than in GA (P<0.05). This study supports the previously reported high frequency of gUMs on another Ashkenazi population in New York. The present findings also support the interethnic variability of CYP2D6 genetic polymorphism in the overall Argentine population.

Pharmacogenetics of drug-metabolizing enzymes in Italian populations

Drug-metabolizing enzymes play a major role in the biotransformation and subsequent elimination of most drugs and xenobiotics from the body. Both phase I and phase II enzymes are highly polymorphic. Inter-individual differences in genes coding for drug-metabolizing enzymes are important for understanding variability in drug response and for individualization of drug prescription. The prevalence of genetic polymorphisms in drug metabolism varies widely with ethnicity, and marked differences in the distribution of allelic variants of genes encoding drug-metabolizing enzymes have been documented in populations of different racial origin. This review aimed to summarize the available studies on genetic polymorphisms associated with drug metabolism conducted in Italian populations and to compare the frequency of the various metabolizer phenotypes and most common variant alleles (and resulting genotypes) with corresponding values from other populations. Notably, published data are not extensive, and most studies were performed on relatively low numbers of individuals. In general, the frequency of polymorphisms in the cytochrome P450 (CYP) genes as well as in the investigated phase II enzymes in the Italian population was similar to values reported for other Caucasian populations. However, the prevalence of CYP2D6 gene duplication among Italians was found to be very high, confirming the higher frequency of CYP2D6 ultrarapid metabolizers in the Mediterranean area compared to Northern Europe. It is worth noting that a geographic gradient in the flavin-containing monooxygenase 3 polymorphism distribution was also seen, the Italian population showing higher similarity to other Mediterranean populations than to North Europeans.

Long-Read Single Molecule Real-Time Full Gene Sequencing of Cytochrome P450-2D6

The cytochrome P450-2D6 (CYP2D6) enzyme metabolizes ∼25% of common medications, yet homologous pseudogenes and copy number variants (CNVs) make interrogating the polymorphic CYP2D6 gene with short-read sequencing challenging. Therefore, we developed a novel long-read, full gene CYP2D6 single molecule real-time (SMRT) sequencing method using the Pacific Biosciences platform. Long-range PCR and CYP2D6 SMRT sequencing of 10 previously genotyped controls identified expected star (*) alleles, but also enabled suballele resolution, diplotype refinement, and discovery of novel alleles. Coupled with an optimized variant-calling pipeline, CYP2D6 SMRT sequencing was highly reproducible as triplicate intra- and inter-run nonreference genotype results were completely concordant. Importantly, targeted SMRT sequencing of upstream and downstream CYP2D6 gene copies characterized the duplicated allele in 15 control samples with CYP2D6 CNVs. The utility of CYP2D6 SMRT sequencing was further underscored by identifying the diplotypes of 14 samples with discordant or unclear CYP2D6 configurations from previous targeted genotyping, which again included suballele resolution, duplicated allele characterization, and discovery of a novel allele and tandem arrangement. Taken together, long-read CYP2D6 SMRT sequencing is an innovative, reproducible, and validated method for full-gene characterization, duplication allele-specific analysis, and novel allele discovery, which will likely improve CYP2D6 metabolizer phenotype prediction for both research and clinical testing applications.

Gene copy number variation analysis reveals dosage-insensitive expression of CYP2E1

Gene copy number variants (CNVs) of CYP2E1 have been described but not functionally characterized. Here we investigated effects of CNVs on hepatic and lymphoblastoid CYP2E1 expression. Using available single-nuleotide polymorphism microarray data and quantitative PCR, CYP2E1 gene duplication and deletion carriers were identified. CYP2E1 mRNA, protein and enzyme activity (chlorzoxazone-6-hydroxylation) phenotypes of CYP2E1 were not associated with gene copy number. Analysis of gene expression in lymphoblastoid cell lines in relation to CNV confirmed this finding in an extrahepatic tissue and for other ethnicities. Further analyses identified a linked haplotype cluster with possible influence on gene expression. In summary, our data suggest a homeostatic, gene dosage-insensitive regulation of CYP2E1 expression by unknown gene dosage compensation mechanisms. This is in striking contrast to well-known structural variations of CYP2A6 and CYP2D6 that have a strong impact on expression and activity. These findings are important in the context of pharmacogenetic prediction.

Prenatal antidepressant exposure associated with CYP2E1 DNA methylation change in neonates

Some but not all neonates are affected by prenatal exposure to serotonin reuptake inhibitor antidepressants (SRI) and maternal mood disturbances. Distinguishing the impact of these 2 exposures is challenging and raises critical questions about whether pharmacological, genetic, or epigenetic factors can explain the spectrum of reported outcomes. Using unbiased DNA methylation array measurements followed by a detailed candidate gene approach, we examined whether prenatal SRI exposure was associated with neonatal DNA methylation changes and whether such changes were associated with differences in birth outcomes. Prenatal SRI exposure was first associated with increased DNA methylation status primarily at CYP2E1(β(Non-exposed) = 0.06, β(SRI-exposed) = 0.30, FDR = 0); however, this finding could not be distinguished from the potential impact of prenatal maternal depressed mood. Then, using pyrosequencing of CYP2E1 regulatory regions in an expanded cohort, higher DNA methylation status--both the mean across 16 CpG sites (P < 0.01) and at each specific CpG site (P < 0.05)--was associated with exposure to lower 3rd trimester maternal depressed mood symptoms only in the SRI-exposed neonates, indicating a maternal mood x SRI exposure interaction. In addition, higher DNA methylation levels at CpG2 (P = 0.04), CpG9 (P = 0.04) and CpG10 (P = 0.02), in the interrogated CYP2E1 region, were associated with increased birth weight independently of prenatal maternal mood, SRI drug exposure, or gestational age at birth. Prenatal SRI antidepressant exposure and maternal depressed mood were associated with altered neonatal CYP2E1 DNA methylation status, which, in turn, appeared to be associated with birth weight.

Single tube genotyping of CYP2A6 gene deletion based on copy number determination by quantitative real-time PCR

The CYP2A6*4 allele, characterized as the whole deletion of this gene, is closely associated with nicotine dependence, cancer susceptibility, and drug responsiveness. It has long been a significant challenge for pharmacogenetics scientists to develop a reliable method to detect this molecular variant due to its high homology with its homologous genes CYP2A6 and CYP2A3 in the clinical setting. Here, we introduce a quantitative real-time PCR assay that specifically amplifies CYP2A6 by designing a specific set of primers and the probe, which effectively prevent the amplification of the CYP2A7 and CYP2A13 alleles. CYP2A6 gene copy numbers were normalized to albumin (ALB) which was co-amplified simultaneously in a single-tube duplex reaction and at a setting as the internal reference gene. The established assay was validated with a selection of previously genotyped DNA samples, which harbored none, one or two CYP2A6 gene copies. The results were in complete concordance with previously published data and no overlap between the three groups was observed. Further analysis of a cohort of 120 samples revealed high specificity and sensitivity of this assay as demonstrated by the agreement of determined gene copy numbers in all of the cases. In conclusion, this novel assay allows reliable and sensitive detection of the CYP2A6 gene deletion, which will be useful for pharmacogenetics studies and routine clinical settings.

Pharmacogenetics of cytochrome P450 2B6 (CYP2B6): advances on polymorphisms, mechanisms, and clinical relevance

Cytochrome P450 2B6 (CYP2B6) belongs to the minor drug metabolizing P450s in human liver. Expression is highly variable both between individuals and within individuals, owing to non-genetic factors, genetic polymorphisms, inducibility, and irreversible inhibition by many compounds. Drugs metabolized mainly by CYP2B6 include artemisinin, bupropion, cyclophosphamide, efavirenz, ketamine, and methadone. CYP2B6 is one of the most polymorphic CYP genes in humans and variants have been shown to affect transcriptional regulation, splicing, mRNA and protein expression, and catalytic activity. Some variants appear to affect several functional levels simultaneously, thus, combined in haplotypes, leading to complex interactions between substrate-dependent and -independent mechanisms. The most common functionally deficient allele is CYP2B6*6 [Q172H, K262R], which occurs at frequencies of 15 to over 60% in different populations. The allele leads to lower expression in liver due to erroneous splicing. Recent investigations suggest that the amino acid changes contribute complex substrate-dependent effects at the activity level, although data from recombinant systems used by different researchers are not well in agreement with each other. Another important variant, CYP2B6*18 [I328T], occurs predominantly in Africans (4-12%) and does not express functional protein. A large number of uncharacterized variants are currently emerging from different ethnicities in the course of the 1000 Genomes Project. The CYP2B6 polymorphism is clinically relevant for HIV-infected patients treated with the reverse transcriptase inhibitor efavirenz, but it is increasingly being recognized for other drug substrates. This review summarizes recent advances on the functional and clinical significance of CYP2B6 and its genetic polymorphism, with particular emphasis on the comparison of kinetic data obtained with different substrates for variants expressed in different recombinant expression systems.