Pharmacogenomic assessment of Mexican and Peruvian populations

Association of CYP3A4-392A/G, CYP3A5-6986A/G, and ABCB1-3435C/T Polymorphisms with Tacrolimus Dose, Serum Concentration, and Biochemical Parameters in Mexican Patients with Kidney Transplant

Tacrolimus (TAC) is an immunosuppressant drug that prevents organ rejection after transplantation. This drug is transported from cells via P-glycoprotein (ABCB1) and is a metabolic substrate for cytochrome P450 (CYP) 3A enzymes, particularly CYP3A4 and CYP3A5. Several single-nucleotide polymorphisms (SNPs) have been identified in the genes encoding CYP3A4, CYP3A5, and ABCB1, including CYP3A4-392A/G (rs2740574), CYP3A5 6986A/G (rs776746), and ABCB1 3435C/T (rs1045642). This study aims to evaluate the association among CYP3A4-392A/G, CYP3A5-6986A/G, and ABCB1-3435C/T polymorphisms and TAC, serum concentration, and biochemical parameters that may affect TAC pharmacokinetics in Mexican kidney transplant (KT) patients.

METHODS

Forty-six kidney transplant recipients (KTR) receiving immunosuppressive treatment with TAC in different combinations were included. CYP3A4, CYP3A5, and ABCB1 gene polymorphisms were genotyped using qPCR TaqMan. Serum TAC concentration (as measured) and intervening variables were assessed. Logistic regression analyses were performed at baseline and after one month to assess the extent of the association between the polymorphisms, intervening variables, and TAC concentration.

RESULTS

The GG genotype of CYP3A5-6986 A/G polymorphism is associated with TAC pharmacokinetic variability OR 4.35 (95%CI: 1.13-21.9; p = 0.0458) at one month of evolution; in multivariate logistic regression, CYP3A5-6986GG genotype OR 9.32 (95%CI: 1.54-93.08; p = 0.028) and the use of medications or drugs that increase serum TAC concentration OR 9.52 (95%CI: 1.79-88.23; p = 0.018) were strongly associated with TAC pharmacokinetic variability.

CONCLUSION

The findings of this study of the Mexican population showed that CYP3A5-6986 A/G GG genotype is associated with a four-fold increase in the likelihood of encountering a TAC concentration of more than 15 ng/dL. The co-occurrence of the CYP3A5-6986GG genotype and the use of drugs that increase TAC concentration correlates with a nine-fold increased risk of experiencing a TAC at a level above 15 ng/mL. Therefore, these patients have an increased susceptibility to TAC-associated toxicity.

CYP2D6 and CYP2C19 Genes Associated with Tricontinental and Latin American Ancestry of Pe-ruvians

Precision medicine seeks to individualize the dose from the beginning of phar-macological therapy based on the characteristics of each patient, genes involved in the metabolic phenotype, ethnicity or miscegenation, with the purpose to minimize adverse effects and optimize drug efficacy. The objective was to re-view studies that describe the association of the CYP2D6 and CYP2C19 genes with the tricontinental and Latin American ancestry of Peruvians. A biblio-graphic search was carried out in PubMed/Medline and SciELO, with various descriptors in Spanish and English. The results of this review confirm that the ethnic origin of Peruvians is triconti-nental due to European (mainly Spanish), African and Asian migration, in addi-tion to Latin American migration, being 60.2% mixed, 25.8% Amerindian, 5.9% white, 3.6% African descent, 1.2% Chinese and Japanese descent, and 3.3% unspecified. Studies on CYP2C19*3, CYP2D6*2, *3 and *6 have been reported in Peruvians, and the frequency is similar to that studied in Ecuadori-ans and Colombians. The CYP2C19*3, CYP2D6*3, and CYP2D6*6 alleles found in Peruvians are common in Europeans, Africans, and Asians; while CYP2D6*4 in Africans and CYP2D6*2 related to Asians. In some studies, the ethnic/gene association has not been demonstrated; while others have shown a significant association, which is why further investigation is warranted. It is concluded that the studies on CYP2D6 and CYP2C19 genes associated with the tricontinental and Latin American ancestry of Peruvians are little, and ac-cording to what has been investigated, the CYP2C19*3, CYP2D6*2, *3, *4 and *6 alleles have more related to their ancestry.

