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

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.

Interrogating Pharmacogenetics Using Next-Generation Sequencing

BACKGROUND

Pharmacogenetics or pharmacogenomics (PGx) is the study of the role of inherited or acquired sequence change in drug response. With the rapid evolution of molecular techniques, bioinformatic tools, and increased throughput of functional genomic studies, the discovery of PGx associations and clinical implementation of PGx test results have now moved beyond a handful variants in single pharmacogenes and multi-gene panels that interrogate a few pharmacogenes to whole-exome and whole-genome scales. Although some laboratories have adopted next-generation sequencing (NGS) as a testing platform for PGx and other molecular tests, most clinical laboratories that offer PGx tests still use targeted genotyping approaches.

CONTENT

This article discusses primarily the technical considerations for clinical laboratories to develop NGS-based PGx tests including whole-genome and whole-exome sequencing analyses and highlights the challenges and opportunities in test design, content selection, bioinformatic pipeline for PGx allele and diplotype assignment, rare variant classification, reporting, and briefly touches a few additional areas that are important for successful clinical implementation of PGx results.

SUMMARY

The accelerated speed of technology development associated with continuous cost reduction and enhanced ability to interrogate complex genome regions makes it inevitable for most, if not all, clinical laboratories to transition PGx testing to an NGS-based platform in the near future. It is important for laboratories and relevant professional societies to recognize both the potential and limitations of NGS-based PGx profiling, and to work together to develop a standard and consistent practice to maximize the variant or allele detection rate and utility of PGx testing.

CYP2D6 genotyping and the clinical impact on outcomes in breast cancer tamoxifen-treated patients

Aims: To report the distribution of allele frequencies of CYP2D6 gene and to evaluate their influence on the clinical outcomes of a group of breast cancer patients receiving adjuvant tamoxifen treatment from Uruguay. Patients & methods: 199 samples were genotyped through real-time polymerase chain reaction assays. Metabolization profiles were inferred from the genotypes. Correlations were evaluated using Pearson's χ2 test. Results: Phenotype frequencies were 0.65 normal (NM), 0.30 intermediate (IM) and 0.05 poor metabolizers (PM). Similar clinical outcomes between NM and (PM + IM) patient groups (odds ratio = 1.011, 95% CI = 0.2703-3.7826; p = 0.987) were found. Conclusion: CYP2D6 allele frequencies were analyzed for the first time in a cohort from Uruguay. Results did not support any impact of CYP2D6 gene polymorphisms on clinical outcomes.

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.

Molecular Variants in Genes related to the Response to Ocular Hypotensive Drugs in an Afro-Colombian Population

Aims:

This study aimed to conduct an exploratory analysis of the pharmacogenomic variants involved in ocular hypotensive drugs to understand the individual differential response in an Afro-descendant population.

Background:

Glaucoma is the leading cause of irreversible blindness worldwide. The pharmacologic treatment available consists of lowering intraocular pressure by administering topical drugs. In Asian and Caucasian people, pharmacogenomic variants associated with the efficacy of these treatments have been identified. However, in Afro-descendant populations, there is a profound gap in this knowledge.

Objective:

This study identified the pharmacogenomic variants related to ocular hypotensive efficacy treatment in Afro-descendant individuals from the Archipelago of San Andres and Providence, Colombia.

Methods:

An analysis of whole-exome sequencings (WES), functional annotation, and clinical significance was performed for pharmacogenomic variants reported in PharmGKB databases; in turn, an in silico available prediction analysis was carried out for the novel variants.

Results:

We identified six out of 18 non-synonymous variants with a clinical annotation in PharmGKB. Five were classified as level three evidence for the hypotensive drugs; rs1801252 and rs1801253 in the ADRB1 gene and rs1042714 in the ADRB2 gene. These pharmacogenomic variants have been involved in a lack of efficacy of topical beta-blockers and higher systolic and diastolic pressure under treatment with ophthalmic timolol drug. The rs1045642 in the ABCB1 gene was associated with greater efficacy of treatments with latanoprost drug. Also, we found the haplotypes *17 for CYP2D6 and *10 for CYP2C19; both related to reducing the enzyme activity to timolol drug metabolization. In addition, we observed 50 novel potentially actionable variants; 36 synonymous, two insertion variants that caused frameshift mutations, and 12 non-synonymous, where five were predicted to be pathogenic based on several pathogenicity predictions.

Conclusion:

Our results suggested that the pharmacogenomic variants were found to decrease the ocular hypotensive efficacy treatment in a Colombian Afro-descendant population and revealed a significant proportion of novel variants with a potential to influence drug response.

Simple and Robust Detection of CYP2D6 Gene Deletions and Duplications Using CYP2D8P as Reference

Genotyping of the CYP2D6 gene is the most commonly applied pharmacogenetic test globally. Significant economic interests have led to the development of a plurality of assays, available for almost any genotyping platform or DNA detection chemistry. Of all the genetic variants, copy number variations are particular difficult to detect by polymerase chain reaction. Here, we present two simple novel approaches for the identification of samples carrying either deletions or duplications of the CYP2D6 gene; by relative quantification using a singleplex 5′nuclease real-time PCR assay, and by high-resolution melting of PCR products. These methods make use of universal primers, targeting both the CYP2D6 and the reference gene CYP2D8P, which is necessary for the analysis. The assays were validated against a reference method using a large set of samples. The singleplex nature of the 5′nuclease real-time PCR ensures that the primers anneal with equal affinity to both the sequence of the CYP2D6 and the reference gene. This facilitates robust identification of gene deletions and duplications based on the cycle threshold value. In contrast, the high-resolution melting assay is an end-point PCR, where the identification relies on variations between the amount of product generated from each of the two genes.

