Citations

Abstract Purpose: Cardiac toxicity is a serious potential complication of HER2-directed therapies and anthracyclines. HER2 codon 655 and SLC28A3 gene polymorphisms have been reported to be associated with cardiac toxicity from anti-HER2 and anthracycline therapy, respectively. Association of the polymorphism at HER2 codon 655 with prognosis has also been reported. Experimental Design: Whole blood samples from patients treated on a randomized adjuvant breast cancer trial (BCIRG-006) that compared chemotherapy with or without trastuzumab plus either anthracycline or nonanthracycline chemotherapy were tested for genetic polymorphisms in HER2 codon 655 and SLC28A3. Genotypes were correlated with cardiac function and disease-free survival (DFS) outcomes. Results: Of 3,222 patients enrolled in BCIRG-006, 662 patient samples were successfully genotyped for the rs1136201 allele in HER2 (codon 655): 424 (64%) were AA, 30 (4.5%) were GG, and 208 (31%) were AG genotype. In addition, 665 patient samples were successfully genotyped for the rs7853758 allele in the SLC28A3 gene: 19 (3%) were AA, 475 (71%) were GG, and 171 (26%) were AG genotype. Follow-up time was 10 years. No correlation between DFS, cardiac event rate, or mean left ventricular ejection fraction (LVEF) and rs1136201 genotype was seen in the trastuzumab-treated or non–trastuzumab-treated patients. Moreover, mean LVEF and cardiac event rates were similar in all rs7853758 genotype groups treated with anthracycline-based therapy. Conclusions: In the largest study to date to evaluate whether two polymorphisms are associated with DFS and/or cardiac toxicity in HER2-positive breast cancer treated with trastuzumab and/or anthracyclines, we observed no correlation.

Major depressive disorder (MDD), the number one psychiatric disease and a major cause of medical disability worldwide, is a heterogeneous disease. Efforts to identify biomarkers for antidepressant therapy response by the use of genome‐wide association studies (GWAS) have generally yielded disappointing results, e.g. only a few markers have been replicated across studies. By applying a “pharmacometabolomics‐informed pharmacogenomic” research strategy, we previously identified a genome‐wide significant single‐nucleotide polymorphism (SNP) signal in the glutamate‐rich 3 ( ERICH3 ) gene that was associated with plasma serotonin concentrations which were themselves associated with selective serotonin reuptake inhibitor (SSRI) response in MDD patients enrolled in the Mayo Clinic PGRN‐AMPS SSRI trial. The same ERICH3 SNPs were associated with SSRI response in other GWAS, including the STAR*D and ISPC studies. However, the function of ERICH3 and the molecular mechanism(s) by which it might be associated with plasma serotonin concentrations and with SSRI response remain unknown. Functional characterization of ERICH3 and its encoded protein might make it possible to identify novel mechanisms of antidepressant treatment response. In the present study, ERICH3 mRNA and protein were studied in cell lines from human peripheral blood and brain. ERICH3 cDNA constructs and specific antibodies were generated to make it possible to study its function. Two ERICH3 mRNA splice variants and protein isoforms, one “full‐length” and the other encoding a C‐terminal protein, were identified in a variety of human cell lines and in human brain tissue. ERICH3 mRNA expression and immunofluorescent staining suggested both cytosolic and vesicular localizations for ERICH3 as well as co‐localization with serotonin and dopamine in both mouse and human neurons. Co‐immunoprecipitation studies identified several ERICH3 interacting proteins, including PIK3C2A, CLTC and AP2A2, proteins known to play a role in vesicular trafficking, as well as ALK and CUX1, both of which display genetic polymorphisms that have been associated with SSRI response in GWAS studies of antidepressant treatment response in MDD patients. This series of experiments suggests that ERICH3 may play a role in vesicular trafficking and, as a result, may influence neurotransmitter actions in the brain that are associated with antidepressant treatment response. Support or Funding Information This work was supported by the NIH grants U19 GM61388, R01 GM28157 and T32 GM008685. This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal .

