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.

Introduction

Pharmacogenetic studies in Mexican populations have been directed towards the identification of markers previously reported with functional consequences for drug safety and efficacy (Contreras et al., 2011; Bonifaz-Pena et al., 2014; Marsh et al., 2015; Fricke-Galindo et al., 2016). The advent of faster and more accessible technologies has made feasible to investigate a broader swath of the pharmacogenome including the presence of novel variants. In most of the developed world pharmacogenetic testing is being implemented with the aid of consortia and institutions such as the PGRN, the PGKB, CPIC, and the FDA. These offer a list of genetic variants to guide drug selection, dose optimization, and to reduce the risk of adverse drug reactions. Moreover, implementation programs such as The electronic Medical Records and Genomics (eMERGE) initiative have started preemptive pharmacogenomics testing in over 10,000 patients, on the basis of the benefits of pharmacogenetic information, despite the so-called lack of cost–benefit evidence (van Driest et al., 2014). Nevertheless, in developing countries broad pharmacogenetic implementation is even more incipient, a curated collection of local variants has not been yet defined, and for many relevant markers there are significant differences for minor allele frequencies (MAF) when comparing to major continental groups. For example, Campos et al. reported lower CYP2C19*17 rs1224856 allele frequencies in 346 Mexican Americans compared to CEU, MAF: 0.14 vs. 0.22, which may be indicative of lower frequency of bleeding with clopidogrel in the former (Claudio-Campos et al., 2015). In Mexicans higher allele frequencies have been reported for CYP3A4 rs2750574 and NQO1 rs1800566, which may be differentially affecting the efficacy of tacrolimus, fluorouracil, and anthracyclines, when compared to Europeans.

Mexico is home to 68 genetically different ethnic groups which suggests urgency for an adequate collection, classification, and characterization of variants on genes that affect drug efficacy and safety. Published studies collect the identification and determination of allele frequencies of at least all PGKB level 1 variants in Mestizos (Cuautle-Rodriguez et al., 2014; Fricke-Galindo et al., 2014; Fricke-Galindo et al., 2016; Gonzalez-Covarrubias et al., 2017). Fewer studies have included Natives, NGs techniques, or genotype-phenotype assessments. Results accede that for several makers, allele frequencies and its clinical impact differ among Natives, Mestizos and, major continental groups (Bonifaz-Pena et al., 2014; Fricke-Galindo et al., 2016; de Andrés et al., 2017) In addition, little is known of the presence and frequency of local private variation, which can only be investigated though sequencing (Moreno-Estrada et al., 2014). NGS studies have also shown that it is the interplay of several variants, rather than one maker, that explains drug response variability (Johnson et al., 2011; Beekman et al., 2013; Katsila and Patrinos, 2015; Yang et al., 2016). In this regard, Han S.M. et al. identified dozens of rare variants with functional consequences to drug response, but more importantly, it revealed that targeted sequencing enabled profiling of actionable and rare variants, some of yet unknown functional consequences (Han et al., 2017). Furthermore, in vitro studies of major pharmacogenes have confirmed that many rare variants do have a functional impact confirming NGS bioinformatic predictions (Matimba et al., 2009; Han et al., 2017). More limited in scope, our own work has documented the importance of variants, derived from NGS, to increase the predictive value of genotypes on the pharmacokinetics of atorvastatin (Cruz-Correa et al., 2017) and the pharmacodynamics of cumarins (Gonzalez-Covarrubias et al., 2017). Nevertheless, the collection of a comprehensive pharmacogenetic set of variation is still underway, mostly for understudied populations.

Here, we identified, and classified pharmacogene variation on 17 genes listed by the PGKB in 1284 Mestizos and 94 Natives from several institutional cohorts at the National Institute for Genomic Medicine in Mexico, from which NGS data were available. Allele frequency comparisons showed that pharmacogenetic variation is significantly different between Natives and Mestizos for several PGKB variants. We observed twice as many actionable variants in Mestizos vs. Natives, but 4× more novel variation per individual in Natives suggesting that the pharmacogenome in the latter is far from complete, and that current pharmacogenetic guidelines may prove to be more beneficial for certain populations within the country.

Methods

Participants and Sequencing Files. We analyzed WES data from 1,284 Mexican Mestizos published at the ExAC platform (Lek et al., 2016; Karczewski et al., 2017) and three institutional cohorts sequenced with custom probes using different Haloplex v2, Agilent Technologies protocols, thus variant calling was performed independently for each cohort. In addition, sequencing files for 94 Natives were procured by the 100 Genome Consortium (INMEGEN, to be published) and the project Metabolic Analysis in an Indigenous Sample (MAIS) (Contreras-Cubas et al., 2016). All participants provided a blood sample after signing an informed consent and all research projects were approved by the Ethics and Research Committees at INMEGEN. Participants were self-reported as Natives and confirmed by AIMS analyses (to be published). DNA from Natives was sequenced by BGISEQ-500 (Cambridge, MA, USA). Samples from Natives belonged to 36 ethnic groups thus, further stratification and multiple comparisons within Natives or with Mestizos was not possible. Similarly, for Mestizos >90% of all DNA samples came from individuals from the center of the country limiting further stratification. In summary, vcf data came from five institutional cohorts: WGS (Natives = 94, to be published), WES (Mestizos = 968) (Flannick et al., 2019), NGS-targeted 1 (Mestizos = 110) (Gonzalez-Covarrubias et al., 2017), NGS-targeted 12 (Mestizos = 146) (Gonzalez-Covarrubias et al., 2016), and NGS-targeted 3 (Mestizos = 60) (Cruz-Correa et al., 2017). All data were processed according to the Broad Institute recommended best practices workflow and the Genome Analysis ToolKit (GATK) (Auwera et al., 2013).