Whole-Exome Sequencing in Patients Affected by Stevens-Johnson Syndrome and Toxic Epidermal Necrolysis Reveals New Variants Potentially Contributing to the Phenotype

Background

Adverse drug reactions (ADRs) are frequent occurring events that can essentially be defined as harmful or unpleasant symptoms secondary to the use of a medicinal product. ADRs involve a wide spectrum of clinical manifestations ranging from minor itching and rash to life-threatening reactions. Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN) are rare ADRs. SJS-TEN may be considered a polygenic pathology due to additive/epistatic effects caused by sequence variants in numerous genes. Next-generation sequencing (NGS) represents a potentially interesting exploration tool in such scenario as it facilitates the simultaneous analysis of large genomic regions and genes at affordable cost.

Methods

The present study has involved using whole-exome sequencing (WES) for the first time on SJS-TEN patients. It involved robust and innovative multistep bioinformatics analysis focusing on 313 candidate genes potentially participating in the disease’s aetiology, specific drugs’ metabolism and gene regulation.

Results

We identified combinations of frequently occurring and rare variants that may contribute to the disease’s pathogenesis. Depending on the specific drug being taken, different variants (and alleles) in NAT2, CYP2D8, CYP2B6, ABCC2, UGT2B7 and TCF3 were identified as coherent candidates representing potential future markers for SJS-TEN.

Conclusion

The present study proposed and has described (for the first time) a large-scale genomic analysis of patients affected by SJS-TEN. The genes and variants identified represent relevant candidates potentially participating in the disease’s pathogenesis. Corroborating that proposed by others, we found that complex combinations of frequently occurring and rare variants participating in particular drug metabolism molecular cascades could be associated with the phenotype. TCF3 TF may be considered a coherent candidate for SJS-TEN that should be analysed in new cohorts of patients having ADRs.

Potential of whole-genome sequencing-based pharmacogenetic profiling

Pharmacogenetics represents a major driver of precision medicine, promising individualized drug selection and dosing. Traditionally, pharmacogenetic profiling has been performed using targeted genotyping that focuses on common/known variants. Recently, whole-genome sequencing (WGS) is emerging as a more comprehensive short-read next-generation sequencing approach, enabling both gene diagnostics and pharmacogenetic profiling, including rare/novel variants, in a single assay. Using the example of the pharmacogene CYP2D6, we demonstrate the potential of WGS-based pharmacogenetic profiling as well as emphasize the limitations of short-read next-generation sequencing. In the near future, we envision a shift toward long-read sequencing as the predominant method for gene diagnostics and pharmacogenetic profiling, providing unprecedented data quality and improving patient care.

High prevalence of CYP2D6 ultrarapid metabolizers in a mestizo Colombian population in relation to Hispanic mestizo populations

Background: Interethnic differences in CYP2D6 allele frequency have been demonstrated across Latin-American countries. Only one previous study describing CYP2D6 genotypes in Colombian population has been performed. Thus, this study aimed to evaluate the CYP2D6 genetic variability in a mestizo Colombian population, as well as the similarities and differences concerning other Hispanic mestizo (HM) populations. Methodology: Two hundred and twelve unrelated healthy Colombian subjects were studied, in which different CYP2D6 polymorphisms were analyzed by extra long-PCR and real-time PCR. Results & discussion: A high percentage of ultrarapid metabolizers (18.4%) was found, representing the highest frequency calculated within the HM populations studied. However, the percentage of poor metabolizers (4.7%) was similar to those previously reported in HM populations.

A Review of the Important Role of CYP2D6 in Pharmacogenomics

Cytochrome P450 2D6 (CYP2D6) is a critical pharmacogene involved in the metabolism of ~20% of commonly used drugs across a broad spectrum of medical disciplines including psychiatry, pain management, oncology and cardiology. Nevertheless, CYP2D6 is highly polymorphic with single-nucleotide polymorphisms, small insertions/deletions and larger structural variants including multiplications, deletions, tandem arrangements, and hybridisations with non-functional CYP2D7 pseudogenes. The frequency of these variants differs across populations, and they significantly influence the drug-metabolising enzymatic function of CYP2D6. Importantly, altered CYP2D6 function has been associated with both adverse drug reactions and reduced drug efficacy, and there is growing recognition of the clinical and economic burdens associated with suboptimal drug utilisation. To date, pharmacogenomic clinical guidelines for at least 48 CYP2D6-substrate drugs have been developed by prominent pharmacogenomics societies, which contain therapeutic recommendations based on CYP2D6-predicted categories of metaboliser phenotype. Novel algorithms to interpret CYP2D6 function from sequencing data that consider structural variants, and machine learning approaches to characterise the functional impact of novel variants, are being developed. However, CYP2D6 genotyping is yet to be implemented broadly into clinical practice, and so further effort and initiatives are required to overcome the implementation challenges and deliver the potential benefits to the bedside.