BACKGROUND Most functionally‐significant single nucleotide polymorphisms (SNPs) are outside of open reading frames. Mechanisms by which these SNPs contribute to gene regulation are poorly understood. Aryl hydrocarbon receptor (AHR) is a nuclear receptor that regulates drug metabolism, xenobiotic response and immune response. Our previous data indicate that SNPs outside response elements can impact transcription factor binding and gene regulation. We hypothesized that SNPs near but outside AHR response elements (AhREs) might impact AHR binding and gene regulation. METHODS We have collected dense SNP‐genotypes across 300 lymphoblastoid cell lines (LCLs) – a powerful tool for generating and testing pharmacogenomic hypotheses. Eight LCLs were used to generate hypotheses about AHR SNP‐dependent transcription utilizing RNA‐ and ChIP‐seq, which were tested in additional LCLs using qRT‐PCR and ChIP after AHR agonist or antagonist treatment. RESULTS 69 genes were highly significantly differentially expressed after agonist or antagonist treatment and 80% of those genes had an AhRE within 50kb. Five SNP‐gene pairs have been functionally validated in additional LCLs (Table 1). CONCLUSION SNPs that are outside of AhREs can alter AHR binding and gene transcription. This is a general phenomenon that may have large‐scale influence on gene transcriptional regulation and pharmacogenomics. This abstract is from the Experimental Biology 2018 Meeting. There is no full text article associated with this abstract published in The FASEB Journal .