To assess the functional impact of variants in coding and non-coding regions, we selected three algorithms. We utilized PROVEAN, VEP, and FATHMM which are well published and have been used by the scientific community (Devarajan et al., 2019). PROVEAN predicts whether a SNV affects the function of a protein by giving a score based on sequence alignment; the lower the score, the more dissimilar is the substitution compared to the reference (Mizzi et al., 2014; Wendt et al., 2018). VEP (Variant Effect Predictor) seeks for genes/transcripts of a variant to determine its effect at the amino acid level (e.g. stop gained, missense, stop lost, altered splicing, frameshift, stop loss, and start loss), a high impact predicted by VEP is considered as deleterious (Oetting et al., 2016; Lakiotaki et al., 2017). FATHMM was utilized to predict functional consequences of non-coding variants, it integrates functional annotations from ENCODE and calculates a score (0–1). Scores >0.5 are indicative of deleterious variation (McInnes et al., 2019; Zhou et al., 2018), and although non-coding variants may not have a direct impact on a protein the overall regulation may render faulty. These tools have been utilized in PGx research (Devarajan et al., 2019), and we have reported their use to filter variants with potential effects on atorvastatin pharmacokinetics with acceptable results. PGx-oriented algorithms for the prediction of the functional impact of novel variants are being developed and validated, these new tools include some of the algorithms utilized here (Zhou et al., 2018; Zhou et al., 2019).

Genotypes of PGKB variants in Table 2 were in complete concordance with previous genotyping experiments by RTPCR ex. CYP2C9/VKORC1 (Villegas-Torres et al., 2015), DMET microarray (Gonzalez-Covarrubias et al., 2016) and by NGS (Gonzalez-Covarrubias et al., 2017).

Pharmacogenetic Analyses. Statistical analyses focused on 55 genes since these were identified in all samples. In-depth analyses were undertaken for variants listed by the PharmGKB with a level of evidence 1 and 2 (Relling and Klein, 2011). PGKB level 1 (1A/1B) represent variation with a strong evidence of PK/PD alteration, and these have been published as markers in CPIC dosing guidelines. Level 2 refers to variants with moderate evidence to be linked to a drug’s safety or efficacy. We computed allele frequencies and performed comparisons for several populations, CEU, Northern Europeans from Utah (Caucasians); MXL, Mexicans from Los Angeles; YRI, Yoruba in Ibadan, Nigeria; CHB, Chinese Han from Beijing. Data analyses, descriptive statistics, allele frequency calculations, and variant inferences were performed with R (Team R Core, 2014), PLINK (Purcell et al., 2007), and PGKB resources (Sangkuhl et al., 2008).

Results

NGS Summary. Sequencing files from 94 Natives and 1,284 Mestizos showed a depth-coverage between 50× and 600×, tNGS aimed for 100× coverage while WES and WGS for at least 50×. Depth coverage varied per genomic region regardless of the theoretical depth calculated per sequencing run. For example, tNGS showed a range of depth and coverage of 40×–600× and 80–100% for the selected genes, but after quality control some pharmacogenetic relevant regions including that of VKORC1 rs9323231 and CYP2C19*17 rs12248560 were not covered.

In Natives, preliminary analyses identified the largest number of variants as data came from WGS, NGS-targeted and WES in Mestizos yielded a lower count. Together, all cohorts gathered 436 pharmacogenes, on which we identified on average 102 novel variants per individual in Natives, a value almost three times higher than in Mestizos. Consistently, non-synonymous SNVs (single nucleotide variants) were five times more abundant in Natives, suggesting a genetic diversity uncovered by NGS, but these 436 genes were not identified in all individuals. Therefore, we directed further analyses on genes identified in all cohorts and genes listed by the PGKB with a validation level of 1 or 2, which are the focus of this report.

Pharmacogenetic Analyses

Fifty-five genes were identified in all in 1,378 individuals including most CYPs, FMOS, members of the UGT1A family, SULT1A, and 11 genes involved in pharmacodynamics. Table 1 lists variation for these 55 genes many of which are listed by the eMERGE initiative and the Pharmacogenomics Research Network (Bush et al., 2016). On these, we observed 3.4× and 2.5× more variants in Mestizos vs. Natives (novel and known), as expected given the larger sample of Mestizos. However, the average number of variants per individual was 3.9× and 3.4× higher (novel and known) in Natives, reflecting the larger diversity unreported in the latter. Also, private variants, i.e., polymorphisms observed as heterozygous in only one individual, were more frequent in Natives (n = 868) vs. Mestizos (n = 3256, Table 1).

TABLE 1

Table 1 Variation in 55 pharmacogenes shared among cohorts.¹

Next, we focused on 17 actionable genes with clinically validated variants to perform comparisons between populations. Genes were selected according to the PGKB classification considering only those with a validation level 1 and 2, included in published CPIC dosing guidelines also, considered actionable by the FDA. These 17 genes were, ABCB1, ACE, APOE, CBR3, CYP2C19, CYP2C9, CYP2B6, CYP2D6, CYP3A5, CYP3A4, COMT, DPYD, NAT2, SLCO1B1, TPMT, UGT1A, and VKORC1. On these, we observed 2,387 non-private, novel, and known variants, 34 of which were listed by the PGKB. The largest number of known variants were observed on ACE (251), UGT1A (318), ABCB1 (182), and COMT (235), (Figure 1). Genes with the lowest number of variants were NAT2 (33) and APOE (24), again with a higher count in Mestizos.

FIGURE 1

Figure 1 Novel variation and its functional impact identified in 16 PGKB genes in Natives (gray) and Mestizos (blank).

Next, we determined minor allele frequencies for these 17 PGKB-listed variants and compared them between Mexicans and those reported for major continental populations. In Natives, we observed a higher MAF, compared to Mestizos, for VKORC1 rs8050894 (1.4×), CYP2B6 rs2279343 (1.5×), and CYP3A5 rs776746 (2×), with potential implications for the differential management in Natives for cumarins, efavirenz, nevirapine, methadone, and tacrolimus am ong several others. For most of these variants allele frequencies in Natives were higher compared to Mestizos, also when compared to CEU and CHB (Table 2).