C677T and A1298C MTHFR gene polymorphisms and response to fluoropyrimidine-based chemotherapy in Mestizo patients with metastatic colorectal cancer

OBJECTIVE

To assess the association between C677T and A1298C methylenetetrahydrofolate reductase (MTHFR) single-nucleotide polymorphisms (SNPs) and response to first-line fluoropyrimidine-based chemotherapy for metastatic colorectal adenocarcinoma.

METHODS

A total of 68 patients were prospectively followed up in San Juan de Dios Hospital (San José, Costa Rica) from January 2019 to November 2020. Patients received first-line therapy with capecitabine or 5-fluorouracil in combination with oxaliplatin or irinotecan. Germline and somatic DNA was extracted from blood samples and paraffin-embedded tissue, respectively. Overall response rate (partial response + complete response) was assessed according to RECIST 1.1 criteria. Cox regression models were performed to identify the effect of MTHFR C677T and A1298C SNPs on progression-free survival (PFS) and overall survival (OS) (NCT registration number: 03852290).

RESULTS

Patients harboring one or both T alleles of the MTHFR C677T SNP had better overall response than homozygous wild-type individuals [odds ratio (OR): 3.21; 95% confidence interval (CI), 1.05-9.81; P = 0.03]. No association was found between the MTHFR A1298C genotypes and overall response (OR: 0.75; 95% CI, 0.26-2.20; P = 0.60). Patients with the MTHFR 677 TT and CT genotypes had longer PFS than CC individuals (hazard ratio: 0.53; 95% CI, 0.28-0.98; P = 0.045), even after adjustment for confounders (hazard ratio: 0.50; 95% CI, 0.25-0.98; P = 0.04). We found no association between the MTHFR A1298C SNP and PFS (hazard ratio: 1.35; 95% CI, 0.72-2.55; P = 0.34). None of the SNPs was associated with OS.

CONCLUSION

Patients carrying at least one mutant allele of the MTHFR C677T SNP had a better overall response and longer PFS than wild-type homozygous patients.

Comparison of Eight Technologies to Determine Genotype at the UGT1A1 (TA)n Repeat Polymorphism: Potential Clinical Consequences of Genotyping Errors?

To ensure accuracy of UGT1A1 (TA)_n (rs3064744) genotyping for use in pharmacogenomics-based irinotecan dosing, we tested the concordance of several commonly used genotyping technologies. Heuristic genotype groupings and principal component analysis demonstrated concordance for Illumina sequencing, fragment analysis, and fluorescent PCR. However, Illumina sequencing and fragment analysis returned a range of fragment sizes, likely arising due to PCR "slippage". Direct sequencing was accurate, but this method led to ambiguous electrophoregrams, hampering interpretation of heterozygotes. Gel sizing, pyrosequencing, and array-based technologies were less concordant. Pharmacoscan genotyping was concordant, but it does not ascertain (TA)₈ genotypes that are common in African populations. Method-based genotyping differences were also observed in the publication record (p < 0.0046), although fragment analysis and direct sequencing were concordant (p = 0.11). Genotyping errors can have significant consequences in a clinical setting. At the present time, we recommend that all genotyping for this allele be conducted with fluorescent PCR (fPCR).

Variation in Actionable Pharmacogenetic Markers in Natives and Mestizos From Mexico

The identification and characterization of pharmacogenetic variants in Latin American populations is still an ongoing endeavor. Here, we investigated SNVs on genes listed by the Pharmacogenomics Knowledge Base in 1284 Mestizos and 94 Natives from Mexico. Five institutional cohorts with NGS data were retrieved from different research projects at INMEGEN, sequencing files were filtered for 55 pharmacogenes present in all cohorts to identify novel and known variation. Bioinformatic tools VEP, PROVEAN, and FATHMM were used to assess, in silico, the functional impact of this variation. Next, we focused on 17 genes with actionable variants that have been clinically implemented. Allele frequencies were compared with major continental groups and differences discussed in the scope of a pharmacogenomic impact. We observed a wide genetic variability for known and novel SNVs, the largest variation was on UGT1A > ACE > COMT > ABCB1 and the lowest on APOE and NAT2. Although with allele frequencies around 1%, novel variation was observed in 16 of 17 PGKB genes. In Natives we identified 59 variants and 58 in Mestizos. Several genes did not show novel variation, on CYP2B6, CYP2D6, and CYP3A4 in Natives; and APOE, UGT1A, and VKORC1 in Mestizos. Similarities in allele frequency, comparing major continental groups for VIP pharmacogenes, hint towards a comparable PGx for drugs metabolized by UGT1A1, DPYD, ABCB1, CBR3, COMT, and TPMT; in contrast to variants on CYP3A5 and CYP2B6 for which significant MAF differences were identified. Our observations offer some discernment into the extent of pharmacogenetic variation registered up-to-date in Mexicans and contribute to quantitatively dissect actionable pharmacogenetic variants in Natives and Mestizos.