Abstract Background: Androgen receptor (AR) is expressed in 60%-90% of breast cancers. The role of AR in breast cancer largely depends on estrogen receptor (ER) status and remains controversial. In ER+ cancers, AR expression is associated with improved prognosis. Enzalutamide (ENZ), an AR antagonist, impairs AR signaling, inhibits ER+/AR+ breast cancer cell proliferation and has been shown to mitigate resistance to anti-estrogen therapies. A recent study identified that RAD140, an AR agonist, suppressed the growth of ER+/AR+ breast cancer via stimulating AR signaling, resulting in down regulation of ERα. Both AR agonists and antagonists are effective in treating ER+ breast cancers. Identifying the subgroup of patients who most likely will benefit from an individual drug is important in clinical practice. Here, we sought to characterize the relative AR to ER levels and their relationships to drug response in ER+ breast cancers.Methods: We evaluated AR and ER expression levels among 68 samples with ER+ breast cancers. On each tumor, IHC staining for AR and ER were performed and results were scored by multiplying the percentage of positive cells by the intensity. TCGA-RPPA (reverse phase protein array) dataset was used to verify AR protein distribution in ER+ cancers. We created MCF7 and T47D cells with doxycycline-inducible AR and ER expression system to manipulate the relative AR to ER ratios and determined the antitumor activity of ENZ and RAD140 by cell proliferation assays. The levels of AR and ER in cells were measured by western blot. Linear regression was used to test the association between dose-response area under curve (AUC) and AR, ER levels. Drug preference was modeled against AR/ER expression levels (AR/ER ratios) using logistic regression.Results: Among 68 ER+ breast cancers, 69.12% were AR-positive. More than 2/3 cases (48 of 68) had more ER than AR expressed. The AR to ER ratios varied from 0 to 6. In the TCGA cohort, 84.15% of 347 ER+ patients had AR/ER ratios less than 1.The range of AR to ER ratios in TCGA dataset were comparable with our data (0.004 to 4.210). In our cell line models, the AR/ER ratios were controlled between 0.19 to 4.05. We found that the AUC of RAD140 was negatively associated with AR protein levels (P=0.008) and AR/ER ratios (P=0.013), and not significantly associated with ER expression levels. On the other hand, the AUC of ENZ was significantly negatively associated with ER protein (P=0.0016) but postiviely associated with AR/ER ratio (p=0.037). Preferred treatment comparing efficacy of the two drug can be best determined at extremes of AR/ER ratios. Using clinically relevant dosages, our model predicted that ENZ (10µM) would be a preferred treatment choice and have a better treatment efficacy compared with RAD140 when the AR/ER ratio was ≤ 0.42, whereas RAD140 (1µM) would be the preferred choice with a better treatment efficacy compared with ENZ when AR/ER ratios were ≥ 3.1 The efficacy preference of the two drugs are equivocal for AR/ER ratios between 0.42 and 3.10. Conclusion: We developed preclinical models using AR and ER expression levels to predict AR-targeting drug response. The results support the use of RAD140 in AR high patients and those with an AR/ER ratio >3.10, and enzalutamide in AR low patients and those with an AR/ER ratio <0.42. Equipoise on choice of drug was found for AR/ER ratios of 0.42-3.10. RAD140 and enzalutamide are compelling candidates for monotherapy or combination with anti-estrogen therapies in ER+/AR+ breast cancer. Future clinical validation of the models and therapeutic effect is warranted. Citation Format: Lixuan Wei, Huanyao Gao, Jia Yu, Duan Liu, Huan Zhang, Thanh Nguyen, Marie R Passow, Jodi M Carter, Richard M Weinshilboum, James N Ingle, Liewei Wang. Computational modeling of androgen receptor (AR) and estrogen receptor as predictive biomarkers of response to AR agonists and antagonists [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr PS5-36.

Several antiepileptic drugs (AEDs) have US Food and Drug Administration (FDA) approval for use as mood stabilizers in bipolar disorder (BD), but not all BD patients respond to these AED mood stabilizers (AED‐MSs). To identify genetic polymorphisms that contribute to the variability in AED‐MS response, we performed a discovery genome‐wide association study (GWAS) of 199 BD patients from the Mayo Clinic Bipolar Disorder Biobank. Most of these patients had been treated with the AED‐MS valproate/divalproex and/or lamotrigine. AED‐MS response was assessed using the Alda scale, which quantifies clinical improvement while accounting for potential confounding factors. We identified two genome‐wide significant single‐nucleotide polymorphism (SNP) signals that mapped to the THSD7A (rs78835388, P = 7.1E‐09) and SLC35F3 (rs114872993, P = 3.2E‐08) genes. We also identified two genes with statistically significant gene‐level associations: ABCC1 (P = 6.7E‐07; top SNP rs875740, P = 2.0E‐6), and DISP1 (P = 8.9E‐07; top SNP rs34701716, P = 8.9E‐07). THSD7A SNPs were previously found to be associated with risk for several psychiatric disorders, including BD. Both THSD7A and SLC35F3 are expressed in excitatory/glutamatergic and inhibitory/γ‐aminobutyric acidergic (GABAergic) neurons, which are targets of AED‐MSs. ABCC1 is involved in the transport of valproate and lamotrigine metabolites, and the SNPs in ABCC1 and DISP1 with the strongest evidence of association in our GWAS are strong splicing quantitative trait loci in the human gut, suggesting a possible influence on drug absorption. In conclusion, our pharmacogenomic study identified novel genetic loci that appear to contribute to AED‐MS treatment response, and may facilitate precision medicine in BD.

Single nucleotide variants in the open reading frames (ORFs) of pharmacogenes are important causes of interindividual variability in drug response. The functional characterization of variants of unknown significance within ORFs remains a major challenge for pharmacogenomics. Deep mutational scanning (DMS) is a high‐throughput technique that makes it possible to analyze the functional effect of hundreds of variants in a parallel and scalable fashion. We adapted a “landing pad” DMS system to study the function of missense variants in the ORFs of cytochrome P450 family 2 subfamily C member 9 ( CYP2C9 ) and cytochrome P450 family 2 subfamily C member 19 ( CYP2C19 ). We studied 230 observed missense variants in the CYP2C9 and CYP2C19 ORFs and found that 19 of 109 CYP2C9 and 36 of 121 CYP2C19 variants displayed less than ~ 25% of the wild‐type protein expression, a level that may have clinical relevance. Our results support DMS as an efficient method for the identification of damaging ORF variants that might have potential clinical pharmacogenomic application.

We previously reported that testis‐specific Y‐encoded‐like protein (TSPYLs) are transcription regulators for CYP3A4, CYP2C9, and CYP2C19. Here, we observed dual roles for TSPYLs in mediating serotonin transport and the metabolism of selective serotonin reuptake inhibitors (SSRIs) in patients with major depressive disorder (MDD). The widely prescribed SSRIs, citalopram, and escitalopram are metabolized mainly by CYP2C19. The TSPYL1 rs3828743 single nucleotide polymorphism (SNP), which decreases its suppression of CYP2C19 expression, was associated with rapid escitalopram metabolism and worse treatment response in the Mayo PGRN‐AMPS clinical trial. We also found that TSPYLs can regulate expression of the serotonin transporter protein, SLC6A4, and, in turn, serotonin transport into cells. The SNPs in tight linkage disequilibrium with the TSPYL1 rs10223646 SNP were significantly correlated with baseline severity of depression in patients with MDD in the Sequenced Treatment Alternatives to Relieve Depression and International SSRI Pharmacogenomics Consortium clinical trials. Our findings suggest that genetic variation in TSPYL genes may be novel indicators for baseline severity of depression and SSRI poor response.

Anastrozole is a widely prescribed aromatase inhibitor for the therapy of estrogen receptor positive (ER+) breast cancer. We performed a genome‐wide association study (GWAS) for plasma anastrozole concentrations in 687 postmenopausal women with ER+ breast cancer. The top single‐nucleotide polymorphism (SNP) signal mapped across SLC38A7 (rs11648166, P = 2.3E‐08), which we showed to encode an anastrozole influx transporter. The second most significant signal (rs28845026, P = 5.4E‐08) mapped near ALPPL2 and displayed epistasis with the SLC38A7 signal. Both of these SNPs were cis expression quantitative trait loci (eQTL)s for these genes, and patients homozygous for variant genotypes for both SNPs had the highest drug concentrations, the highest SLC38A7 expression, and the lowest ALPPL2 expression. In summary, our GWAS identified a novel gene encoding an anastrozole transporter, SLC38A7, as well as epistatic interaction between SNPs in that gene and SNPs near ALPPL2 that influenced both the expression of the transporter and anastrozole plasma concentrations.

The aryl hydrocarbon receptor (AHR) is a nuclear receptor that modulates the response to environmental stimuli. It was recognized historically for its role in toxicology but, in recent decades, it has been increasingly recognized as an important modulator of disease—especially for its role in modulating immune and inflammatory responses. AHR has been implicated in many diseases that are driven by immune/inflammatory processes, including major depressive disorder, multiple sclerosis, rheumatoid arthritis, asthma, and allergic responses, among others. The mechanisms by which AHR has been suggested to impact immune/inflammatory diseases include targeted gene expression and altered immune differentiation. It has been suggested that single nucleotide polymorphisms (SNPs) that are near AHR-regulated genes may contribute to AHR-dependent disease mechanisms/pathways. Further, we have found that SNPs that are outside of nuclear receptor binding sites (i.e., outside of AHR response elements (AHREs)) may contribute to AHR-dependent gene regulation in a SNP- and ligand-dependent manner. This review will discuss the evidence and mechanisms of AHR contributions to immune/inflammatory diseases and will consider the possibility that SNPs that are outside of AHR binding sites might contribute to AHR ligand-dependent inter-individual variation in disease pathophysiology and response to pharmacotherapeutics.

AbstractBackgroundRenal dysfunction occurs commonly after heart transplantation (HTx) with wide inter‐individual variability but whether a genetic predisposition exists in these patients is unknown. Genomewide association studies (GWAS) have not been performed to assess the association of genetic variation with change in renal function after HTx.MethodsClinical and demographic data of patients who underwent HTx and provided blood samples and consent for genetic analysis were included. Genotyping was performed using Illumina Infinium Human CoreExome v1.0 analysis kit. A GWAS utilizing linear regression models was performed with estimated glomerular filtration rate (eGFR) at 1 year as the phenotype after adjusting for baseline eGFR prior to HTx and conversion from calcineurin inhibitor to sirolimus as primary immunosuppression therapy.ResultsA total of 251 HTx recipients were genotyped for 314,903 single nucleotide polymorphisms (SNPs). The mean (SD) age was 50 (12.5) years; most patients were of European origin (n = 243, 96.8%) and males (n = 179, 71.3%). After adjustment for potential confounders, two variants, rs17033285 (P = 4.3 × 10−7) and rs4917601 (P = 6.46 × 10−7), in a long non‐coding RNA (lncRNA) gene LINC01121 and a pseudogene BTBD7P2, were identified to have a significant association with change in GFR at 1 year after HTx.ConclusionsOur first of its kind GWAS demonstrates that genetic variation affects renal function after HTx independent of other risk factors. Agnostic genetic approaches such as these may lead to identification of novel biological pathways such as the role of lncRNAs in the development of renal dysfunction post‐HTx.