TABLE 2

Table 2 Allele frequencies for 17 actionable pharmacogenetic variants.

The opposite was also observed, i.e., lower allele frequencies in Natives were identified for CYP2C19 rs4244285 (2×), CYP2C9 rs1799853 (4.5×), rs1057910 (2×), NAT2 rs179930 (2×), SLCO1B1 rs4149015, and APOE rs7412. The latter two were not found in Natives. CYP2C19 and CYP2C9 lower allele frequencies in Natives may be indicative of a lower proportion of poor metabolizers potentially affecting the pharmacokinetics of over 18% of all drugs. For the remaining 26 variants on 8 genes, allele frequencies were comparable between Natives and Mestizos. A summary of variation on these 17 genes is presented in Table 2.

Novel Variants on Actionable PGKB Genes

Using NGS data we sought for novel variation, and identified 116 not-previously reported variants on 16 of the 17 PGKB genes (MAF ≥1%). We found some novel variation in all genes studied in Mestizos or Natives except for SLCO1B1. The functional impact of these variants was assessed by the algorithms, VEP, Provean, and FATHMM (Dong et al., 2015). Table 3 lists novel variants that were predicted as deleterious by more than two algorithms, observed on ABCB1 (4), CYP2C19 (1), and CYP3A4 (1), and although their frequency does not exceed 2% we believe their strongly predicted functional impact pinpoints them as worth investigating. A full list of novel variants annotated by the three algorithms is listed in Supplementary Table 1.

TABLE 3

Table 3 Selected novel variants predicted as deleterious.¹

Lack of novel variation was detected for APOE, UGT1A1, and VKORC1 in Mestizos, and for CYP2D6, CYP3A4, CYP2B6, and COMT in Natives (Figure 1). The largest count of novel variants in Natives was observed for ABCB1 (25) and DPYD (15), and in Mestizos on CYP2C19 (10) and CYP3A4/5 (13). Novel variation was on average comparable between population groups, but differed largely per gene Figure 1 and Supplementary Table 1.

Discussion

As genetic variation is better defined within populations it is relevant for each country to assess whether current pharmacogenetic platforms can offer direct benefits for their people. There are 68 Native groups in Mexico representing around 10% of the population, but admixture proportions vary greatly throughout the country. Thus, we sought to determine and compare differences in gene variation related to pharmacokinetics and pharmacodynamics in Natives and Mestizos by consolidating different institutional cohorts. Sequencing data for all 1,378 DNA samples were obtained, by various groups utilizing different NGS approaches, and our analyses focused on 17 PGKB genes. Despite the sample size difference between Natives (94) and Mestizos (1284) we were able to capture representative pharmacogenetic variation for both groups. Not unexpectedly, Natives showed almost 10× more variants, which hints towards a larger proportion of unaccounted variation (novel and known).

Comparative analyses among Mestizos, Natives, and continental populations, highlight the extent of our incomplete registry of pharmacogenetic variability within the country, but also indicate a closer completion for Mestizos than for Natives, because these have been much less studied (Cid-Soto et al., 2018; Sánchez-Pozos et al., 2018). For instance, Natives showed 2.4 functional novel variants on average per individual compared to 0.7 in Mestizos, which is in agreement with Romero-Hidalgo et al. (2017) showing that damaging variation is 2× higher for unreported variants in Native Americans. These observations indicate that genetics may underlie drug disposition differences not only between Mexico and major continental groups, but also within the country’s populations.

PGx Markers in Guidelines for Clinical Implementation