Geographic distribution of the 3435C>T polymorphism of the MDR1 gene in Peruvian populations

Background The MDR1 gene presents several genetic polymorphisms with pharmacological implications. Therefore, the aim of the present study is to establish the genotype and allele frequencies of 3435C>T polymorphism of MDR1 gene into Peruvian populations (Coastal, Andean and Amazonian ecoregions), even considering the altitude (lowland <2500 m and highland >2500 m). Methods The polymorphism was analyzed by TaqMan genotyping assays in a group of 181 healthy unrelated Peruvian individuals. The comparison of genotype and allele frequencies of 3435C>T polymorphism was made with the Pearson test (X2), and, to calculate the genotype distributions, the Hardy-Weinberg equilibrium (HWE) was used. Results In all populations evaluated in this study, the genotype frequency distributions met HWE assumptions. The comparison between genotype and allele frequencies showed significant differences (p < 0.05), when the Andean, Coastal and Amazonian populations were compared. Also, significant differences (p < 0.05) were obtained when these populations were compared considering their altitudes. Likewise, in comparison with countries like USA, Finland, Nigeria and Kenya, the results showed significant differences (p < 0.05). Conclusions This investigation allowed us to establish the genotype and allele frequencies of 3435C>T polymorphism in different Peruvian populations, considering the geographic localization and even the altitude.

Population-scale genomics-Enabling precision public health

The current excitement for affordable genomics technologies and national precision medicine initiatives marks a turning point in worldwide healthcare practices. The last decade of global population sequencing efforts has defined the enormous extent of genetic variation in the human population resulting in insights into differential disease burden and response to therapy within and between populations. Population-scale pharmacogenomics helps to provide insights into the choice of optimal therapies and an opportunity to estimate, predict and minimize adverse events. Such an approach can potentially empower countries to formulate national selection and dosing policies for therapeutic agents thereby promoting public health with precision. We review the breadth and depth of worldwide population-scale sequencing efforts and its implications for the implementation of clinical pharmacogenetics toward making precision medicine a reality.

High frequency of mutant thiopurine S-methyltransferase genotypes in Mexican patients with systemic lupus erythematosus and rheumatoid arthritis

Systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA) are treated with immunosuppressive purine analogs, 6-mercaptopurine/6-thioguanine/azathiopurine, which are inactivated by thiopurine S-methyltransferase (TPMT). Non-synonymous polymorphisms in TPMT are associated with increased risk of adverse effects in patients treated with thiopurines. This study aimed to determine the frequency of the most common mutant TPMT alleles in Mexican patients with SLE (a prototype autoimmune disease) and RA (one of the most common autoimmune diseases in Mexico). Five hundred fifty-three consecutive patients from Central Mexico with SLE (178) and RA (375) were included. Subjects were genotyped to identify TPMT*2 (rs1800462), TPMT*3A (rs1800460 and rs1142345), TPMT*3B (rs1800460), and TPMT*3C (rs1142345) mutant alleles. DNA samples were assayed with the 5' exonuclease technique and TaqMan probes. Mutant alleles were detected in 6.2 and 5.2% of SLE and RA cases, respectively. Of note, 12.4% of SLE cases and 10.1% of RA cases carried mutant genotypes. Among those, the null genotype (TPMT*2/*3A, 0.3%) and the TPMT*3B (0.5%) and TPMT*3C (1.0%) alleles were found in RA, but not SLE cases. Mexican SLE cases displayed the highest frequency of mutant TPMT genotypes worldwide. TPMT genotyping should be performed for Mexican patients with SLE and RA before prescribing purine analogs.

Interethnic variability of pharmacogenetic biomarkers in Mexican healthy volunteers: a report from the RIBEF (Ibero-American Network of Pharmacogenetics and Pharmacogenomics)

