1. Introduction
Adverse drug reaction (ADR) may be defined as “any response to a drug which is noxious and unintended, and which occurs at doses normally used in man for prophylaxis, diagnosis, or therapy of disease, or for the modification of physiological function” [1].
ADRs are important cause of patient morbidity and mortality, as they are responsible for increased health costs and lead to lack of patient adherence to medical treatment worldwide. It is estimated they may occur in over 30–60% of drug prescriptions [2].
There are several factors that trigger the appearance of ADRs in different types of treatment schemes, including patient age, disease status, physiological conditions, drug dosage, and drug–drug interaction [3]. Additionally, interindividual genetic variability has also being demonstrated to play an essential role to this condition [4,5].
Pharmacogenomics—the study of genetic variants that may influence drug efficacy, response, and/or toxicity—accounts for around 80% of the variability in drug pharmacokinetics and pharmacodynamics and has been associated with over 60% of ADRs events nowadays [2]. Given this close relationship, an assessment about the pharmacogenomic profile of patients prior to drug prescription or dose adjustment is a pivotal strategy.
The goal of pharmacogenomics is the establishment of precision medicine, increasing drug efficacy and safety while minimizing events of resistance, toxicity, and adverse effects associated to drugs [6]. To achieve this purpose, it is necessary to obtain accurate information about actionable pharmacogenetic variants. The interethnic variability is an important issue to be considered in pharmacogenomic studies since there is an enormous variability in the distribution and frequency of pharmacogenetics variants among different populations, with most of them being rare or specific to a given biogeographical population [3,7,8].
Most of the information acquired in the pharmacogenomic field, especially that employing large-scale sequencing data and specific pharmacogenomic protocols, has been centered on more homogeneous populations, specifically those of European origins [9,10], which cannot be fully applicable to other ethnic groups
Genetically distinct populations need specific algorithms for drug dosing. This is especially relevant when we consider Native American populations, which have a unique genetic profile due to a long process of geographic isolation and other genetic and evolutionary processes, such as genetic drift, founding effect, natural selection, and inbreeding [11].
These populations are extremely underrepresented in genetic research. In parallel, some studies have already demonstrated that Amerindian groups have distinct genetic profiles concerning pharmacogenomic biomarkers widely used in clinical practice [12,13], which further reinforces the need for their inclusion in pharmacogenomic studies.
Therefore, the purpose of this study is to describe the pharmacogenetic variability of Native American populations regarding 160 pharmacogenes involved in absorption, distribution, metabolism, and excretion (ADME) processes and biological pathways of different therapies based on data obtained through complete exome sequencing of 64 individuals from different Amerindian groups of the Brazilian Amazon.
2. Materials and Methods
2.1. Study Population
The study population is composed of 64 Amerindians. The Amerindians represent 12 different Amazonian ethnic groups that were grouped together as the Native American (NAM) group. Ancestry analysis was performed as described by Ramos et al. 2016 [14], using 61 autosomal ancestry informative markers (AIMs). The individual proportions of European, African, and Native Americans’ genetic ancestries in the NAM group were 0.022, 0.014, and 0.964, respectively. Details such as name, location, and number of individuals in each ethnic group are described in Supplementary Table S1. The present study was approved by the National Committee for Ethics in Research (CONEP), identified by Nos. 1062/2006 and 123/98, and the Research Ethics Committee of the UFPA Tropical Medicine Center, under CAAE number 20654313.6.0000.5172. All participants signed a free-informed consent as well as the tribe leaders when necessary. We compared our results with those of populations from other continents obtained from the Exome Aggregation Consortium (ExAC) (available at
2.2. Exome Library Preparation
The exome library preparation of the initial 58 individuals of the NAM group was performed as previously described in [15]. Afterward, six additional individuals were added in this group set, and the exome library preparation is described in [16].
2.3. Read Calling and Processing
The quality of the FASTQ reads was analyzed (FastQC v.0.11—
2.4. Selection of PGx Biomarkers
From the generated exome data, 160 pharmacogenes associated with the ADME process or the pharmacodynamic changes of different drug classes were selected for analysis (Supplementary Table S2). The genes were chosen because they are classified as very important pharmacogenes (VIPs) genes by the PharmGKB database [17] and/or because they are described with clinical significance by the Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines [18]. Additionally, some of these markers are already recommended by drug regulatory agencies as predictors of toxicity and response to therapeutic regimens routinely used in clinical practice.
2.5. Bioinformatic and Statistical Analyses
The analysis of the variant annotations was run in the ViVa® (Viewer of Variants) software developed by the Federal University of Rio Grande do Norte (UFRN) bioinformatics team. ViVa® performs genomic annotation for each variant, including chromosomal location and position, reference and variant allele(s), and quality control analysis of nucleotide identification according to Phred score (Q score) (ILLUMINA, 2011) and sample coverage. According to the different National Center for Biotechnology Information (NCBI) databases (
Additionally, annotations also include variant impact according to SnpEFF (
The allele frequencies of the NAM populations were obtained directly by gene counting and compared with the other study populations (AFR, EUR, AMR, EAS, and SAS). The difference in frequencies between the populations were analyzed by Fisher’s exact test, and results were considered significant when Bonferroni-adjusted p-value ≤ 0.00015. The principal component analysis (PCA) was performed using the FactoMineR, version 2.4 package [19]. For all analyses, R package v. 4.1.3 [20] and Python 3.10.4 [21] were used.
3. Results
Of the total of 160 genes analyzed, our investigation found a total of 3311 variants in the study subjects. Of this total, 167 are exclusive to Amerindian populations. The general characteristics of the variants are described in Table 1.
The great majority are single-nucleotide variants (SNV). There are six gene regions where the variants were located; most of them were in intronic regions (1567), followed by exonic regions, specifically in the coding sequence DNA (CDS) (1183). Regarding the in silico functional impact prediction, most of the variants received modifier status, which is defined as “generally non-coding variants or variants that affect non-coding genes, where predictions are difficult or there is no evidence of impact” [22].
Among the exclusive variants found in the Amerindian group, we can highlight 21 that presented a population frequency above 0.01 (Table 2). Of these, two are missense-like variants, which have a moderate impact: the variant at position 19:39248563 (MAF = 0.0161) found in the IFNL4 gene and the variant at position 1:97098596 (MAF = 0.0313) in the DPYD gene. In addition, we could also find seven variants with high allelic frequencies (MAF > 0.5) in the Amerindian populations: two in the GSTT1 gene, three in the CYP2D6 gene, one in the MTHFR gene, and one in the TYMS gene. All seven have a modifying impact and are present in intronic regions or 5’UTR regions, which are known to regulate gene expression.
Due to the large number of variants found, we narrowed down variants based in three main criteria: (i) a minimum of 10 reads of coverage; (ii) those which show a high (a disruptive impact in the protein, probably causing protein truncation, loss of function, or triggering nonsense mediated decay) or moderate (a non-disruptive variant that might change protein effectiveness) impact; and (iii) variants that show allelic frequency >1% in at least of two other continental populations. Therefore, in the posterior analyses, the results are based on 332 variants (Supplementary Table S3).
Plotting a principal component analysis based on the allelic frequencies of the 332 variants found in each of the six populations (AFR, AMR, EAS, EUR, NAM, and SAS) could explain nearly 66% of the total variation found in our group of variants (Figure 1).
The NAM population formed a cluster isolated from all other study populations, which allows us to corroborate the differentiated genetic profile of this group concerning the variants analyzed in this study. Furthermore, we can observe that the population with the most similar genetic profile of the NAM population is the AMR population, composed of Latin American population groups, followed by the EAS population, formed by population groups from East Asia.
Additionally, the 332 pharmacogenomic variants were evaluated in pairwise comparisons between the NAM population and each of the Exome Aggregation Consortium (ExAC) population, individually, using Fisher’s exact test (Supplementary Table S4). The largest number of variants with significant differences was found in the comparison between NAM and AFR (144 variants in total), followed by the comparison NAM × EUR (133). The smallest number of variants with significant differences was observed between the comparison of NAM and AMR (75). The number of variants with significant differences between NAM and Asian populations was basically the same (114 compared with EAS and 115 with SAS). These data agree with that shown in Figure 1. Table 3 shows the 42 pharmacogenetic variants that showed significant differences between NAM and all EXAC populations.
Of these, ten variants are classified with a high impact: rs881711 of the CBR3 gene; rs200088269rs35616319, rs753045406, and rs757139064 of the HLA-DRB1 gene; rs199556640 and rs9282026 of the HLA-DQA1 gene; rs74597329 and rs11322783 of the IFNL4 gene; rs3745274 of the CYP2B6 gene; and finally, rs149012039 of the CYP2D6 gene. Additionally, most of variants present higher allele frequencies in the group of Amazonian Amerindians compared to the other world populations except for nineteen variants that present the opposite profile (rs881711, rs1050147, rs130066, rs149012039, rs11322783, rs2032582, rs396991, rs3745274, rs4963, rs707952, rs2308466, rs2523600, rs1042713, rs713031, rs697742, rs1143146, rs1131170, rs2228001, and rs2228570).
4. Discussion
The incorporation of pharmacogenomic assessment prior to treatment would be of great benefit in terms of costs, quality of life, and optimization of therapeutic resources. Given this, one of the main goals in the future of pharmacogenomic research is characterization of gene polymorphisms to elucidate the genetic background underlying differences in drug responses. However, some population groups remain underrepresented in these types of investigations, such as Native American populations. Recently, some efforts have been made to include genetic information about allelic frequencies of Native American populations and their influence on susceptibility to diseases, such as severe forms of COVID-19 [23], psychological disorders [24], oncological diseases [25,26], and also on different types of pharmacological treatments [12,13,27,28]. Most of these studies are focused on one or a few group of genes. Hence, there is still a paucity of large-genome-scale studies on Native American populations, particularly Amazonian Amerindian groups, regarding the response pattern and toxicity of drugs that are known to have pharmacogenetic influences. The focus of this work was to evaluate genes of pharmacogenomic importance as recommended by global drug regulatory institutions, such as the Food and Drug Administration from United States, the European Medicine Agency from Europe, and the Clinical Pharmacogenetics Implementation Consortium (CPIC), in a group of 64 Amerindians from the Amazon region of Brazil. This is the first study to analyze a large set of guided-treatment pharmacogenes by whole-exome sequencing in Amazonian Amerindian populations.
4.1. General Data Found in the Amazonian Amerindian Population
The exome sequencing of the 160 pharmacogenes analyzed revealed more than 3000 variants. Interestingly, most variants that have already been defined as a biomarker-definers of drug indication, such as DPYD*2A, DPYD *13, rs67376798, or rs75017182 of the DPYD [29] and the rs9923231 of the VKORC1 genes [30], were not found in the set of Amerindians of this study. An exception is the rs1272632214, which was previously reported [16] with a high frequency in the Amerindians of this study and in complete linkage disequilibrium with the rs116855232 of the NUDT15 gene, which is a variant linked to thiopurine toxic, included in the recommendations of CPIC guidelines for adjusting starting doses of azathioprine, mercaptopurine, and thioguanine [31].
An important finding of this is study is the large number of novel variants firstly reported in the Amazon Amerindians groups (here, they will be referred by their chromosomal location). Most variants were SNVs, and one-third of them were located in gene coding sequences. Of the 167 novel variants reported, nine must be highlighted: (1) the 19:39248563 of the IFNL4 gene; (2) the 1:97098596 in the DPYD gene; (3) the 22_KI270879v1_alt:278462 and 22_KI270879v1_alt:278129 of the GSTT1; (4) the 1:11803345 in MTHFR; (5) the 18:657685 in TYMS; and (6) the 22_KI270928v1_alt:52910, 22_KI270928v1_alt:52912, and 22_KI270928v1_alt:52919 of the CYP2D6 gene. To validate the results reported here, these variants must be confirmed in studies with a larger number of individuals and analyzed regarding their functional effect. Here, we will discuss the role of these genes to the pharmacogenomic field.
4.1.1. IFNL4
The IFNL4, along with three neighboring genes of the interferon lamba family (IFNL1, IFNL2, and IFNL3) leads to activation of the JAK-STAT signaling pathway and upregulation of numerous interferon-stimulated genes [32]. This gene has been included in the CPIC guidelines in association with peginterferon alfa-2a and peginterferon alfa-2b, two regimen treatments with direct-acting antivirals approved for hepatitis C virus (HCV) genotype 1 infection [33]. Variants in the IFNL4 have shown to be a better predictor of response to PegIFN in patients with either African or European ancestry [34]. Because there is a strong relationship with sustained virologic response to HCV treatments, a polymorphism (rs12979860) of the IFNL4 gene is included as an actionable variant in drug labels for sofosbuvir, ombitasvir, paritaprevir, and ritonavir by the Swissmedic agency (
4.1.2. DPYD
The DPYD gene encodes the dihydropyrimidine dehydrogenase (DPD), a critical enzyme responsible for catabolism of approximately 85% of the 5-fluorouracil (5-FU) to its inactive metabolic form [35]. The 5-flurouracil is the active agent of fluoropyrimidine chemotherapy, frequently prescribed for treatment of a variety of cancers, including colorectal, upper gastrointestinal, breast, and head and neck cancers. Despite their large use in clinical practice, the fluoropyrimidine treatment has a narrow therapeutic index; even at standard dose, a substantial proportion of patients may develop grade 3 or higher toxicities [36]. The guidelines of the Clinical Pharmacogenetics Implementation Consortium (CPIC) describe four decreased-function DPYD variants, namely rs3918290 (also known as DPYD*2A) rs55886062 (DPYD *13), rs67376798 (p.D949V), and rs75017182 (HapB3), and strongly recommend the use of alternative drugs or the reduction (25% to 50%) of the standard dose of fluoropyrimidines for patients who are classified as DPYD intermediate or poor metabolizer, respectively [29]. A study evaluating Amerindian populations from the Amazon reported a total of nine polymorphic variants of the DPYD; two of them are variants responsible for decreasing the DPD enzymatic activity (included in the CPIC guidelines), the rs3918290 and rs55886062, with 1% and 2% of frequency, respectively; and three of them (rs17116806, rs1801159, and rs4970722) had frequencies higher than 40%, contrasting with frequencies much lower found in populations with African, European, or Asian ancestry; these variants show lower frequencies [12].
4.1.3. GSTT1
The GSTT1 belong to the super family of glutathione S-transferases. These enzymes are implicated in the second-phase metabolism of several compounds, such as anticancer and xenobiotic substances [37]. Polymorphisms in the GSTT1 gene show significant frequency differences among distinct ethnic groups: East Asians present the highest frequencies worldwide for GSTT1 deletions [38], which could implicate in higher interindividual variability in pharmacotherapy responses and susceptibility to various diseases, especially those of hepatological disorders [39]. The results regarding the role of GSTT1 conflict depending on the type of treatment evaluated. Regarding breast cancer therapy, the GSTT1 double-null genotype has been significantly associated with an increased tumor response [40], while for patients with chronic myeloid leukemia, this genotype has been linked to worse treatment outcome to imatinib [41].
4.1.4. MTHFR
The MTHFR gene encodes the 5,10-methylenetetrahydrofolates enzyme, a substrate for de novo purine synthesis and therefore crucial to produce amino acids. Additionally, this enzyme is involved in the conversion of homocysteine to methionine, which is then converted to the universal methyl donor, S-adenosylmethionine, used for methylation of DNA and proteins [42]. Given the role of MTHFR in both DNA synthesis and methylation status, this gene participates of several pathways commonly used as targets for chemotherapeutic antineoplastic and antirheumatic drugs, such as 5-fluorouracil and methotrexate [43,44]. The MTHFR gene has two main variants: the rs1801133 and rs1801131. A recent meta-analysis evaluating 34 studies showed that both variants are statistically associated with toxicities during high-dose methotrexate treatment; the rs1801133 was positively associated with increased risk of hepatotoxicity, mucositis, and renal toxicity, while the rs1801131 showed a decreased risk of renal toxicity during the therapy [45]. Concerning fluoropyrimidine-based treatments, results involving these variants are conflicting: while Campbell and colleagues found significant associations of rs1801133 and rs1801131 with a higher frequency of nonhematologic toxicity (nausea/vomiting), a protective effect for neutropenia, and global toxicity [46], Zhong and collaborators found no significant association between both polymorphisms and the clinical response to fluoropyrimidine-based chemotherapy under any of the three genetic models (allele model, dominant model, and recessive model) [47]. Curiously, a study evaluating patients with gastrointestinal or colorectal cancer from the Brazilian Amazon region in which the study population showed means of 31.1% of Amerindian ancestry, reported a significant association with the rs1801133 in the MTHFR and severe toxicity during fluoropyridimine-based treatment [48].
4.1.5. TYMS
The TYMS gene encodes the thymidylate synthase, the enzyme responsible for catalyzing the methylation of deoxyuridylate (dUMP) to deoxythymidylate (dTMP). Since the dTMP is the only source of intracellular thymidine, the TS is a key enzyme for the production of DNA and the principal target of 5-fluorouracil [49]. The expression levels of TS are crucial to determine the effectiveness of 5-FU because the inhibition of TS can disrupt normal DNA synthesis, leading tumor cells to apoptosis. Studies have shown the downregulation of TYMS contributes to 5-FU sensitivity [50], while the overexpression of TS represents a pathway of tumor cells resistance to 5-FU [51]. One of the main reasons for increased TS levels in cancer cells is polymorphisms in the TYMS gene, particularly a triple tandem repeat (TSER *3) detected in the 5′-UTR of the gene; also, TS copy number and different genomic instability statuses, such as chromosome and microsatellite instability, may influence 5-FU responsiveness [52]. The FDA-approved starting dosage for capecitabine monotherapy (a prodrug of 5-FU) is 1250 mg/m2 b.i.d on a 2 weeks on/1 week off regimen; however, a clinical trial conducted by Soo and colleagues suggested that patients with the TSER 3R/3R genotype exhibit enhanced tolerance to capecitabine and may potentially benefit from a higher dose (1500 mg/m2) compared to the standard protocol of 1250 mg/m2 [53]. Another study showed that the overall survival of patients with high TS levels observed in primary tumors with metastasis and in those with lymph node metastasis is shorter than in patients with low TS levels [52]. Together, these findings suggest that the TYMS genotyping assessment may predict efficacy and tolerance of 5-FU-based treatment.
4.1.6. CYP2D6
The cytochrome P450 2D6 (CYP2D6) is responsible for metabolizing approximately 25% of all drugs in the human liver [54,55]. Typical CYP2D6 substrates are lipophilic bases, including some anticancer drugs, antidepressants, antipsychotics, antiarrhythmics, antiemetics, beta-adrenergic receptor antagonists (beta-blockers), and opioids. At least 160 therapeutic targets are known to be metabolized by the CYP2D6 enzyme. Of all the pharmacogenetic biomarkers categorized by the FDA, there are 69 drugs from different pharmacological classes that have relevant pharmacogenomic information regarding the CYP2D6 gene. The information ranges from simple descriptions of clinical pharmacology and drug interactions to even possible dose adjustments depending on the CYP2D6 genotype. The CPIC also formulated six specific guidelines for genetic testing regarding possible dose adjustments based on the metabolism profile of CYP2D6, which include: (i) anesthetics [56], specifically codeine, which is bioactivated in morphine by CYP2D6; (ii) selective serotonin reuptake inhibitor drugs [57], which are the main treatment options for depressive and anxiety disorders; (iii) tricyclic antidepressants [58], for example, nortriptyline, desipramine, etc.; (iv) 5-hydroxytryptamine type 3 (5-HT3) receptor antagonists [59] used for the prevention of nausea and vomiting induced by chemotherapy, radiation, and postoperative care; (v) tamoxifen, a selective estrogen receptor modulator [60] used against breast cancer; and (vi) atomoxetine [61], a non-stimulant medication used to treat attention-deficit/hyperactivity.
4.2. Comparative Analyses between Variants Found in the Amazon Amerindians and the Five Continental Populations from ExAC Database
The origin and expansion of Amerindian populations in the Amazon have been the subject of numerous academic discussions and still remain challenging. Historical evidence indicates that the primitive populations that arrived in the Amazon derive from an ancestral population that originally left Asia by way of the Bering Strait more than 15,000 years ago [62,63]. Other analyses have also shown that some groups in South America, particularly in Brazil, share alleles with indigenous New Guineans, Australians, and Andaman Islanders, which together form a group with a distinct genomic profile known as Australasians [64,65,66].
Recently, Ribeiro-dos Santos et al., 2020, contributed to the reconstruction of the genetic history of Native Brazilian populations, performing whole-exome sequencing of the 58 Amazonian Native Americans included in this study and confirming an occupation model with separate migration waves in the region [15]. These data confirm the unique Amazonian genetic signature in the Amerindian populations, which may also be observed in the PCA generated in this study through the complete segregation of the Amazonian Amerindians from the other populations.
The most closely plotted groups of Amazonian Amerindians were individuals from the AMR group that share origin branches with Native Americans and individuals from East Asia, the ancestral group that engendered the America continent settlement. Still, the PCA also shows an isolated cluster formed by the African population. These findings agree with the history of human evolution: Amerindians are more genetically similar to American and East Asian populations and more genetically distant from African populations [67,68]. This pattern of genetic profile between Amazonian Amerindians and other global population corroborates with other genetic association studies [12,24,25].
5. Conclusions
Our study analyzed genes of great importance in the pharmacogenomics of different types of therapies, based in data generated from whole-exome sequencing of a group of Amazonian Amerindians from northern Brazil. Our results identified many novel variants located in essential genes to cancer (DPYD, TYMS, MTHFT, and GSTT1), infectious diseases (IFNL4), and psychiatric treatments (CYP2D6). This study also demonstrates that the Amazonian Amerindians show very distinct allelic frequencies of pharmacogenomic variants when compared to populations of African, European, Asians, and even Latin origins, which could be observed by statistical and clustering comparative analysis.
These data reinforce the distinct genetic profile concerning biomarkers widely used in clinical practice in Amazonian Amerindians, which further reinforces the need for their inclusion in pharmacogenomic research. Studies with a larger number of individuals and assessment of functional impact of the variants unraveled here are needed to confirm our results. Together, these findings may serve as the scientific basis for the future design of specific treatment protocols for Amazonian Amerindian populations as well as for populations admixed with them, such as the northern Brazilian population.
Conceptualization, J.C.G.R., N.P.C.d.S. and S.S.; methodology, J.C.G.R., G.S.d.A. and A.M.R.-d.-S.; software, S.J.d.S.; formal analysis, J.C.G.R.; resources, M.R.F., Â.R.-d.-S. and S.S.; data curation, J.F.G.; writing—original draft preparation, J.C.G.R.; writing—review and editing, J.C.G.R., N.P.C.d.S. and S.S.; supervision, N.P.C.d.S. and S.S.; project administration, A.M.R.-d.-S., Â.R.-d.-S., M.R.F., N.P.C.d.S. and S.S.; funding acquisition, Â.R.-d.-S., P.P.d.A., N.P.C.d.S. and S.S. All authors have read and agreed to the published version of the manuscript.
This study was conducted according to the guidelines of the Declaration of Helsinki and approved by the National Committee for Ethics in Research (CONEP), identified by Nos. 1062/2006 and 123/98, and the Research Ethics Committee of the UFPA Tropical Medicine Center, under CAAE number 20654313.6.0000.5172.
All participants signed a free-informed consent as well as the ethnic group leaders, when necessary.
The data obtained from the public domain are available at gnomAD (
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Enlarge this image.Figure 1. Principal component analysis of the 332 pharmacogenetic variants analyzed in the six populations of the present study. AFR, African; AMR, American; EAS, East Asian; EUR, European; NAM, Amazonian Amerindian; SAS, South Asian.
Description of all variants found in the exome of 64 Amazonian Amerindians of the study.
| Characteristics | Number |
|---|---|
| All variants | 3.311 |
| New variants | 167 |
| Type | |
| Single-nucleotide variant (SNV) | 2.934 |
| Insertion or deletion (INDEL) | 377 |
| Gene region | |
| 3′UTR | 188 |
| 5′UTR | 90 |
| Exon (CDS) | 1183 |
| Intragenic | 33 |
| Intron | 1567 |
| Others | 250 |
| Consequence | |
| Codon change and codon deletion | 10 |
| Codon change and codon deletion + splice site region | 1 |
| Codon change and codon insertion | 2 |
| Codon deletion | 5 |
| Codon deletion + splice site region | 1 |
| Codon Insertion | 2 |
| Frameshift | 53 |
| Frameshift + splice site region + splice site donor | 1 |
| Intragenic | 33 |
| Intronic | 1567 |
| Protein functional site | 82 |
| Nonsynonymous | 641 |
| Protein structural interaction locus | 30 |
| Splice site acceptor | 3 |
| Splice site donor | 8 |
| Splice region | 102 |
| Start codon | 1 |
| Stop codon | 5 |
| Synonymous | 486 |
| 3′UTR | 188 |
| 5′UTR | 90 |
| SNPeff Impact | |
| High | 101 |
| Moderate | 672 |
| Modifier | 1878 |
| Low | 660 |
Novel variants reported in the Amazonian Amerindians with minor allele frequency greater than 1%.
| Chromosomal Position | Reference | Variant | Gene | Region Detailed | Impact | Minor Allele Frequency |
|---|---|---|---|---|---|---|
| 19:39248563 | C | G | IFNL4 | NON_SYNONYMOUS_CODING | MODERATE | 0.0161 |
| X:134500131 | G | T | HPRT1 | UTR_3_PRIME | MODIFIER | 0.0161 |
| 19:38442623 | G | T | RYR1 | INTRON | MODIFIER | 0.0161 |
| X:38352846 | G | C | OTC | INTRON | MODIFIER | 0.0161 |
| 1:186679301 | G | A | PTGS2 | INTRON | MODIFIER | 0.0161 |
| 7:99778057 | T | C | CYP3A4 | SYNONYMOUS_CODING | LOW | 0.0172 |
| 4:20850609 | C | T | KCNIP4 | SYNONYMOUS_CODING | LOW | 0.0172 |
| 15:74754842 | C | A | CYP1A2 | SYNONYMOUS_CODING | LOW | 0.0172 |
| X:154534297 | AC | A | G6PD | INTRON | MODIFIER | 0.0208 |
| 10:112951639 | TGCCC | T | TCF7L2 | INTRON | MODIFIER | 0.0250 |
| 9:130855098 | TAGGGG | T | ABL1 | SPLICE_SITE_REGION + INTRON | LOW | 0.0250 |
| 1:97098596 | T | C | DPYD | NON_SYNONYMOUS_CODING | MODERATE | 0.0313 |
| 2:233760783 | C | T | UGT1A4 | INTRON | MODIFIER | 0.0357 |
| 2:233760783 | C | T | UGT1A1 | SYNONYMOUS_CODING | LOW | 0.0357 |
| 22_KI270879v1_alt:278462 | C | G | GSTT1 | UTR_5_PRIME | MODIFIER | 0.0500 |
| 1:11803345 | C | A | MTHFR | INTRON | MODIFIER | 0.0556 |
| 18:657685 | G | GGCCTGCCTCCGTCCCGCCGCGCCACTTC | TYMS | UTR_5_PRIME | MODIFIER | 0.1154 |
| 22_KI270928v1_alt:52910 | C | G | CYP2D6 | INTRON | MODIFIER | 0.5400 |
| 22_KI270879v1_alt:278129 | C | T | GSTT1 | INTRON | MODIFIER | 0.5417 |
| 22_KI270928v1_alt:52912 | T | G | CYP2D6 | INTRON | MODIFIER | 0.6230 |
| 22_KI270928v1_alt:52919 | G | C | CYP2D6 | INTRON | MODIFIER | 0.6563 |
Pairwise comparison of variants with a statistically significant difference of allele frequency between Amazonian Amerindians and all continental populations from the ExAC database.
| Pairwise Comparation (p-Value) | |||||||
|---|---|---|---|---|---|---|---|
| Chromosomal |
Gene | dbSNP | NAM × AFR | NAM × AMR | NAM × EAS | NAM × EUR | NAM × SAS |
| chr7:87531302 | ABCB1 | rs2032582 | 2.18306 × 10−56 | 6.53986 × 10−16 | 2.51344 × 10−12 | 1.02973 × 10−16 | 1.00647 × 10−7 |
| chr4:2915035 | ADD1 | rs4963 | 9.45593 × 10−5 | 2.69477 × 10−6 | 4.94438 × 10−17 | 3.89984 × 10−5 | 3.67996 × 10−5 |
| 10:114045297 | ADRB1 | rs1801253 | 3.97318 × 10−15 | 0.000145471 | 2.41152 × 10−8 | 1.2311 × 10−9 | 4.82047 × 10−8 |
| chr5:148826877 | ADRB2 | rs1042713 | 1.19482 × 10−17 | 6.47618 × 10−14 | 1.1541 × 10−20 | 3.6564 × 10−12 | 8.98985 × 10−16 |
| chr11:113400106 | ANKK1 | rs1800497 | 8.91338 × 10−13 | 1.89293 × 10−7 | 7.07468 × 10−10 | 1.75701 × 10−25 | 8.89022 × 10−17 |
| chr11:113396099 | ANKK1 | rs7118900 | 1.57472 × 10−13 | 4.5083 × 10−7 | 1.35542 × 10−9 | 2.61132 × 10−26 | 2.51963 × 10−16 |
| chr22:23285051 | BCR | rs12484731 | 2.18658 × 10−26 | 1.16413 × 10−6 | 9.95701 × 10−26 | 5.38079 × 10−46 | 2.00229 × 10−35 |
| chr22:23290360 | BCR | rs35537221 | 1.8611 × 10−28 | 1.05573 × 10−6 | 4.00474 × 10−21 | 2.14428 × 10−45 | 1.21405 × 10−35 |
| chr21:36135447 | CBR3 | rs881711 | 5.8767 × 10−19 | 2.17402 × 10−19 | 4.67622 × 10−28 | 2.69704 × 10−29 | 1.15895 × 10−27 |
| chr6:31154538 | CCHCR1 | rs130066 | 6.77941 × 10−5 | 1.33884 × 10−6 | 0.000115616 | 1.58769 × 10−6 | 4.19502 × 10−5 |
| chr19:41006936 | CYP2B6 | rs3745274 | 1.64912 × 10−11 | 8.995 × 10−10 | 5.73451 × 10−5 | 9.859 × 10−7 | 1.34727 × 10−12 |
| chr22:42142513 | CYP2D6 | rs149012039 | 1.14215 × 10−11 | 1.59914 × 10−6 | 6.78214 × 10−20 | 3.22875 × 10−9 | 7.70391 × 10−5 |
| chr1:161544752 | FCGR3A | rs396991 | 7.44824 × 10−9 | 1.65936 × 10−5 | 7.84237 × 10−9 | 1.16585 × 10−9 | 2.03133 × 10−9 |
| chr6:29943337 | HLA-A | rs3173420 | 4.31224 × 10−9 | 6.87922 × 10−8 | 4.56562 × 10−7 | 8.57084 × 10−12 | 8.71629 × 10−12 |
| chr6:29942965 | HLA-A | rs281864739rs78306866 | 8.46059 × 10−12 | 4.05707 × 10−17 | 8.97765 × 10−12 | 3.63056 × 10−14 | 4.43796 × 10−13 |
| chr6:29942953 | HLA-A | rs199474436 | 5.79981 × 10−6 | 6.93234 × 10−13 | 1.41018 × 10−11 | 6.78811 × 10−9 | 6.08387 × 10−6 |
| chr6:29942581 | HLA-A | rs1143146 | 2.45192 × 10−16 | 3.07371 × 10−24 | 1.35824 × 10−18 | 1.34049 × 10−18 | 9.97823 × 10−14 |
| chr6:31356226 | HLA-B | rs2308466 | 3.65704 × 10−9 | 3.66621 × 10−9 | 5.81378 × 10−5 | 2.5862 × 10−10 | 7.79316 × 10−12 |
| chr6:31356227 | HLA-B | rs2523600 | 1.9236 × 10−13 | 2.13422 × 10−14 | 8.63059 × 10−15 | 4.81953 × 10−11 | 7.03642 × 10−13 |
| chr6:31356889 | HLA-B | rs713031 | 1.07572 × 10−17 | 1.41509 × 10−16 | 6.48134 × 10−29 | 2.226 × 10−13 | 1.45874 × 10−21 |
| chr6:31356247 | HLA-B | rs697742 | 1.16978 × 10−17 | 4.25927 × 10−20 | 3.31496 × 10−24 | 2.21922 × 10−16 | 3.51264 × 10−21 |
| chr6:31356928 | HLA-B | rs1131170 | 2.28835 × 10−24 | 1.15345 × 10−18 | 2.11729 × 10−29 | 4.93973 × 10−20 | 6.05278 × 10−27 |
| chr6:31269347 | HLA-C | rs1130838 | 6.45396 × 10−16 | 2.25074 × 10−12 | 1.6135 × 10−15 | 4.15458 × 10−9 | 5.96131 × 10−12 |
| chr6:31271999 | HLA-C | rs2074493 | 8.10462 × 10−10 | 9.63528 × 10−11 | 1.02002 × 10−10 | 9.09483 × 10−15 | 4.07089 × 10−20 |
| chr6:31270025 | HLA-C | rs1050147 | 1.18398 × 10−31 | 3.7936 × 10−27 | 7.42085 × 10−32 | 7.33005 × 10−23 | 1.50973 × 10−26 |
| chr6:33080851 | HLA-DPB1 | rs1042136 | 1.6523 × 10−8 | 1.88517 × 10−7 | 2.47857 × 10−10 | 5.23217 × 10−5 | 1.23339 × 10−6 |
| chr6:32641487 | HLA-DQA1 | rs1142333 | 4.373 × 10−28 | 1.9016 × 10−20 | 3.07937 × 10−28 | 1.11648 × 10−30 | 4.88343 × 10−43 |
| chr6:32641535 | HLA-DQA1 | rs1129808 | 7.97358 × 10−92 | 1.19523 × 10−84 | 8.96057 × 10−84 | 8.9906 × 10−104 | 3.2328 × 10−109 |
| chr6:32641477 | HLA-DQA1 | rs1142331 | 1.18138 × 10−25 | 7.36406 × 10−16 | 2.39546 × 10−25 | 2.60828 × 10−27 | 1.4441 × 10−39 |
| chr6:32641521 | HLA-DQA1 | rs199556640 | 9.45834 × 10−19 | 2.68598 × 10−19 | 3.99666 × 10−22 | 6.58968 × 10−19 | 1.29641 × 10−24 |
| chr6:32641519 | HLA-DQA1 | rs9282026 | 6.76015 × 10−12 | 5.23003 × 10−12 | 1.15636 × 10−14 | 8.59051 × 10−12 | 3.46752 × 10−16 |
| chr6:32642029 | HLA-DQA1 | rs707952 | 7.78654 × 10−6 | 4.53726 × 10−8 | 9.34217 × 10−5 | 3.23443 × 10−7 | 2.9079 × 10−7 |
| chr6:32581807 | HLA-DRB1 | rs200088269rs35616319 | 7.84356 × 10−15 | 1.98563 × 10−8 | 7.60999 × 10−12 | 5.08348 × 10−10 | 2.14344 × 10−10 |
| chr6:32581802 | HLA-DRB1 | rs753045406 | 1.2253 × 10−13 | 9.15172 × 10−8 | 4.74343 × 10−11 | 3.01143 × 10−9 | 1.4626 × 10−9 |
| chr6:32581621 | HLA-DRB1 | rs757139064 | 4.04995 × 10−18 | 6.91184 × 10−18 | 2.46801 × 10−16 | 2.08408 × 10−16 | 7.75881 × 10−20 |
| chr19:39248515 | IFNL4 | rs74597329 | 9.56326 × 10−33 | 4.95051 × 10−6 | 1.16957 × 10−5 | 4.53584 × 10−15 | 1.72728 × 10−6 |
| chr19:39248513 | IFNL4 | rs11322783 | 9.56326 × 10−33 | 4.95051 × 10−6 | 1.16957 × 10−5 | 4.53584 × 10−15 | 1.72728 × 10−6 |
| chr6:160540105 | LPA | rs3798220 | 1.75821 × 10−48 | 7.48752 × 10−5 | 8.65473 × 10−19 | 1.63951 × 10−43 | 6.71261 × 10−58 |
| chr12:21178615 | SLCO1B1 | rs4149056 | 9.76195 × 10−23 | 7.37959 × 10−10 | 2.78132 × 10−8 | 2.34092 × 10−6 | 3.06483 × 10−17 |
| chr2:159230143 | TANC1 | rs4664277 | 1.61934 × 10−17 | 5.62787 × 10−8 | 3.42429 × 10−10 | 1.86855 × 10−24 | 5.11983 × 10−15 |
| chr12:47879112 | VDR | rs2228570 | 7.39552 × 10−41 | 6.3868 × 10−19 | 1.78553 × 10−20 | 4.04936 × 10−25 | 1.4743 × 10−38 |
| chr3:14145949 | XPC | rs2228001 | 2.43937 × 10−34 | 1.84077 × 10−32 | 4.51201 × 10−27 | 1.39797 × 10−23 | 1.02515 × 10−27 |
Supplementary Materials
The following supporting information can be downloaded at:
References
1. World Health OrganizationQuality Assurance and Safety of Medicines Team. Safety of Medicines: A Guide to Detecting and Reporting Adverse Drug Reactions: Why Health Professionals Need to Take Action 2002; WHO/EDM/QSM/2002.2 World Health Organization: Geneva, Switzerland, 2002.
2. Cacabelos, R.; Cacabelos, N.; Carril, J.C. The Role of Pharmacogenomics in Adverse Drug Reactions. Expert Rev. Clin. Pharmacol.; 2019; 12, pp. 407-442. [DOI: https://dx.doi.org/10.1080/17512433.2019.1597706] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30916581]
3. Kozyra, M.; Ingelman-Sundberg, M.; Lauschke, V.M. Rare Genetic Variants in Cellular Transporters, Metabolic Enzymes, and Nuclear Receptors Can Be Important Determinants of Interindividual Differences in Drug Response. Genet. Med.; 2017; 19, pp. 20-29. [DOI: https://dx.doi.org/10.1038/gim.2016.33] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27101133]
4. Zhou, Z.W.; Chen, X.W.; Sneed, K.B.; Yang, Y.X.; Zhang, X.; He, Z.X.; Chow, K.; Yang, T.; Duan, W.; Zhou, S.F. Clinical Association between Pharmacogenomics and Adverse Drug Reactions. Drugs; 2015; 75, pp. 589-631. [DOI: https://dx.doi.org/10.1007/s40265-015-0375-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25895462]
5. Malki, M.A.; Pearson, E.R. Drug-Drug-Gene Interactions and Adverse Drug Reactions. Pharm. J.; 2020; 20, pp. 355-366. [DOI: https://dx.doi.org/10.1038/s41397-019-0122-0]
6. Verstegen, R.H.J.; Ito, S. The Future of Precision Medicine. Clin. Pharmacol. Ther.; 2019; 106, pp. 903-906. [DOI: https://dx.doi.org/10.1002/cpt.1622]
7. Huddart, R.; Fohner, A.E.; Whirl-Carrillo, M.; Wojcik, G.L.; Gignoux, C.R.; Popejoy, A.B.; Bustamante, C.D.; Altman, R.B.; Klein, T.E. Standardized Biogeographic Grouping System for Annotating Populations in Pharmacogenetic Research. Clin. Pharmacol. Ther.; 2019; 105, pp. 1256-1262. [DOI: https://dx.doi.org/10.1002/cpt.1322]
8. Suarez-Kurtz, G.; de Araújo, G.S. Pharmacogenetic Differentiation across Latin America. Pharmacogenomics; 2022; 23, pp. 225-233. [DOI: https://dx.doi.org/10.2217/pgs-2021-0152]
9. Maroñas, O.; Latorre, A.; Dopazo, J.; Pirmohamed, M.; Rodríguez-Antona, C.; Siest, G.; Carracedo, Á.; Llerena, A. Progress in Pharmacogenetics: Consortiums and New Strategies. Drug Metab. Pers. Ther.; 2016; 31, pp. 17-23. [DOI: https://dx.doi.org/10.1515/dmpt-2015-0039]
10. Popejoy, A.B. Diversity in Precision Medicine and Pharmacogenetics: Methodological and Conceptual Considerations for Broadening Participation. Pharm. Pers. Med.; 2019; 12, pp. 257-271. [DOI: https://dx.doi.org/10.2147/PGPM.S179742]
11. Amorim, C.E.G.; Wang, S.; Marrero, A.R.; Salzano, F.M.; Ruiz-Linares, A.; Bortolini, M.C. X-Chromosomal Genetic Diversity and Linkage Disequilibrium Patterns in Amerindians and Non-Amerindian Populations. Am. J. Hum. Biol.; 2011; 23, pp. 299-304. [DOI: https://dx.doi.org/10.1002/ajhb.21110]
12. Rodrigues, J.C.G.; Fernandes, M.R.; Guerreiro, J.F.; da Costa da Silva, A.L.; Ribeiro-dos-Santos, Â.; Santos, S.; dos Santos, N.P.C. Polymorphisms of ADME-Related Genes and Their Implications for Drug Safety and Efficacy in Amazonian Amerindians. Sci. Rep.; 2019; 9, 7201. [DOI: https://dx.doi.org/10.1038/s41598-019-43610-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31076604]
13. Fernandes, M.R.; Rodrigues, J.C.G.; Maroñas, O.; Pellicer, A.L.; Cruz, R.; Guerreiro, J.F.; Burbano, R.M.R.; de Assumpção, P.P.; Ribeiro-Dos-santos, A.; dos Santos, S.E.B. et al. Genetic Diversity of Drug-Related Genes in Native Americans of the Brazilian Amazon. Pharm. Pers. Med.; 2021; 14, 117. [DOI: https://dx.doi.org/10.2147/PGPM.S274741] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33519226]
14. de Andrade Ramos, B.R.; Mendes, N.D.; Tanikawa, A.A.; Amador, M.A.T.; dos Santos, N.P.C.; dos Santos, S.E.B.; Castelli, E.C.; Witkin, S.S.; da Silva, M.G. Ancestry Informative Markers and Selected Single Nucleotide Polymorphisms in Immunoregulatory Genes on Preterm Labor and Preterm Premature Rupture of Membranes: A Case Control Study. BMC Pregnancy Childbirth; 2016; 16, 30. [DOI: https://dx.doi.org/10.1186/s12884-016-0823-1]
15. Ribeiro-dos-Santos, A.M.; Vidal, A.F.; Vinasco-Sandoval, T.; Guerreiro, J.; Santos, S.; Ribeiro-dos-Santos, Â.; de Souza, S.J. Exome Sequencing of Native Populations From the Amazon Reveals Patterns on the Peopling of South America. Front. Genet.; 2020; 11, 548507. [DOI: https://dx.doi.org/10.3389/fgene.2020.548507] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33193622]
16. Rodrigues, J.C.G.; de Souza, T.P.; Pastana, L.F.; dos Santos, A.M.R.; Fernandes, M.R.; Pinto, P.; Wanderley, A.V.; de Souza, S.J.; Kroll, J.E.; Pereira, A.L. et al. Identification of NUDT15 Gene Variants in Amazonian Amerindians and Admixed Individuals from Northern Brazil. PLoS ONE; 2020; 15, e0231651. [DOI: https://dx.doi.org/10.1371/journal.pone.0231651]
17. Whirl-Carrillo, M.; Huddart, R.; Gong, L.; Sangkuhl, K.; Thorn, C.F.; Whaley, R.; Klein, T.E. An Evidence-Based Framework for Evaluating Pharmacogenomics Knowledge for Personalized Medicine. Clin. Pharmacol. Ther.; 2021; 110, pp. 563-572. [DOI: https://dx.doi.org/10.1002/cpt.2350]
18. Relling, M.V.; Klein, T.E. CPIC: Clinical Pharmacogenetics Implementation Consortium of the Pharmacogenomics Research Network. Clin. Pharmacol. Ther.; 2011; 89, pp. 464-467. [DOI: https://dx.doi.org/10.1038/clpt.2010.279]
19. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R Package for Multivariate Analysis. J. Stat. Softw.; 2008; 25, pp. 1-18. [DOI: https://dx.doi.org/10.18637/jss.v025.i01]
20. Giorgi, F.M.; Ceraolo, C.; Mercatelli, D. The R Language: An Engine for Bioinformatics and Data Science. Life; 2022; 12, 648. [DOI: https://dx.doi.org/10.3390/life12050648]
21. Dubois, P.F.; Hinsen, K.; Hugunin, J. Numerical Python. Comput. Phys.; 1998; 10, 262. [DOI: https://dx.doi.org/10.1063/1.4822400]
22. Cingolani, P.; Platts, A.; Wang, L.L.; Coon, M.; Nguyen, T.; Wang, L.; Land, S.J.; Lu, X.; Ruden, D.M. A Program for Annotating and Predicting the Effects of Single Nucleotide Polymorphisms, SnpEff: SNPs in the Genome of Drosophila Melanogaster Strain W1118; Iso-2; Iso-3. Fly; 2012; 6, pp. 80-92. [DOI: https://dx.doi.org/10.4161/fly.19695] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22728672]
23. Pastana, L.F.; Silva, T.A.; Gellen, L.P.A.; Vieira, G.M.; de Assunção, L.A.; Leitão, L.P.C.; da Silva, N.M.; Coelho, R.d.C.C.; de Alcântara, A.L.; Vinagre, L.W.M.S. et al. The Genomic Profile Associated with Risk of Severe Forms of COVID-19 in Amazonian Native American Populations. J. Pers. Med.; 2022; 12, 554. [DOI: https://dx.doi.org/10.3390/jpm12040554] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35455670]
24. da Costa, G.E.; Fernandes, G.L.; Rodrigues, J.C.G.; da Leal, D.F.; Pastana, L.F.; Pereira, E.E.B.; Assumpção, P.P.; Burbano, R.M.R.; dos Santos, S.E.B.; Guerreiro, J.F. et al. Exome Evaluation of Autism-Associated Genes in Amazon American Populations. Genes; 2022; 13, 368. [DOI: https://dx.doi.org/10.3390/genes13020368] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35205412]
25. Dobbin, E.A.F.; Medeiros, J.A.G.; Costa, M.S.C.R.; Rodrigues, J.C.G.; Guerreiro, J.F.; Kroll, J.E.; de Souza, S.J.; de Assumpção, P.P.; Ribeiro-Dos-santos, Â.; dos Santos, S.E.B. et al. Identification of Variants (Rs11571707, Rs144848, and Rs11571769) in the BRCA2 Gene Associated with Hereditary Breast Cancer in Indigenous Populations of the Brazilian Amazon. Genes; 2021; 12, 142. [DOI: https://dx.doi.org/10.3390/genes12020142] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33499154]
26. de Carvalho, D.C.; Wanderley, A.V.; dos Santos, A.M.R.; Moreira, F.C.; de Sá, R.B.A.; Fernandes, M.R.; Modesto, A.A.C.; de Souza, T.P.; Cohen-Paes, A.; Leitão, L.P.C. et al. Characterization of Pharmacogenetic Markers Related to Acute Lymphoblastic Leukemia Toxicity in Amazonian Native Americans Population. Sci. Rep.; 2020; 10, 10292. [DOI: https://dx.doi.org/10.1038/s41598-020-67312-y]
27. Fernandes, V.C.; Pretti, M.A.M.; Tsuneto, L.T.; Petzl-Erler, M.L.; Suarez-Kurtz, G. Single Nucleotide Variants as Proxies for HLA-A*31:01 in Native American Populations. Front. Pharmacol.; 2022; 13, 849136. [DOI: https://dx.doi.org/10.3389/fphar.2022.849136]
28. Moreira, R.G.; Saraiva-Duarte, J.M.; Pereira, A.C.; Sosa-Macias, M.; Galaviz-Hernandez, C.; Santolalla, M.L.; Magalhães, W.C.S.; Zolini, C.; Leal, T.P.; Balázs, Z. et al. Population Genetics of PDE4B (Phosphodiesterase-4B) in Neglected Native Americans: Implications for Cancer Pharmacogenetics. Clin. Transl. Sci.; 2022; online version of record [DOI: https://dx.doi.org/10.1111/cts.13266]
29. Amstutz, U.; Henricks, L.M.; Offer, S.M.; Barbarino, J.; Schellens, J.H.M.; Swen, J.J.; Klein, T.E.; McLeod, H.L.; Caudle, K.E.; Diasio, R.B. et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for Dihydropyrimidine Dehydrogenase Genotype and Fluoropyrimidine Dosing: 2017 Update. Clin. Pharmacol. Ther.; 2018; 103, pp. 210-216. [DOI: https://dx.doi.org/10.1002/cpt.911]
30. Johnson, J.A.; Caudle, K.E.; Gong, L.; Whirl-Carrillo, M.; Stein, C.M.; Scott, S.A.; Lee, M.T.; Gage, B.F.; Kimmel, S.E.; Perera, M.A. et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for Pharmacogenetics-Guided Warfarin Dosing: 2017 Update. Clin. Pharmacol. Ther.; 2017; 102, pp. 397-404. [DOI: https://dx.doi.org/10.1002/cpt.668]
31. Relling, M.V.; Schwab, M.; Whirl-Carrillo, M.; Suarez-Kurtz, G.; Pui, C.H.; Stein, C.M.; Moyer, A.M.; Evans, W.E.; Klein, T.E.; Antillon-Klussmann, F.G. et al. Clinical Pharmacogenetics Implementation Consortium Guideline for Thiopurine Dosing Based on TPMT and NUDT15 Genotypes: 2018 Update. Clin. Pharmacol. Ther.; 2019; 105, pp. 1095-1105. [DOI: https://dx.doi.org/10.1002/cpt.1304]
32. Fang, M.Z.; Jackson, S.S.; O’Brien, T.R. IFNL4: Notable Variants and Associated Phenotypes. Gene; 2020; 730, 144289. [DOI: https://dx.doi.org/10.1016/j.gene.2019.144289] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31846709]
33. Muir, A.; Gong, L.; Johnson, S.; Lee, M.; Williams, M.; Klein, T.; Caudle, K.; Nelson, D. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guidelines for IFNL3 (IL28B) Genotype and PEG Interferon-α–Based Regimens. Clin. Pharmacol. Ther.; 2013; 95, pp. 141-146. [DOI: https://dx.doi.org/10.1038/clpt.2013.203] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24096968]
34. Frias, M.; Rivero-Juárez, A.; López-López, P.; Rivero, A. Pharmacogenetics and the Treatment of HIV-/HCV-Coinfected Patients. Pharmacogenomics; 2018; 19, pp. 979-995. [DOI: https://dx.doi.org/10.2217/pgs-2018-0046] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29992850]
35. Amstutz, U.; Froehlich, T.K.; Largiadèr, C.R. Dihydropyrimidine Dehydrogenase Gene as a Major Predictor of Severe 5-Fluorouracil Toxicity. Pharmacogenomics; 2011; 12, pp. 1321-1336. [DOI: https://dx.doi.org/10.2217/pgs.11.72] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21919607]
36. Boisdron-Celle, M.; Capitain, O.; Faroux, R.; Borg, C.; Metges, J.P.; Galais, M.P.; Kaassis, M.; Bennouna, J.; Bouhier-Leporrier, K.; Francois, E. et al. Prevention of 5-Fluorouracil-Induced Early Severe Toxicity by Pre-Therapeutic Dihydropyrimidine Dehydrogenase Deficiency Screening: Assessment of a Multiparametric Approach. Semin. Oncol.; 2017; 44, pp. 13-23. [DOI: https://dx.doi.org/10.1053/j.seminoncol.2017.02.008] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28395758]
37. Kassogue, Y.; Diakite, B.; Kassogue, O.; Konate, I.; Tamboura, K.; Diarra, Z.; Dehbi, H.; Nadifi, S.; Bougadari Traore, C.; Dao, S. et al. Genetic Polymorphism of Drug Metabolism Enzymes (GSTM1, GSTT1 and GSTP1) in the Healthy Malian Population. Mol. Biol. Rep.; 2020; 47, pp. 393-400. [DOI: https://dx.doi.org/10.1007/s11033-019-05143-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31650384]
38. Nakanishi, G.; Bertagnolli, L.S.; Pita-Oliveira, M.; Scudeler, M.M.; Torres-Loureiro, S.; Almeida-Dantas, T.; Alves, M.L.C.; Cirino, H.S.; Rodrigues-Soares, F. GSTM1 and GSTT1 Polymorphisms in Healthy Volunteers—A Worldwide Systematic Review. Drug Metab. Rev.; 2022; 54, pp. 37-45. [DOI: https://dx.doi.org/10.1080/03602532.2022.2036996]
39. Prysyazhnyuk, V.; Bobkovych, K.; Sydorchuk, L.; Dzuryak, V.; Sydorchuk, R.; Prysiazhniuk, I.; Buzdugan, I.; Prysyazhnyuk, P. Glutathione-S-Transferases Genes-Promising Predictors of Hepatic Dysfunction. World J. Hepatol.; 2021; 13, 620. [DOI: https://dx.doi.org/10.4254/wjh.v13.i6.620]
40. Hu, X.-Y.; Huang, X.-Y.; Ma, J.; Zuo, Y.; Luo, N.-B.; Lai, S.-L.; Su, D.-K. GSTT1 and GSTM1 Polymorphisms Predict Treatment Outcome for Breast Cancer: A Systematic Review and Meta-Analysis. Tumour Biol.; 2016; 37, pp. 151-162. [DOI: https://dx.doi.org/10.1007/s13277-015-4401-3]
41. Xiao, Q.; Deng, D.; Li, H.; Ye, F.; Huang, L.; Zhang, B.; Ye, B.; Mo, Z.; Yang, X.; Liu, Z. GSTT1 and GSTM1 Polymorphisms Predict Treatment Outcome for Acute Myeloid Leukemia: A Systematic Review and Meta-Analysis. Ann. Hematol.; 2014; 93, pp. 1381-1390. [DOI: https://dx.doi.org/10.1007/s00277-014-2050-z]
42. Liew, S.C.; Gupta, E. das Methylenetetrahydrofolate Reductase (MTHFR) C677T Polymorphism: Epidemiology, Metabolism and the Associated Diseases. Eur. J. Med. Genet.; 2015; 58, pp. 1-10. [DOI: https://dx.doi.org/10.1016/j.ejmg.2014.10.004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25449138]
43. Giletti, A.; Esperon, P. Genetic Markers in Methotrexate Treatments. Pharm. J.; 2018; 18, pp. 689-703. [DOI: https://dx.doi.org/10.1038/s41397-018-0047-z]
44. Xie, P.; Mo, J.L.; Liu, J.H.; Li, X.; Tan, L.M.; Zhang, W.; Zhou, H.H.; Liu, Z.Q. Pharmacogenomics of 5-Fluorouracil in Colorectal Cancer: Review and Update. Cell Oncol.; 2020; 43, pp. 989-1001. [DOI: https://dx.doi.org/10.1007/s13402-020-00529-1] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32474853]
45. Song, Z.; Hu, Y.; Liu, S.; Jiang, D.; Yi, Z.; Benjamin, M.M.; Zhao, R. The Role of Genetic Polymorphisms in High-Dose Methotrexate Toxicity and Response in Hematological Malignancies: A Systematic Review and Meta-Analysis. Front. Pharmacol.; 2021; 12, 757464. [DOI: https://dx.doi.org/10.3389/fphar.2021.757464] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34744734]
46. Campbell, J.M.; Bateman, E.; Peters, M.D.J.; Bowen, J.M.; Keefe, D.M.; Stephenson, M.D. Fluoropyrimidine and Platinum Toxicity Pharmacogenetics: An Umbrella Review of Systematic Reviews and Meta-Analyses. Pharmacogenomics; 2016; 17, pp. 435-451. [DOI: https://dx.doi.org/10.2217/pgs.15.180] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26894782]
47. Zhong, L.; He, X.; Zhang, Y.; Chuan, J.L.; Chen, M.; Zhu, S.M.; Peng, Q. Relevance of Methylenetetrahydrofolate Reductase Gene Variants C677T and A1298C with Response to Fluoropyrimidine-Based Chemotherapy in Colorectal Cancer: A Systematic Review and Meta-Analysis. Oncotarget; 2018; 9, pp. 31291-31301. [DOI: https://dx.doi.org/10.18632/oncotarget.24933] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30131855]
48. Fernandes, M.R.; Rodrigues, J.C.G.; Dobbin, E.A.F.; Pastana, L.F.; da Costa, D.F.; Barra, W.F.; Modesto, A.A.C.; de Assumpção, P.B.; da Costa Silva, A.L.; dos Santos, S.E.B. et al. Influence of FPGS, ABCC4, SLC29A1, and MTHFR Genes on the Pharmacogenomics of Fluoropyrimidines in Patients with Gastrointestinal Cancer from the Brazilian Amazon. Cancer Chemother. Pharmacol.; 2021; 88, pp. 837-844. [DOI: https://dx.doi.org/10.1007/s00280-021-04327-w] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34331561]
49. Matsusaka, S.; Lenz, H.J. Pharmacogenomics of Fluorouracil-Based Chemotherapy Toxicity. Expert Opin. Drug Metab. Toxicol.; 2015; 11, pp. 811-821. [DOI: https://dx.doi.org/10.1517/17425255.2015.1027684]
50. Peters, G.J.; Smid, K.; Meijer, E.; van Groeningen, C.J.; Leon, L.G. Role of Genomic Factors beyond Thymidylate Synthase in the Prediction of Response to 5-Fluorouracil. Nucleosides Nucleotides Nucleic Acids; 2016; 35, pp. 595-603. [DOI: https://dx.doi.org/10.1080/15257770.2016.1218020]
51. Luo, Y.; Yang, Z.; Chen, Y.; Lu, X.; Quan, Y. Genomic and Immunological Features of Microsatellite Instability in Colon Cancer. Gene; 2021; 781, 145534. [DOI: https://dx.doi.org/10.1016/j.gene.2021.145534]
52. Blondy, S.; David, V.; Verdier, M.; Mathonnet, M.; Perraud, A.; Christou, N. 5-Fluorouracil Resistance Mechanisms in Colorectal Cancer: From Classical Pathways to Promising Processes. Cancer Sci.; 2020; 111, 3142. [DOI: https://dx.doi.org/10.1111/cas.14532] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32536012]
53. Soo, R.A.; Syn, N.; Lee, S.C.; Wang, L.; Lim, X.Y.; Loh, M.; Tan, S.H.; Zee, Y.K.; Wong, A.L.A.; Chuah, B. et al. Pharmacogenetics-Guided Phase I Study of Capecitabine on an Intermittent Schedule in Patients with Advanced or Metastatic Solid Tumours. Sci. Rep.; 2016; 6, 27826. [DOI: https://dx.doi.org/10.1038/srep27826] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27296624]
54. Zhou, S.F. Polymorphism of Human Cytochrome P450 2D6 and Its Clinical Significance: Part II. Clin. Pharmacokinet.; 2009; 48, pp. 761-804. [DOI: https://dx.doi.org/10.2165/11318070-000000000-00000] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19902987]
55. Zanger, U.M.; Schwab, M. Cytochrome P450 Enzymes in Drug Metabolism: Regulation of Gene Expression, Enzyme Activities, and Impact of Genetic Variation. Pharmacol. Ther.; 2013; 138, pp. 103-141. [DOI: https://dx.doi.org/10.1016/j.pharmthera.2012.12.007] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23333322]
56. Crews, K.R.; Gaedigk, A.; Dunnenberger, H.M.; Leeder, J.S.; Klein, T.E.; Caudle, K.E.; Haidar, C.E.; Shen, D.D.; Callaghan, J.T.; Sadhasivam, S. et al. Clinical Pharmacogenetics Implementation Consortium Guidelines for Cytochrome P450 2D6 Genotype and Codeine Therapy: 2014 Update. Clin. Pharmacol. Ther.; 2014; 95, pp. 376-382. [DOI: https://dx.doi.org/10.1038/clpt.2013.254] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24458010]
57. Hicks, J.K.; Bishop, J.R.; Sangkuhl, K.; Muller, D.J.; Ji, Y.; Leckband, S.G.; Leeder, J.S.; Graham, R.L.; Chiulli, D.L.; LLerena, A. et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for CYP2D6 and CYP2C19 Genotypes and Dosing of Selective Serotonin Reuptake Inhibitors. Clin. Pharmacol. Ther.; 2015; 98, pp. 127-134. [DOI: https://dx.doi.org/10.1002/cpt.147]
58. Hicks, J.K.; Sangkuhl, K.; Swen, J.J.; Ellingrod, V.L.; Müller, D.J.; Shimoda, K.; Bishop, J.R.; Kharasch, E.D.; Skaar, T.C.; Gaedigk, A. et al. Clinical Pharmacogenetics Implementation Consortium Guideline (CPIC) for CYP2D6 and CYP2C19 Genotypes and Dosing of Tricyclic Antidepressants: 2016 Update. Clin. Pharmacol. Ther.; 2017; 102, pp. 37-44. [DOI: https://dx.doi.org/10.1002/cpt.597]
59. Bell, G.C.; Caudle, K.E.; Whirl-Carrillo, M.; Gordon, R.J.; Hikino, K.; Prows, C.A.; Gaedigk, A.; Agundez, J.A.G.; Sadhasivam, S.; Klein, T.E. et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for CYP2D6 Genotype and Use of Ondansetron and Tropisetron. Clin. Pharmacol. Ther.; 2017; 102, pp. 213-218. [DOI: https://dx.doi.org/10.1002/cpt.598]
60. Goetz, M.P.; Sangkuhl, K.; Guchelaar, H.-J.; Schwab, M.; Province, M.; Whirl-Carrillo, M.; Symmans, W.F.; Mcleod, H.L.; Ratain, M.J.; Zembutsu, H. et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) Guideline for CYP2D6 and Tamoxifen Therapy. Clin. Pharmacol. Ther.; 2014; 103, pp. 770-777. [DOI: https://dx.doi.org/10.1002/cpt.1007]
61. Brown, J.T.; Bishop, J.R.; Sangkuhl, K.; Nurmi, E.L.; Mueller, D.J.; Dinh, J.C.; Gaedigk, A.; Klein, T.E.; Caudle, K.E.; Mccracken, J.T. et al. Clinical Pharmacogenetics Implementation Consortium Guideline for Cytochrome P450 (CYP)2D6 Genotype and Atomoxetine Therapy. Clin. Pharmacol. Ther.; 2019; 106, pp. 94-102. [DOI: https://dx.doi.org/10.1002/cpt.1409]
62. Hoffecker, J.F.; Elias, S.A.; O’Rourke, D.H.; Scott, G.R.; Bigelow, N.H. Beringia and the Global Dispersal of Modern Humans. Evol. Anthropol.; 2016; 25, pp. 64-78. [DOI: https://dx.doi.org/10.1002/evan.21478] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27061035]
63. Reich, D.; Patterson, N.; Campbell, D.; Tandon, A.; Mazieres, S.; Ray, N.; Parra, M.V.; Rojas, W.; Duque, C.; Mesa, N. et al. Reconstructing Native American Population History. Nature; 2012; 488, pp. 370-374. [DOI: https://dx.doi.org/10.1038/nature11258] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22801491]
64. Raghavan, M.; Steinrücken, M.; Harris, K.; Schiffels, S.; Rasmussen, S.; DeGiorgio, M.; Albrechtsen, A.; Valdiosera, C.; Ávila-Arcos, M.C.; Malaspinas, A.S. et al. POPULATION GENETICS. Genomic Evidence for the Pleistocene and Recent Population History of Native Americans. Science; 2015; 349, aab3884. [DOI: https://dx.doi.org/10.1126/science.aab3884] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26198033]
65. Skoglund, P.; Mallick, S.; Bortolini, M.C.; Chennagiri, N.; Hünemeier, T.; Petzl-Erler, M.L.; Salzano, F.M.; Patterson, N.; Reich, D. Genetic Evidence for Two Founding Populations of the Americas. Nature; 2015; 525, pp. 104-108. [DOI: https://dx.doi.org/10.1038/nature14895] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26196601]
66. Moreno-Mayar, J.V.; Vinner, L.; de Barros Damgaard, P.; de La Fuente, C.; Chan, J.; Spence, J.P.; Allentoft, M.E.; Vimala, T.; Racimo, F.; Pinotti, T. et al. Early Human Dispersals within the Americas. Science; 2018; 362, eaav2621. [DOI: https://dx.doi.org/10.1126/science.aav2621]
67. Mendes, M.; Alvim, I.; Borda, V.; Tarazona-Santos, E. The History behind the Mosaic of the Americas. Curr. Opin. Genet. Dev.; 2020; 62, pp. 72-77. [DOI: https://dx.doi.org/10.1016/j.gde.2020.06.007]
68. Willerslev, E.; Meltzer, D.J. Peopling of the Americas as Inferred from Ancient Genomics. Nature; 2021; 594, pp. 356-364. [DOI: https://dx.doi.org/10.1038/s41586-021-03499-y]
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Marianne Rodrigues Fernandes 1
; Ribeiro-dos-Santos, André Maurício 2 ; Gilderlanio Santana de Araújo 2
; de Souza, Sandro José 3 ; João Farias Guerreiro 2
; Ribeiro-dos-Santos, Ândrea 4
; Paulo Pimentel de Assumpção 1 ; Ney Pereira Carneiro dos Santos 1 ; Santos, Sidney 4
1 Núcleo de Pesquisa em Oncologia, Universidade Federal do Pará, Belem 66073-000, Brazil;
2 Laboratório de Genética Humana e Médica, Universidade Federal do Pará, Belem 66075-110, Brazil;
3 BioME, Universidade Federal do Rio Grande do Norte, Natal 59078-400, Brazil;
4 Núcleo de Pesquisa em Oncologia, Universidade Federal do Pará, Belem 66073-000, Brazil;