CPIC dosing guidelines follow data curation by the PGKB and represent the most useful pharmacogenetic resource around the world. Therefore, local adoption of pharmacogenetics in developing countries requires identifying the frequency of already published actionable markers, and to assess if their presence warrants its application to local populations prior to implementation. Here, we compared allele frequencies for 34 PGKB variants in 17 genes with a validation level 1 and 2 to assess if they were similar or different within the studied populations, and were also compared to the major ancestral population (Table 2). First, we discuss variants with higher allele frequency in Natives vs. Mestizos vs. CEU. These were observed on VKORC1, CYP3A5, and UGT1A1. Cumarin-sensitivity variants VKORC1 rs8050894, rs994438, and 9923231, all in LD, were 25–40% more frequent in Natives hinting for higher cumarin sensitivity. However, only 1.6% of Natives showed variants on CYP2C9 vs. 5–10% of Mestizo and CEU. In Natives VKORC1 variation may reflect higher sensitivity to cumarins, but the low CYP2C9*2/*3 frequencies a higher metabolism this particular VKORC1–CYP2C9 interplay may abrogate cumarin dosing adjustment in Natives. It is relevant to consider that this is mostly the case for warfarin since other CYPs such as CYP2C8 are responsible for R-acenocumarol metabolism, the isomer to which the anticoagulant effect is ascribed.

CYP3A5*3, rs776746 affects the disposition of over 20 drugs, but it is particularly relevant for cyclosporine and tacrolimus. In Europeans the T allele has a frequency of 4%, significantly lower than the one observed in Natives (MAF: 0.303) and Mestizos (MAF: 0.189). It is possible that these differences may lead to more dose adjustments for tacrolimus in Mexicans compared to CEU (Zhu et al., 2011), but similar to that in Asians (MAF: 0.311) (Niioka et al., 2012).

Variants on UGT1A1 significantly influence the pharmacokinetics of several drugs including irinotecan and atazanavir, lower enzyme levels are a consequence of several polymorphisms including, rs887829 and rs4148323, the former showed similar allele frequencies in Natives, Mestizos, MXL, and CEU, but rs4148323 was 4× and 3× more frequent in Natives and Mestizos compared to CEU, this may be indicative of low glucuronidation in a higher proportion of Mexicans. The overall impact of these differences cannot be conclusive, but it is relevant for many other drugs, for which its clearance rate is limited by UGT1A1. Our NGS data did not identify other key variants such as UGT1A1*28 rs3064774 or rs4148323 which have been shown to be differently distributed in other Latin American countries (Marsh et al., 2015).

We identified 4 of the 24 CYP2B6 variants listed by the PGKB, only rs4803419 (Level 2B) showed differential allele frequencies Natives > Mestizos > MXL > CEU, suggesting that Natives and Mestizos are 2× more likely to have decreased CYP2B6 metabolism due to rs4803419 which is in LD with rs3745274 a Level 1 variant. This is mostly relevant for efavirenz elimination, since CYP2B6 is its major metabolizing enzyme. A recently published guideline considers rs3745274 as the pharmacogenetic marker for efavirenz dosing assessment (Desta et al., 2019).

We identified five variants with lower allele frequency in Natives vs. Mestizos, CEU or MXL on, CYP2C19*2 rs4244285, NAT2 rs1799930, CYP2D6 rs28371725, SLCO1B1 rs41419015, and ABCB1 rs2032582. In addition to the already discussed VKORC1–CYP2C9*2/*3 interplay, our results on CYP2C9 validate previous reports showing that Mestizos are 3× and 4× more likely to show impaired CYP2C9 metabolism than Natives due to CYP2C9*2 and *3 (Villegas-Torres et al., 2015). CYP2C9 is the second most expressed CYP in liver after CYP3A4, and impacts 15% of all drugs including losartan, NSAIDs, phenytoin, and hypoglycemic agents, (Sconce et al., 2005;Van Booven et al., 2010). Here, we hypothesize that Natives will show a lower probability of impaired CYP2C9 activity compared to Mestizos and both of them lower to CEU. These comparisons have been previously confirmed for CYP2C9, ABCB1, CYP3A5, and CYP2C19 (Vargas-Alarcón et al., 2014; Zhou et al., 2017). For CYP2C19, we can infer that a higher proportion of Natives are normal/extensive metabolizers. Clopidrogrel dosing guideline lists around 40 variants on CYP2C19, we identified 3 of the 10 variants reported for the Americas, CYP2C19*3 and *4 were rare (<1%) in Mestizos and CEU, and we did not observe them in Natives. The CPIC reported these variants associated to ADRs and are relevant for dose estimation for clopidogrel, escitalopram, and voriconazole. It is possible that ADR or dose adjustments would be seen in a lower proportion in Mexicans. Similarly, for SLCO1B1 rs4149015 and ABCB1 rs2032582 frequencies were lower in Natives, from we could infer a higher efficacy of drugs such as statins, atazanavir, or sunitinib.

Similarities among Natives, Mestizos, MXL, and CEU were observed for several variants on UGT1A1, DPYD, ABCB1, CBR3, CYP2B6, COMT, and TPMT hinting towards a comparable pharmacokinetics between populations for the role variants observed. Maybe these variants could smoothly transition into clinical implementation, but these inferences should not only be validated, but should also account for unreported, rare, and novel variation to determine the overall pharmacogenetic impact for each gene-drug combination.

Novel Variants on PGKB Genes

Pharmacogenetic research in populations from Mexico has been actively increasing, however the number of publications by 2018 barely reached 0.3% of the >22,000 NCBI pharmacogenetic/pharmacogenomics reports. Here, we identified novel variants on 16 of the 17 studied PGKB genes, in Mestizos (58 variants) and Natives (67 variants, Figures 1 and 2, Table 3 and Supplemental Table 2). For some genes novel variants were identified only in Natives, APOE (1), UGT1A1 (3), and VKORC1 (3) or only in Mestizos, COMT (4), CYP2B6 (1), CYP2D6 (3), and CYP3A4 (6). Although these counts are influenced by sample size, it may offer an estimate of the completion in gene variation. Interesting differences in the number of novel variants arose for ABCB1, ACE, DPYD, UGT1A1, and VKORC1. For example, the transporter gene ABCB1 showed 25 novel variants in Natives, and only 5 in Mestizos, the functional impact of most of these were confirmed by independent algorithms and 10 were validated as deleterious. Similar results were observed for variants on DPYD suggesting that variation on these genes is far from complete in Natives. In the previous section, we mentioned that the allele frequency of several PGKB variants in Natives were similar to that in other populations however, novel variation indicates that it is likely that the overall pharmacogenetic impact has not been fully described. This is in agreement with a few reports on DPYD and ABCB1 novel variation in underrepresented populations, indicating that our current information of pharmacogenetic predictors remains to be thoroughly depicted (Mukonzo et al., 2009; Del Re et al., 2015; Elraiyah et al., 2017).

FIGURE 2

Figure 2 Functional impact of novel variation in Mestizos (gray) and Natives (blank).

Variation on CYP2C19 did not show novel variants in Natives, we may infer that most common variation has been recorded for this gene, and that further phenotypic variability may be imputed to rare mutations. In Mestizos, we found 10 novel variants, 2 were missense deleterious (POS.10:96522531 and POS 10:96609809) which may be relevant for dozens of drugs including, proton pump inhibitors, antiepileptics, antiplatelets, and antidepressants. We anticipated fewer novel variation on this gene since intense research has been done given its high importance in multiple drug-drug interactions and drug metabolism. Nevertheless, 10 novel variants were identified in Mestizos, which may accord with a couple of recent NGS studies, reporting novel, common, and rare variation on this gene in Asians and Africans (Matimba et al., 2009; Han et al., 2017; Dai et al., 2015).

Interestingly, we detected 14 novel variants with a MAF>3% on (4 SNVs MAF:8%), DPYD (1 SNV, MAF:3%), CYP2C9 (1 SNV, 3%), ACE (1 SNV, MAF:5%), APOE (1 SNV, MAF:5%), CYP3A4/5 (5 SNVs MAF 3%), and CYP2D6 (1 SNV, MAF:40%), although none of these were predicted as deleterious (Supplemental Table 1).

We are aware of the different sequencing techniques and size of the study, potentially affecting the number of variants and the inferences made. Also, PGKB variants not reported here reflect either lack of sequencing coverage or an allele frequency lower than 1%. Nevertheless, several of our observations have been confirmed by previous reports (Villegas-Torres et al., 2015) or paralleled those from Latin American countries such as Brazil with which we share pharmacogenetic similarities and differences ex. MAF for CYP2B6 rs3745274, CYP3A5 rs776746, VKORC1 rs8050894, and ABCB1 rs2032582, supporting the notion that pharmacogenetic diversity across the Americas ought to be consolidated.

Conclusions

Our observations summarize variation in 55 pharmacogenes in 1,378 Natives and Mestizos from Mexico, focusing on 17 PGKB genes. This is one of the largest collections of genetic variability related to pharmacogenomics in Mexicans. Our report offers a collection of variants in core pharmacogenes, confirming previous knowledge and contributing to the list of novel variants that can be further investigated and may become a part of a preliminary catalogue for PK/PD, and phenotype-genotype correlations. These results may also complement genotyping platforms with relevant pharmacogenetic variants with specific population background. Future studies will seek to validate this variation and to confirm its potential application in pharmacogenomics.

Data Availability Statement

The datasets analyzed for this study can be found at the databases: dbGAP: WES (Mestizos=968) (Flannick et al. 2019) [https://www.ncbi.nlm.nih.gov/gap/ phs001393 and phs001099], and at database, European Variation Archive (EVA) under project PRJEB343334 for: NGS-targeted 1 (Mestizos =110) (Gonzalez-Covarrubias et al. 2017) [EVA INMEGEN-ACOAG analysis: ERZ1079024], NGS-targeted 12 (Mestizos =146) (Gonzalez-Covarrubias et al. 2016) [INMEGEN-SCZ analysis: ERZ1079026], and NGS-targeted 3 128 (Mestizos = 60) (Cruz-Correa et al. 2017) [INMEGEN-ATV analysis: ERZ1079025]. WGS (Natives=94) Access to this dataset can be requested via email to Dr. Cristobal Fresno of 100G-Consortium, INMEGEN [Y2ZyZXNub0Bpbm1lZ2VuLmdvYi5teA==].

Ethics Statement

The studies involving human participants were reviewed and approved by Comite de Etica e Investigacion, Instituto Nacional de Medicina Genomica #25/2016/I. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

VG-C conceived the analyses and wrote the paper. MM-F performed bioinformatics and statistical analyses. OC-C contributed to data collection. AM-H, HG-O, and FB-O contributed to patient recruitment and genetic experiments. AG-M contributed to patient recruitment and genetic experiments. JM-M recruited patients and performed genetic experiments. HN designed project, contributed to patient recruitment. LO conceived, coordinated, and executed the project. XS conceived and coordinated the project.

Funding

Funding from Conacyt grant no. 252952 to XS, Conacyt grant no. 233970 to LO, Conacyt grant no. 272795 to VG-C and grant INMEGEN-No.14/2014/I.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Abbreviations

NGS, next generation sequencing; WGS, whole genome sequencing; WES, whole exome sequencing; MAF, minor allele frequency; PGKB, Pharmacogenomics Knowledge Base; SNVs, single nucleotide variant; PK/PD, pharmacokinetics/pharmacodynamics.

Acknowledgments

We gratefully acknowledge the procurement of sequencing data by the 100 Genomes Project and the MAIS cohort at INMEGEN Mexico.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2019.01169/full#supplementary-material

References

Auton, A., Abecasis, G. R., Altshuler, D. M., Durbin, R. M., Abecasis, G. R., Bentley, D. R., et al. (2015). A global reference for human genetic variation. Nature 526 (7571), 68–74. doi: 10.1038/nature15393

PubMed Abstract | CrossRef Full Text | Google Scholar

Auwera, G. A. V. D., Carneiro, M. O., Poplin, C. H. R., del Angel, G., Levy-Moonshine, A., Jordan, T., et al. (2013). From FastQ data to high confidence variant calls: the genome analysis toolkit best practices pipeline. Curr. Protoc. Bioinf. 11 (1110), 11.10.1–11.10.33. doi: 10.1002/0471250953.bi1110s43

CrossRef Full Text | Google Scholar

Beekman, M., Blanche, H., Perola, M., Hervonen, A., Bezrukov, V., Sikora, E., et al. (2013). Genome-wide linkage analysis for human longevity: genetics of healthy aging study. Aging Cell 12 (2), 184–193. doi: 10.1111/acel.12039

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonifaz-Pena, V., Contreras, A. V., Struchiner, C. J., Roela, R. A., Furuya-Mazzotti, T. K., Chammas, R., et al. (2014). Exploring the distribution of genetic markers of pharmacogenomics relevance in Brazilian and Mexican Populations. PLoS ONE 9 (11), e112640. doi: 10.1371/journal.pone.0112640

PubMed Abstract | CrossRef Full Text | Google Scholar

Bush, W. S., Crosslin, D. R., Owusu-Obeng, A., Wallace, J., Almoguera, B., Basford, M. A., et al. (2016). Genetic variation among 82 pharmacogenes: the PGRNseq Data from the EMERGE Network. Clin. Pharmacol. Ther. 100 (2), 160–169. doi: 10.1002/cpt.350

PubMed Abstract | CrossRef Full Text | Google Scholar

Cid-Soto, M. A., Martínez-Hernández, A., García-Ortíz, H., Córdova, E. J., Barajas-Olmos, F., Centeno-Cruz, F., et al. (2018). Gene variants in AKT1, GCKR and SOCS3 are differentially associated with metabolic traits in Mexican Amerindians and Mestizos. Gene 679, 160–171 doi: 10.1016/j.gene.2018.08.076

PubMed Abstract | CrossRef Full Text | Google Scholar

Claudio-Campos, K., Duconge, J., Cadilla, C. L., Ruaño, G. (2015). Pharmacogenetics of drug-metabolizing enzymes in US Hispanics. Drug Metab. Pers. Ther. 30 (2), 87–105. doi: 10.1515/dmdi-2014-0023

PubMed Abstract | CrossRef Full Text | Google Scholar

Contreras, A. V., Monge-Cazares, T., Alfaro-Ruiz, L., Hernandez-Morales, S., Miranda-Ortiz, H., Carrillo-Sanchez, K., et al. (2011). Resequencing, haplotype construction and identification of novel variants of CYP2D6 in Mexican Mestizos. Pharmacogenomics 12 (5), 745–756. doi: 10.2217/pgs.11.8

PubMed Abstract | CrossRef Full Text | Google Scholar

Contreras-Cubas, C., Sánchez-Hernández, B. E., García-Ortiz, H., Martínez-Hernández, A., Barajas-Olmos, F., Cid, M., et al. (2016). Heterogenous distribution of MTHFR gene variants among Mestizos and diverse Amerindian groups from Mexico. PLoS ONE. 11 (9), e0163248 doi: 10.1371/journal.pone.0163248

PubMed Abstract | CrossRef Full Text | Google Scholar

Cruz-Correa, O. F., Baltazar Reyes León-Cachón, R., Barrera-Saldaña, H. A., Soberón, X. (2017). Prediction of atorvastatin plasmatic concentrations in healthy volunteers using integrated pharmacogenetics sequencing. Pharmacogenomics 18 (2), 121–131. doi: 10.2217/pgs-2016-0072

PubMed Abstract | CrossRef Full Text | Google Scholar

Cuautle-Rodriguez, P., Llerena, A., Molina-Guarneros, J. (2014). Present status and perspective of pharmacogenetics in Mexico. Drug Metab. Drug Interact. 29, 37. doi: 10.1515/dmdi-2013-0019

CrossRef Full Text | Google Scholar

Dai, D.-P., Hu, L.-M., Geng, P.-W., Wang, S.-H., Cai, J., Hu, G.-X., et al. (2015). In-vitro functional analysis of 24 novel CYP2C19 variants recently found in the Chinese Han population. Xenobiotica 45 (11), 1030–1035. doi: 10.3109/00498254.2015.1028512

PubMed Abstract | CrossRef Full Text | Google Scholar

de Andrés, F., Sosa-Macías, M., Ramos, B. P. L., G. Naranjo, M-, LLerena, A. (2017). CYP450 genotype/phenotype concordance in Mexican Amerindian indigenous populations—where to from here for global precision medicine? OMICS: J. Integr. Biol. 21 (9), 509–519. doi: 10.1089/omi.2017.0101

CrossRef Full Text | Google Scholar

Del Re, M., Michelucci, A., Di Leo, A., Cantore, M., Bordonaro, R., Simi, P., et al. (2015). Discovery of novel mutations in the dihydropyrimidine dehydrogenase gene associated with toxicity of fluoropyrimidines and viewpoint on preemptive pharmacogenetic screening in patients. EPMA J. 6 (1), 17. doi: 10.1186/s13167-015-0039-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Desta, Z., Gammal, R. S., Gong, L., Whirl-Carrillo, M., Gaur, A. H., Sukasem, C., et al. (2019). Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for CYP2B6 and efavirenz-containing antiretroviral therapy. Clin. Pharmacol. Ther. doi: 10.1002/cpt.1477

CrossRef Full Text | Google Scholar

Devarajan, S., Moon, I., Ho, M.-F., Larson, N. B., Neavin, D. R., Moyer, A. M., et al. (2019). Pharmacogenomic next-generation DNA sequencing: lessons from the identification and functional characterization of variants of unknown significance in CYP2C9 and CYP2C19. Drug Metab. Dispos. 47 (4), 425–435. doi: 10.1124/dmd.118.084269

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, C., Wei, P., Jian, X., Gibbs, R., Boerwinkle, E., Wang, K., et al. (2015). Comparison and integration of deleteriousness prediction methods for nonsynonymous SNVs in whole exome sequencing studies. Hum. Mol. Genet. 24 (8), 2125–2137. doi: 10.1093/hmg/ddu733

PubMed Abstract | CrossRef Full Text | Google Scholar

Elraiyah, T., Jerde, C. R., Shrestha, S., Wu, R., Nie, Q., Giama, N. H., et al. (2017). Novel deleterious dihydropyrimidine dehydrogenase variants may contribute to 5-fluorouracil sensitivity in an East African population. Clin. Pharmacol. Ther. 101 (3), 382–390. doi: 10.1002/cpt.531

PubMed Abstract | CrossRef Full Text | Google Scholar

Flannick, J., Mercader, J. M., Fuchsberger, C., Udler, M. S., Mahajan, A., Wessel, J., et al. (2019). Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls. Nature 1. 570 (7759), 71–76. doi: 10.1038/s41586-019-1231-2

CrossRef Full Text | Google Scholar

Fricke-Galindo, I., Jung-Cook, H., LLerena, A., López-López, M., Scott, S., Stein, C. M. (2016). Interethnic variability of pharmacogenetic biomarkers in Mexican healthy volunteers: a report from the RIBEF (Ibero-American Network of Pharmacogenetics and Pharmacogenomics). Drug Metab. Pers. Ther. 31 (2), 625–629. doi: 10.1515/dmpt-2015-0030

CrossRef Full Text | Google Scholar

Fricke-Galindo, I., Martínez-Juárez, I. E., Monroy-Jaramillo, N., Jung-Cook, H., Falfán-Valencia, R., Ortega-Vázquez, A., et al. (2014). HLA-A*02:01:01/-B*35:01:01/-C*04:01:01 haplotype associated with lamotrigine-induced maculopapular exanthema in Mexican Mestizo Patients. Pharmacogenomics 15 (15), 1881–1891. doi: 10.2217/pgs.14.135

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzalez-Covarrubias, V., MartÃnez-MagaÃ, J. J., Coronado-Sosa, R., Villegas-Torres, B., Genis-Mendoza, A. D., Canales-Herrerias, P., et al. (2016). Exploring variation in known pharmacogenetic variants and its association with drug response in different Mexican populations. Pharm. Res. 33 (11), 2644–2652. doi: 10.1007/s11095-016-1990-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzalez-Covarrubias, V., Urena-Carrion, J., Villegas-Torres, B., Cossío-Aranda, J. E., Trevethan-Cravioto, S., Izaguirre-Avila, R., et al. (2017). Pharmacogenetic variation in over 100 genes in patients receiving acenocumarol. Front. Pharmacol. 8 (NOV). 8, 863–870. doi: 10.3389/fphar.2017.00863

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, S. M., Park, J., Lee, J. H., Lee, S. S., Kim, H., Han, H., et al. (2017). Targeted next-generation sequencing for comprehensive genetic profiling of pharmacogenes. Clin. Pharmacol. Ther. 101 (3), 396–405. doi: 10.1002/cpt.532

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, J. A., Cavallari, L. H., Beitelshees, A. L., Lewis, J. P., Shuldiner, A. R., Roden, D. M. (2011). Pharmacogenomics: application to the management of cardiovascular disease. Clin. Pharmacol. Ther. 90 (4), 519–531. doi: 10.1038/clpt.2011.179

PubMed Abstract | CrossRef Full Text | Google Scholar

Karczewski, K. J., Weisburd, B., Thomas, B., Solomonson, M., Ruderfer, D. M., Kavanagh, D., et al. (2017). The ExAC browser: displaying reference data information from over 60 000 exomes. Nucleic Acids Res. 45 (D1), D840–D845. doi: 10.1093/nar/gkw971

PubMed Abstract | CrossRef Full Text | Google Scholar

Katsila, T., Patrinos, G. P. (2015). Whole genome sequencing in pharmacogenomics. Front. Pharmacol. 6 (MAR). 26 (6), 61–69 doi: 10.3389/fphar.2015.00061

CrossRef Full Text | Google Scholar

Lakiotaki, K., Kanterakis, A., Kartsaki, E., Katsila, T., Patrinos, G. P., Potamias, G. (2017). Exploring public genomics data for population pharmacogenomics. Edited by Alvaro. Galli. PLOS ONE 12 (8), e0182138. doi: 10.1371/journal.pone.0182138

CrossRef Full Text | Google Scholar

Lek, M., Karczewski, K. J., Minikel, E. V., Samocha, K. E., Banks, E., Fennell, T., et al. (2016). Analysis of protein-coding genetic variation in 60,706 humans. Nature 536 (7616), 285–291. doi: 10.1038/nature19057

PubMed Abstract | CrossRef Full Text | Google Scholar

Marsh, S., King, C. R., Van Booven, D. J., Revollo, J. Y., Gilman, R. H., McLeod, H. L. (2015). Pharmacogenomic assessment of Mexican and Peruvian populations. Pharmacogenomics 16 (5), 441–448. doi: 10.2217/pgs.15.10

PubMed Abstract | CrossRef Full Text | Google Scholar

Matimba, A., Del-Favero, J., Van Broeckhoven, C., Masimirembwa, C. (2009). Novel variants of major drug-metabolising enzyme genes in diverse African populations and their predicted functional effects. Hum. Genomics 3 (2), 169. doi: 10.1186/1479-7364-3-2-169

PubMed Abstract | CrossRef Full Text | Google Scholar

McInnes, G., Dalton, R., Sangkuhl, K., Whirl-Carrillo, M., Lee, S.-B., Altman, R. B., et al. (2019). Hubble2D6: a deep learning approach for predicting drug metabolic activity. BioRxiv June, 684357. doi: 10.1101/684357

CrossRef Full Text | Google Scholar

Mizzi, C., Peters, B., Mitropoulou, C., Mitropoulos, K., Katsila, T., Agarwal, M. R., et al. (2014). Personalized pharmacogenomics profiling using whole-genome sequencing. Pharmacogenomics 15 (9), 1223–1234. doi: 10.2217/pgs.14.102

PubMed Abstract | CrossRef Full Text | Google Scholar

Moreno-Estrada, A., Gignoux, C. R., Fernendez-Lopez, J. C., Zakharia, F., Sikora, M., Contreras, A. V., et al. (2014). The genetics of Mexico recapitulates Native American substructure and affects biomedical traits. Science 344 (6189), 1280–1285. doi: 10.1126/science.1251688

PubMed Abstract | CrossRef Full Text | Google Scholar

Mukonzo, J. K., Röshammar, D., Waako, P., Andersson, M., Fukasawa, T., Milani, L., et al. (2009). A novel polymorphism in ABCB1 gene, CYP2B6*6 and sex predict single-dose efavirenz population pharmacokinetics in Ugandans. Br. J. Clin. Pharmacol. 68 (5), 690–699. doi: 10.1111/j.1365-2125.2009.03516.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Niioka, T., Satoh, S., Kagaya, H., Numakura, K., Inoue, T., Saito, M., et al. (2012). Comparison of pharmacokinetics and pharmacogenetics of once- and twice-daily tacrolimus in the early stage after renal transplantation. Transplantation 94 (10), 1013–1019. doi: 10.1097/TP.0b013e31826bc400

PubMed Abstract | CrossRef Full Text | Google Scholar

Oetting, W. S., Schladt, D. P., Guan, W., Miller, M. B., Remmel, R. P., Dorr, C., et al. (2016). Genomewide association study of tacrolimus concentrations in African American kidney transplant recipients identifies multiple CYP3A5 alleles. Am. J. Transplant. 16 (2), 574–582 doi: 10.1111/ajt.13495

PubMed Abstract | CrossRef Full Text | Google Scholar

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81 (3), 559–575. doi: 10.1086/519795

PubMed Abstract | CrossRef Full Text | Google Scholar

Relling, M. V., Klein, T. E. (2011). CPIC: clinical pharmacogenetics implementation consortium of the pharmacogenomics research network. Clin. Pharmacol. Ther. 89 (3), 464–467. doi: 10.1038/clpt.2010.279

PubMed Abstract | CrossRef Full Text | Google Scholar

Romero-Hidalgo, S., Ochoa-Leyva, A., Garcíarrubio, A., Acuña-Alonzo, V., Antúnez-Argüelles, E., Balcazar-Quintero, M., et al. (2017). Demographic history and biologically relevant genetic variation of Native Mexicans inferred from whole-genome sequencing. Nat. Commun. 8 (1), 1005. doi: 10.1038/s41467-017-01194-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Sánchez-Pozos, K., Ortíz-López, M. G., Peña-Espinoza, B. I., de Los Ángeles Granados-Silvestre, M., Jiménez-Jacinto, V., Verleyen, J., et al. (2018). Whole-exome sequencing in Maya indigenous families: variant in PPP1R3A is associated with type 2 diabetes. Mol. Genet. Genomics. 293 (5), 1205–1216 doi: 10.1007/s00438-018-1453-2

CrossRef Full Text | Google Scholar

Sangkuhl, K., Berlin, D. S., Altman, R. B., Klein, T. E. (2008). PharmGKB: understanding the effects of individual genetic variants. Drug Metab. Rev. 40. 40 (4), 539–551. doi: 10.1080/03602530802413338

PubMed Abstract | CrossRef Full Text | Google Scholar

Sconce, E. A., Khan, T. I., Wynne, H. A., Avery, P., Monkhouse, L., King, B. P., et al. (2005). The impact of CYP2C9 and VKORC1 genetic polymorphism and patient characteristics upon warfarin dose requirements: proposal for a new dosing regimen. Blood 106 (7), 2329–2333. doi: 10.1182/blood-2005-03-1108

PubMed Abstract | CrossRef Full Text | Google Scholar

Team R Core. (2014). A language and environment for statistical computing, Http://Www.R-Project.Org.”.

Google Scholar

Van Booven, D., Marsh, S., McLeod, H., Carrillo, M. W., Sangkuhl, K., Klein, T. E., et al. (2010). Cytochrome P450 2C9-CYP2C9. Pharmacogenet. Genomics 20 (4), 277–281. doi: 10.1097/FPC.0b013e3283349e84

PubMed Abstract | CrossRef Full Text | Google Scholar

van Driest, S. L., Shi, Y., Bowton, E. A., Schildcrout, J. S., Peterson, J. F., Pulley, J., et al. (2014). Clinically actionable genotypes among 10,000 patients with preemptive pharmacogenomic testing. Clin. Pharmacol. Ther. 95 (4), 423–431. doi: 10.1038/clpt.2013.229

PubMed Abstract | CrossRef Full Text | Google Scholar

Vargas-Alarcón, G., Ramírez-Bello, J., de la Peña, A., Calderón-Cruz, B., Peña-Duque, M. A., Martínez-Ríos, M. A., et al. (2014). Distribution of ABCB1, CYP3A5, CYP2C19, and P2RY12 gene polymorphisms in a Mexican Mestizos population. Mol. Biol. Rep. 41 (10), 7023–7029. doi: 10.1007/s11033-014-3590-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Villegas-Torres, B., Sanchez-Giron, F., Jaramillo-Villafuerte, K., Soberon, X., Gonzalez-Covarrubias, V. (2015). Genotype frequencies of VKORC1 and CYP2C9 in Native and Mestizo populations from Mexico, potential impact for coumarin dosing. Gene 558 (2), 235–240. doi: 10.1016/j.gene.2014.12.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Wendt, F. R., Pathak, G., Sajantila, A., Chakraborty, R., Budowle, B. (2018). Global genetic variation of select opiate metabolism genes in self-reported healthy individuals. Pharmacogenomics J. 18 (2), 281–294. doi: 10.1038/tpj.2017.13

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Lu, Y. L., Wang, H. J., Zhou, W. H. (2016). Pharmacogenomics study of 620 whole-exome sequencing: focusing on aspirin application. Zhonghua Er Ke Za Zhi = Chin. J. Pediatr. 54 (5), 332–336. doi: 10.3760/cma.j.issn.0578-1310.2016.05.005

CrossRef Full Text | Google Scholar

Zhou, Y., Fujikura, K., Mkrtchian, S., Lauschke, V. M. (2018). Computational methods for the pharmacogenetic interpretation of next generation sequencing data. Front. Pharmacol. 9 (December), 1437. doi: 10.3389/fphar.2018.01437

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Ingelman-Sundberg, M., Lauschke, V. M. (2017). Worldwide distribution of cytochrome P450 alleles: a meta-analysis of population-scale sequencing projects. Clin. Pharmacol. Ther. 102 (4), 688–700. doi: 10.1002/cpt.690

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Mkrtchian, S., Kumondai, M., Hiratsuka, M., Lauschke, V. M. (2019). An optimized prediction framework to assess the functional impact of pharmacogenetic variants. Pharmacogenomics J. 19 (2), 115–126. doi: 10.1038/s41397-018-0044-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, H. J., Yuan, S. H., Fang, Y., Sun, X. Z., Kong, H., Ge, W. H. (2011). The effect of CYP3A5 polymorphism on dose-adjusted cyclosporine concentration in renal transplant recipients: a meta-analysis. Pharmacogenomics J. 11 (3), 237–246. doi: 10.1038/tpj.2010.26

PubMed Abstract | CrossRef Full Text | Google Scholar