Mexico presents a complex population diversity integrated by Mexican indigenous (MI) (7% of Mexico's population) and Mexican mestizos (MMs). This composition highlights the importance of pharmacogenetic studies in Mexican populations. The aims of this study were to analyze the reported frequencies of the most relevant pharmacogenetic biomarkers and metabolic phenotypes in healthy volunteers from Mexican populations and to assess its interethnic variability across MI and MM populations. After a literature search in PubMed, and according to previously defined inclusion criteria, 63 pharmacogenetic studies performed in Mexican healthy volunteers up to date were selected. These reports comprised 56,292 healthy volunteers (71.58% MM). Allele frequencies in 31 pharmacogenetic biomarkers, from 121 searched, are described. Nine of these biomarkers presented variation within MM and MI groups. The frequencies of CYP2D6*3, *4, *5, *10, *17, *35 and *41 alleles in the MM group were different from those reported in the MI group. CYP2C9*2 and *3 alleles were more frequent in MM than in MI populations (χ2 test, p<0.05). CYP2C19*3 allele was not found in the MM or MI populations reported. For UGT1A1*28, only one study was found. HLA-A*31:01 and HLA-B*15:02 were present in some MM and MI populations. Poor metabolizers for CYP2D6 and CYP2C9 were more frequent in MM than in MI groups (χ2 test, p<0.05). Only 26% of the relevant pharmacogenetic biomarkers searched have been studied in Mexican healthy volunteers; therefore, further studies are warranted. The frequency variation of biomarkers in MM and MI populations could be important for the clinical implementation of pharmacogenetics in Mexico.

IL1A rs1800587 associates with chronic noncrisis pain in sickle cell disease

AIM

Pain is prevalent in sickle cell disease (SCD) patients who display great heterogeneity in pain severity and frequency. Hypothesizing that inflammatory factors are involved in the pathogenesis of SCD pain, we focused on the IL1A C/T polymorphism rs1800587 that is an SNP located in a cis-transcriptional regulatory region.

METHODS

We genotyped IL1A rs1800587 and performed association studies with phenotype data obtained by a multidimensional pain assessment tool using the PAINReportIt^® Questionnaire.

RESULTS

Each T allele was associated with a 3.9 increase in composite pain index score (p = 0.04) as determined by multiple linear regression.

CONCLUSION

IL1A rs1800587 may influence chronic pain in SCD.

Genetic variability of CYP2C9*2 and CYP2C9*3 in seven indigenous groups from Mexico

Aim: CYP2C9 is one of the major drug metabolizing enzymes, however, little is known about polymorphisms in CYP2C9 gene and pharmacological implications in Mexican indigenous populations. Thus, frequencies of CYP2C9*2 and CYP2C9*3 alleles were evaluated in indigenous groups located in northwest (Cora), center (Mazahua and Teenek), south (Chatino and Mixteco) and southeast (Chontal and Maya) regions Mexico. Materials & methods: Allelic discrimination was performed by real-time PCR. Results: CYP2C9*2 allele was found only in Chontal and Maya groups, despite the low contribution of Caucasian component in these populations. CYP2C9*3 allele was present in all populations except in Mazahua, showing a wide genetic variability in the studied populations. Interestingly, we found significant differences between indigenous groups in CYP2C9*3 allele, even in groups located at the same region and belonging to the same linguistic family. Conclusion: These results contribute to laying the pharmacogenetic bases in Mexico, in addition to improving treatment, taking into account the genetic interethnic differences.

Prevalence of pain-related single nucleotide polymorphisms in patients of African origin with sickle cell disease

BACKGROUND

Prospective pain genetics research is hindered by a lack of data on the prevalence of polymorphisms in pain-relevant genes for patients with sickle cell disease (SCD). For African-Americans in general, limited information is available in public databases.

METHODS

We prioritized and examined the genotype and allele frequencies of 115 SNPs from 49 candidate pain genes in 199 adult African-Americans and pediatric patients of African origin with SCD. Analyses were performed and compared with available data from public databases.

RESULTS

Genotype and allele frequencies of a number of SNPs were found to be different between our cohort and those from the databases and between adult and pediatric subjects.

CONCLUSION

As pain therapy is inadequate in a significant percentage of patients with SCD, candidate pain genetic studies may aid in designing precision pain medicine. We provide prevalence data as a reference for prospective genetic studies in this population.