Selective serotonin reuptake inhibitors (SSRIs) are standard of care for major depressive disorder (MDD) pharmacotherapy, but only approximately half of these patients remit on SSRI therapy. Our ...previous genome-wide association study identified a single-nucleotide polymorphism (SNP) signal across the glutamate-rich 3 (ERICH3) gene that was nearly genome-wide significantly associated with plasma serotonin (5-HT) concentrations, which were themselves associated with SSRI response for MDD patients enrolled in the Mayo Clinic PGRN-AMPS SSRI trial. In this study, we performed a meta-analysis which demonstrated that those SNPs were significantly associated with SSRI treatment outcomes in four independent MDD trials. However, the function of ERICH3 and molecular mechanism(s) by which it might be associated with plasma 5-HT concentrations and SSRI clinical response remained unclear. Therefore, we characterized the human ERICH3 gene functionally and identified ERICH3 mRNA transcripts and protein isoforms that are highly expressed in central nervous system cells. Coimmunoprecipitation identified a series of ERICH3 interacting proteins including clathrin heavy chain which are known to play a role in vesicular function. Immunofluorescence showed ERICH3 colocalization with 5-HT in vesicle-like structures, and ERICH3 knock-out dramatically decreased 5-HT staining in SK-N-SH cells as well as 5-HT concentrations in the culture media and cell lysates without changing the expression of 5-HT synthesizing or metabolizing enzymes. Finally, immunofluorescence also showed ERICH3 colocalization with dopamine in human iPSC-derived neurons. These results suggest that ERICH3 may play a significant role in vesicular function in serotonergic and other neuronal cell types, which might help explain its association with antidepressant treatment response.
We set out to determine whether machine learning–based algorithms that included functionally validated pharmacogenomic biomarkers joined with clinical measures could predict selective serotonin ...reuptake inhibitor (SSRI) remission/response in patients with major depressive disorder (MDD). We studied 1,030 white outpatients with MDD treated with citalopram/escitalopram in the Mayo Clinic Pharmacogenomics Research Network Antidepressant Medication Pharmacogenomic Study (PGRN‐AMPS; n = 398), Sequenced Treatment Alternatives to Relieve Depression (STAR*D; n = 467), and International SSRI Pharmacogenomics Consortium (ISPC; n = 165) trials. A genomewide association study for PGRN‐AMPS plasma metabolites associated with SSRI response (serotonin) and baseline MDD severity (kynurenine) identified single nucleotide polymorphisms (SNPs) in DEFB1, ERICH3, AHR, and TSPAN5 that we tested as predictors. Supervised machine‐learning methods trained using SNPs and total baseline depression scores predicted remission and response at 8 weeks with area under the receiver operating curve (AUC) > 0.7 (P < 0.04) in PGRN‐AMPS patients, with comparable prediction accuracies > 69% (P ≤ 0.07) in STAR*D and ISPC. These results demonstrate that machine learning can achieve accurate and, importantly, replicable prediction of SSRI therapy response using total baseline depression severity combined with pharmacogenomic biomarkers.
Tamoxifen is biotransformed to the potent anti-estrogen, endoxifen, by the cytochrome P450 (CYP) 2D6 enzyme. CYP2D6 genetic variation and inhibitors of the enzyme markedly reduce endoxifen plasma ...concentrations in tamoxifen-treated patients. Using a North Central Cancer Treatment Group adjuvant tamoxifen trial, we performed a comprehensive evaluation of CYP2D6 metabolism by assessing the combined effect of genetic variation and inhibition of the enzyme system on breast cancer recurrence and death.
Medical records were reviewed at each randomizing site to determine whether CYP2D6 inhibitors were co-prescribed with tamoxifen. Extensive metabolizers were defined as patients without a *4 allele (i.e., wt/wt) who were not co-prescribed a CYP2D6 inhibitor. Patients with decreased CYP2D6 metabolism were classified as intermediate or poor metabolizers (PM) based on the presence of one or two CYP2D6*4 alleles or the co-administration of a moderate or potent CYP2D6 inhibitor. The association between CYP2D6 metabolism and clinical outcome was assessed using Cox modeling.
Medication history was available in 225/256 eligible patients and CYP2D6*4 genotype in 190 patients. Thirteen patients (6%) were co-prescribed a CYP2D6 inhibitor potent (n = 3), moderate (n = 10), resulting in the following CYP2D6 metabolism: extensive (n = 115) and decreased (n = 65). In the multivariate analysis, patients with decreased metabolism had significantly shorter time to recurrence (p = 0.034; adj HR = 1.91; 95% CI 1.05-3.45) and worse relapse-free survival (RFS) (p = 0.017; adj HR = 1.74; 1.10-2.74); relative to patients with extensive metabolism. Cox' modeling demonstrated that compared to extensive metabolizers, PM had the most significant risk of breast cancer relapse (HR 3.12, p = 0.007).
CYP2D6 metabolism, as measured by genetic variation and enzyme inhibition, is an independent predictor of breast cancer outcome in post-menopausal women receiving tamoxifen for early breast cancer. Determination of CYP2D6 genotype may be of value in selecting adjuvant hormonal therapy and it appears CYP2D6 inhibitors should be avoided in tamoxifen-treated women.
Acamprosate is an anti-craving drug used in alcohol use disorder (AUD) pharmacotherapy. However, only a subset of patients achieves optimal treatment outcomes. The identification of predictive ...biomarkers of acamprosate treatment response in patients with AUD would be a substantial advance in addiction medicine. We designed this study to use proteomics data as a quantitative biological trait as a step toward identifying inflammatory modulators that might be associated with acamprosate treatment outcomes. The NIAAA-funded Mayo Clinic Center for the Individualized Treatment of Alcoholism study had previously recruited 442 AUD patients who received 3 months of acamprosate treatment. However, only 267 subjects returned for the 3-month follow-up visit and, as a result, had treatment outcome information available. Baseline alcohol craving intensity was the most significant predictor of acamprosate treatment outcomes. We performed plasma proteomics using the Olink target 96 inflammation panel and identified that baseline plasma TNF superfamily member 10 (TNFSF10) concentration was associated with alcohol craving intensity and variation in acamprosate treatment outcomes among AUD patients. We also performed RNA sequencing using baseline peripheral blood mononuclear cells from AUD patients with known acamprosate treatment outcomes which revealed that inflammation-related pathways were highly associated with relapse to alcohol use during the 3 months of acamprosate treatment. These observations represent an important step toward advancing our understanding of the pathophysiology of AUD and molecular mechanisms associated with acamprosate treatment response. In conclusion, applying omics-based approaches may be a practical approach for identifying biologic markers that could potentially predict alcohol craving intensity and acamprosate treatment response.
We previously discovered that the single nucleotide polymorphisms (SNP) rs9940645 in the ZNF423 gene regulate ZNF423 expression and serve as a potential biomarker for response to selective estrogen ...receptor modulators (SERMs). Here we explored pathways involved in ZNF423-mediated SERMs response and drugs that potentially sensitize SERMs.
RNA sequencing and label-free quantitative proteomics were performed to identify genes and pathways that are regulated by ZNF423 and the ZNF423 SNP. Both cultured cells and mouse xenograft models with different ZNF423 SNP genotypes were used to study the cellular responses to metformin.
We identified ribosome and AMP-activated protein kinase (AMPK) signaling as potential pathways regulated by ZNF423 or ZNF423 rs9940645 SNP. Moreover, using clustered regularly interspaced short palindromic repeats/Cas9-engineered ZR75-1 breast cancer cells with different ZNF423 SNP genotypes, striking differences in cellular responses to metformin, either alone or in the combination of tamoxifen, were observed in both cell culture and the mouse xenograft model.
We found that AMPK signaling is modulated by the ZNF423 rs9940645 SNP in estrogen and SERM-dependent fashion. The ZNF423 rs9940645 SNP affects metformin response in breast cancer and could be a potential biomarker for tailoring the metformin treatment.
Aims
Citalopram (CT) and escitalopram (S‐CT) are among the most widely prescribed selective serotonin reuptake inhibitors used to treat major depressive disorder (MDD). We applied a genome‐wide ...association study to identify genetic factors that contribute to variation in plasma concentrations of CT or S‐CT and their metabolites in MDD patients treated with CT or S‐CT.
Methods
Our genome‐wide association study was performed using samples from 435 MDD patients. Linear mixed models were used to account for within‐subject correlations of longitudinal measures of plasma drug/metabolite concentrations (4 and 8 weeks after the initiation of drug therapy), and single‐nucleotide polymorphisms (SNPs) were modelled as additive allelic effects.
Results
Genome‐wide significant associations were observed for S‐CT concentration with SNPs in or near the CYP2C19 gene on chromosome 10 (rs1074145, P = 4.1 × 10−9) and with S‐didesmethylcitalopram concentration for SNPs near the CYP2D6 locus on chromosome 22 (rs1065852, P = 2.0 × 10−16), supporting the important role of these cytochrome P450 (CYP) enzymes in biotransformation of citalopram. After adjustment for the effect of CYP2C19 functional alleles, the analyses also identified novel loci that will require future replication and functional validation.
Conclusions
In vitro and in vivo studies have suggested that the biotransformation of CT to monodesmethylcitalopram and didesmethylcitalopram is mediated by CYP isozymes. The results of our genome‐wide association study performed in MDD patients treated with CT or S‐CT have confirmed those observations but also identified novel genomic loci that might play a role in variation in plasma levels of CT or its metabolites during the treatment of MDD patients with these selective serotonin reuptake inhibitors.
Metabolomics provides valuable tools for the study of drug effects, unraveling the mechanism of action and variation in response due to treatment. In this study we used electrochemistry-based ...targeted metabolomics to gain insights into the mechanisms of action of escitalopram/citalopram focusing on a set of 31 metabolites from neurotransmitter-related pathways. Overall, 290 unipolar patients with major depressive disorder were profiled at baseline, after 4 and 8 weeks of drug treatment. The 17-item Hamilton Depression Rating Scale (HRSD
) scores gauged depressive symptom severity. More significant metabolic changes were found after 8 weeks than 4 weeks post baseline. Within the tryptophan pathway, we noted significant reductions in serotonin (5HT) and increases in indoles that are known to be influenced by human gut microbial cometabolism. 5HT, 5-hydroxyindoleacetate (5HIAA), and the ratio of 5HIAA/5HT showed significant correlations to temporal changes in HRSD
scores. In the tyrosine pathway, changes were observed in the end products of the catecholamines, 3-methoxy-4-hydroxyphenylethyleneglycol and vinylmandelic acid. Furthermore, two phenolic acids, 4-hydroxyphenylacetic acid and 4-hydroxybenzoic acid, produced through noncanconical pathways, were increased with drug exposure. In the purine pathway, significant reductions in hypoxanthine and xanthine levels were observed. Examination of metabolite interactions through differential partial correlation networks revealed changes in guanosine-homogentisic acid and methionine-tyrosine interactions associated with HRSD
. Genetic association studies using the ratios of these interacting pairs of metabolites highlighted two genetic loci harboring genes previously linked to depression, neurotransmission, or neurodegeneration. Overall, exposure to escitalopram/citalopram results in shifts in metabolism through noncanonical pathways, which suggest possible roles for the gut microbiome, oxidative stress, and inflammation-related mechanisms.
We performed a case-control genome-wide association study (GWAS) to identify single nucleotide polymorphisms (SNPs) associated with musculoskeletal adverse events (MS-AEs) in women treated with ...aromatase inhibitors (AIs) for early breast cancer.
A nested case-control design was used to select patients enrolled onto the MA.27 phase III trial comparing anastrozole with exemestane. Cases were matched to two controls and were defined as patients with grade 3 or 4 MS-AEs (according to the National Cancer Institute's Common Terminology Criteria for Adverse Events v3.0) or those who discontinued treatment for any grade of MS-AE within the first 2 years. Genotyping was performed with the Illumina Human610-Quad BeadChip.
The GWAS included 293 cases and 585 controls. A total of 551,358 SNPs were analyzed, followed by imputation and fine mapping of a region of interest on chromosome 14. Four SNPs on chromosome 14 had the lowest P values (2.23E-06 to 6.67E-07). T-cell leukemia 1A (TCL1A) was the gene closest (926-7000 bp) to the four SNPs. Functional genomic studies revealed that one of these SNPs (rs11849538) created an estrogen response element and that TCL1A expression was estrogen dependent, was associated with the variant SNP genotypes in estradiol-treated lymphoblastoid cells transfected with estrogen receptor alpha and was directly related to interleukin 17 receptor A (IL17RA) expression.
This GWAS identified SNPs associated with MS-AEs in women treated with AIs and with a gene (TCL1A) which, in turn, was related to a cytokine (IL17). These findings provide a focus for further research to identify patients at risk for MS-AEs and to explore the mechanisms for these adverse events.
Despite potential clinical benefits, implementation of pharmacogenomics (PGx) faces many technical and clinical challenges. These challenges can be overcome with a comprehensive and systematic ...implementation model.
The development and implementation of PGx were organized into eight interdependent components addressing resources, governance, clinical practice, education, testing, knowledge translation, clinical decision support (CDS), and maintenance. Several aspects of implementation were assessed, including adherence to the model, production of PGx-CDS interventions, and access to educational resources.
Between August 2012 and June 2015, 21 specific drug–gene interactions were reviewed and 18 of them were implemented in the electronic medical record as PGx-CDS interventions. There was complete adherence to the model with variable production time (98–392 days) and delay time (0–148 days). The implementation impacted approximately 1,247 unique providers and 3,788 unique patients. A total of 11 educational resources complementary to the drug–gene interactions and 5 modules specific for pharmacists were developed and implemented.
A comprehensive operational model can support PGx implementation in routine prescribing. Institutions can use this model as a roadmap to support similar efforts. However, we also identified challenges that will require major multidisciplinary and multi-institutional efforts to make PGx a universal reality.
Genet Med19 4, 421–429.