Design and implementation of a novel pharmacogenetic assay for the identification of the CYP2D6*10 genetic variant
BMC Research Notes volume 15, Article number: 104 (2022)
Tamoxifen is considered to be the most widely used adjuvant therapy for hormone receptor positive breast cancer in premenopausal women. However, it is reported that nearly 30% of patients receiving tamoxifen therapy have shown reduced or no benefits. This may be due to the high inter-individual variations in the CYP2D6 gene that is involved in tamoxifen metabolism. The CYP2D6*10 gene variant (rs1065852C>T) is reported to be commonly found in Asian and South Asian populations. The present study was undertaken to design a novel pharmacogenetic assay (Single step-Tetra Arms Polymerase Chain Reaction) for the identification of the CYP2D6*10 variant and implement the designed assay by genotyping a cohort of breast cancer patients.
The novel assay was successfully designed, optimized and validated using Sanger sequencing. Blood samples from 70 patients were genotyped. The following bands were observed in the gel image: Control band at 454 bp; band for C allele at 195 bp; band for T allele at 300 bp. The genotype frequencies for the CYP2D6*10 (rs1065852C>T) variant were: CC-24.28% (17/70), CT-75.71% (53/70), TT-0% (0/70). The allele frequencies were: T-allele-37.86% and C-allele-62.14%.
Breast cancer is the most common cancer among women and overall, it is the second most common cancer in the world . In Sri Lanka, it is the most common cancer among females with a prevalence rate of 27% . Although tamoxifen has been the drug of choice to treat hormone receptor positive breast cancer in premenopausal women, nearly a one third of patients treated with adjuvant tamoxifen therapy do not obtain the desired benefits. It is reported that there is a substantial inter-individual variability in the response to this drug and the beneficial and adverse effects also appear to be variable and unpredictable for individual patients . Variants in the CYP2D6 gene have been reported to affect the metabolism of tamoxifen in breast cancer patients resulting in variable responses to the drug due to reduced enzyme activity. Tamoxifen is converted to its more active form endoxifen by the CYP2D6 enzyme . However, due to marked variability in the CYP2D6 genotypes among population groups, it has been difficult to formulate a common guideline for using the genotype as a determinant factor in prescribing tamoxifen. Although many pharmacogenomic studies have been conducted to evaluate the possibility of incorporating CYP2D6 genotyping into the therapeutic management of hormone receptor positive breast cancer, conflicting results have posed as a barrier to achieving this. In order to support the establishment of pharmacogenomic testing for tamoxifen therapy in breast cancer patients, it is fundamental to establish prevalence data of the most common variants in under-represented populations. Existing data reports CYP2D6*10:rs1065852C > T as one of the most common variants in the South Asian and Sri Lankan populations that is known to affect tamoxifen metabolism. Knowledge of the prevalence of this variant would be of importance to gain the maximum benefit associated with the use of this drug . The current study was undertaken to design a novel pharmacogenetic assay for easy detection of the CYP2D6*10:rs1065852C>T variant and implement the designed assay by genotyping a cohort of Sri Lankan breast cancer patients.
The study was conducted at the Human Genetics Unit, Faculty of Medicine, University of Colombo. It was an experimental study, where a novel assay was designed for the targeted variant, and a cohort of hormone receptor positive breast cancer patients were genotyped for the CYP2D6*10 variant using the optimized and validated assay. An existing resource of stored venous blood samples obtained from 70 hormone receptor positive breast cancer patients which were collected into EDTA tubes and stored at − 20 °C were used for the present study.
The primers for the CYP2D6*10 rs1065852 variant were newly designed using web-based tools. CYP2D6 gene sequence was obtained from the NCBI (CYP2D6 ID: 1565). The Primer1 is a special software designed for obtaining primer tetrads that are specific for T-ARMS PCR technique . The Primer1 software was operated through the web interface where an input of a target DNA sequence is enabled. The maximum size of the DNA sequence is up to 1000 bases in the 5′ to 3′ direction. Serial cloner enabled cutting a 1000 bp fragment (which included the variant in the middle) that was inserted to Primer1. Other parameters were also considered for primer design: the desired minimum, maximum and average melting temperatures (°C) for a primer oligo; the millimolar concentration of salt and the nanomolar concentration of annealing oligos in the PCR (this is used in the calculation of Tm—melting temperature); the % GC content (the minimum and maximum allowable for any primer); the primer lengths and the minimum and maximum complementarity and product sizes . In par with the selected parameters, the program generated primer sequence sets. To select the optimum primer sets the melting temperatures were considered for each primer set (inner and outer). NCBI primer blast was used to check the specificity of the primers. The details of the designed set of primers, and their product sizes are given in Table 1.
Genomic DNA was extracted using QIAamp® DNA mini kits (Qiagen Ltd., UK) according to the manufacturer’s protocol. The Thermal Cycler (Bio-Rad) was used to conduct gradient PCR for the targeted variant CYP2D6*10. The primer annealing temperatures mentioned in the primer vials by the manufacturer were considered for the annealing temperature gradient in the first attempt at PCR optimization. After a series of optimization attempts, the final optimized cyclic conditions for this variant consisted of an initial denaturation of 3 min at 94 °C (initial denaturation), followed by 35 PCR cycles of 1 min at 94 °C (denaturation), 30 s at 65 °C (annealing temperature), 1 min at 72 °C and the final extension of 5 min at 72 °C. The optimized PCR master mixture contained 5.0 µL of 5 × Green Go-Taq Flexi Buffer, 1.0 µL of MgCl2 (25 mM), 1.0 µL dNTP mixture (2.5 mM each), 1.5 µL of 12.5 mM CYP2D6*10-CF, 1.5 µL of 12.5 mM CYP2D6*10-CR, 1.0 µL of 12.5 mM CYP2D6*10-F (C) and 0.75 µL of 12.5 mM CYP2D6*10-R (T), 1.5 µL of 5 mM Betaine and dH2O. 8 μL of the PCR products of each protocol were separated based on 2% agarose gel at 80 V for 1 h and visualized by ethidium bromide under a UV illuminator.
Four randomly selected samples were validated by Sanger sequencing. The Sanger validated samples were used as controls to compare the genotyping results of the remaining samples [under Additional file information, Additional file 1: Figure S1A (homozygous wild type), Figure S1B (heterozygous variant), Figure S1C (homozygous wild type) and Figure S1D (homozygous wild type) show the Sanger sequencing electropherogram images of the variant].
The bands for each genotype were observed in the gel images: Control band at 454 bp; band for C allele at 195 bp; band for T allele at 300 bp. The genotype frequencies were determined using the data (allele counts) observed from the gel images. They were systematically recorded, and genotype frequencies for homozygous wild-type (CC), heterozygous variant (CT), and homozygous variant (TT) genotypes and the allele frequencies for C (wild type) and T (variant) alleles were calculated using an online calculator (available from: https://wpcalc.com/en/equilibrium-hardy-weinberg/) and tested for Hardy‑Weinberg equilibrium. Standard descriptive statistics were used to analyze the demographic data.
All patients in the sample cohort that was genotyped were females with hormone receptor positive breast cancer. Their mean age was 60.12 ± 6.56 years.
T-ARMS PCR protocol was successfully designed and implemented for the identification of the CYP2D6*10 variant. The expected gel band patterns were observed in the gel image results (Control band: 454 bp; band for C allele: 195 bp; band for T allele: 300 bp). Under Additional file information, Additional file 2: Figure S2 shows the expected gel band patterns for the targeted variant.
Figure 1 shows an actual gel image of the genotype results.
The gel images of all 70 patient samples (with the sample numbers indicated in table) are given in Additional file 3: Figure S3.
The genotype and allele frequencies are shown in Table 2.
Out of the 70 samples which were genotyped, homozygous wild type (CC) and heterozygous (CT) samples were detected but no homozygous samples were detected for the variant T allele (Table 2). The predominant genotype was CT (75.7%) and the minor allele frequency of the CYP2D6*10 variant in this cohort was 37.9%. The results of all the genotyped samples are shown in Additional file 4: Table S1.
The T-ARMS PCR is a technique which many researchers are beginning to adopt as it is a cost-effective method compared to other techniques that are used to detect SNVs [8, 9]. However, the optimization step can be time-consuming and laborious . The presence of GC rich regions near the SNV of interest might be a limiting factor for the use of this technique . One of the common problems faced with T-ARMS PCR is the specificity, where the MgCl2 concentration, primer ratio, and dNTP ratio can play a vital part in the optimization process . Different gradients of MgCl2, dNTP and annealing temperatures were used during the optimization process in this study. The use of Betaine  and Dimethyl sulfoxide (DMSO)  is also commonly reported across studies that have used this technique. In this study, Betaine was used in a lower concentration and volume during the optimization process.
The CYP2D6*10 variant has been studied in different populations and the present study is an addition to this data from the South Asian region. The genotyping results obtained were comparable to data available from a previous Sri Lankan study conducted by Tharanga et al. in 2013, where the CYP2D6*10 minor allele frequency in a Sri Lankan population (n = 75) was reported to be 39%, similar to the 37.9% obtained in the cohort of 70 hormone receptor positive breast cancer patients genotyped in this study. The study by Tharanga et al. further identified that the minor allele frequency was highest among the Moors and 50% of the participants in that study were males . The current study did not enumerate the ethnic specific variations since the ethnicity of the breast cancer patients was not available for subgroup analysis.
A novel T-ARMS PCR assay was successfully designed, optimized and validated to genotype the CYP2D6*10 variant and could be implemented as a cost-effective technique. Homozygous wild type (CC) and heterozygous variant (CT) alleles in CYP2D6*10 were found in the genotyped samples. Homozygotes for the variant allele were not identified in this cohort.
The T-ARMS PCR assay development is challenging specially when the target region has a high GC percentage which could lead to prolongation of the optimization process.
With the use of Betaine, the intensity of the control and other bands may reduce, this could cause difficulty in detecting the relevant bands on the gel image.
Primer dimer formation.
The small sample size used for genotyping in this study. A larger study would be ideal for frequency estimation of the targeted variant in the population.
Availability of data and materials
The following databases were accessed to obtain the details and FASTA sequence related to the gene of interest: NCBI—Gene—CYP2D6—https://www.ncbi.nlm.nih.gov/gene/1565; NCBI RefSeq—CYP2D6—https://www.ncbi.nlm.nih.gov/nuccore/NG_008376.4?report=fasta; SNPedia—CYP2D6 gene—https://www.snpedia.com/index.php/CYP2D6.
- CYP2D6 :
Cytochrome P450 2D6
- T-ARMS PCR:
Tetra Primer Amplification Refractory Mutation System based Polymerase Chain Reaction
Polymerase Chain Reaction
Single Nucleotide Variant
National Center for Biotechnology Information
Basic Local Alignment Search Tool
- F (C):
Forward primer for C allele
- R (T):
Reverse primer for T allele
Ferlay J, Colombet M, Soerjomataram I, Mathers C, Parkin DM, Piñeros M, Znaor A, Bray F. Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int J Cancer. 2019;144(8):1941–53.
Cancer incidence data, 2011; National Cancer Control Programme; Ministry of Health, Nutrition and Indigenous medicine Sri Lanka; Sri Lanka Cancer Registry, 13 publication at : https://www.nccp.health.gov.lk/images/PDF_PUBLICATIONS/Cancer_Incidence_Data_2011.pdf. Accessed 23 Dec 2021.
Ferraldeschi R, Newman WG. The impact of CYP2D6 genotyping on tamoxifen treatment. Pharmaceuticals. 2010;3(4):1122–38.
Saghafi F, Salehifar E, Janbabai G, Zaboli E, Hedayatizadeh-Omran A, Amjadi O, Moradi S. CYP2D6* 3 (A2549del),* 4 (G1846A),* 10 (C100T) and* 17 (C1023T) genetic polymorphisms in Iranian breast cancer patients treated with adjuvant tamoxifen. Biomedical reports. 2018;9(5):446–52.
Dean L. Tamoxifen therapy and CYP2D6 genotype. Medical Genetics Summaries. 2014.
Collins A, Ke X. Primer1: primer design web service for tetra-primer ARMS-PCR. Open Bioinform J. 2012;6(1):55.
Owen RP, Sangkuhl K, Klein TE, Altman RB. Cytochrome P450 2D6. Pharmacogenet Genomics. 2009;19(7):559.
Jin C, Li Z, Zheng X, Shen K, Chao J, Dong Y, Huang Q, Yin Q, Deng Y, Zhu W. Development and validation of T-ARMS-PCR to detect CYP2C19* 17 allele. J Clin Lab Anal. 2020;34(1):e23005.
Ehnert S, Linnemann C, Braun B, Botsch J, Leibiger K, Hemmann P, Nussler AK. One-step ARMS-PCR for the DETECTION of SNPs—Using the example of the PADI4 gene. Methods Protoc. 2019;2(3):63.
Medrano RF, de Oliveira CA. Guidelines for the tetra-primer ARMS–PCR technique development. Mol Biotechnol. 2014;56(7):599–608.
Alyethodi RR, Singh U, Kumar S, Alex R, Deb R, Sengar GS, Raja TV, Prakash B. T-ARMS PCR genotyping of SNP rs445709131 using thermostable strand displacement polymerase. BMC Res Notes. 2018;11(1):1–5.
Tharanga TP, Jinadasa CM, Risama MF, Galappatthy P, Jayakody RL, Dissanayake VH. Genetic variants in the cytochrome P450 2D6 gene in the Sri Lankan population. Indian J Hum Genet. 2013;19(4):392.
The authors extend their sincere appreciation to all staff at the Human Genetics Unit, Faculty of Medicine, University of Colombo for their support.
Ethics approval and consent to participate
An existing resource of stored venous blood samples obtained from hormone receptor positive breast cancer patients who had given informed written consent at the time of sample collection for their samples to be used for future genetic studies at the Human Genetics Unit along with the approval of the Ethics Review Committee, Faculty of Medicine, University of Colombo (EC-15–082) was used for the investigations in the present study.
Consent to publication
The authors declare that there are no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
, S1B, S1C and S1D: Sanger confirmation of CYP2D6*10 homozygous wild type (CC) and heterozygous variant (CT) samples. S1A Homozygous wild type sample. S1B Heterozygous variant type sample. S1C Homozygous wild type sample. S1D Homozygous wild type sample
Gel band patterns of CYP2D6*10 variant
Gel images of all genotyped samples for CYP2D6*10 variant
Genotyping results and age of the hormone receptor positive breast cancer patient cohort
About this article
Cite this article
Ranadeva, N.D.K., Sirisena, N.D., Wetthasinghe, T.K. et al. Design and implementation of a novel pharmacogenetic assay for the identification of the CYP2D6*10 genetic variant. BMC Res Notes 15, 104 (2022). https://doi.org/10.1186/s13104-022-05993-6