Search

Supplemental Digital Content is available in the text., Objectives 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. Methods 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. Results 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. Conclusions 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.

SLCO1B1 (solute carrier organic anion transporter family member 1B1) is an important transmembrane hepatic uptake transporter. Genetic variants in the SLCO1B1 gene have been associated with altered protein folding, resulting in protein degradation and decreased transporter activity. Next-generation sequencing (NGS) of pharmacogenes is being applied increasingly to associate variation in drug response with genetic sequence variants. However, it is difficult to link variants of unknown significance with functional phenotypes using “one-at-a-time” functional systems. Deep mutational scanning (DMS) using a “landing pad cell–based system” is a high-throughput technique designed to analyze hundreds of gene open reading frame (ORF) missense variants in a parallel and scalable fashion. We have applied DMS to analyze 137 missense variants in the SLCO1B1 ORF obtained from the Exome Aggregation Consortium project. ORFs containing these variants were fused to green fluorescent protein and were integrated into “landing pad” cells. Florescence-activated cell sorting was performed to separate the cells into four groups based on fluorescence readout indicating protein expression at the single cell level. NGS was then performed and SLCO1B1 variant frequencies were used to determine protein abundance. We found that six variants not previously characterized functionally displayed less than 25% and another 12 displayed approximately 50% of wild-type protein expression. These results were then functionally validated by transporter studies. Severely damaging variants identified by DMS may have clinical relevance for SLCO1B1-dependent drug transport, but we need to exercise caution since the relatively small number of severely damaging variants identified raise questions with regard to the application of DMS to intrinsic membrane proteins such as organic anion transporter protein 1B1. Significance Statement The functional implications of a large numbers of open reading frame (ORF) “variants of unknown significance” (VUS) in transporter genes have not been characterized. This study applied deep mutational scanning to determine the functional effects of VUS that have been observed in the ORF of SLCO1B1(solute carrier organic anion transporter family member 1B1). Several severely damaging variants were identified, studied, and validated. These observations have implications for both the application of deep mutational scanning to intrinsic membrane proteins and for the clinical effect of drugs and endogenous compounds transported by SLCO1B1.

Background: Breast cancer is the most frequently diagnosed cancer and the leading cause of cancer-related death among women worldwide. Triple negative breast cancer (TNBC) represents 15-20% of all breast cancers and is an aggressive subtype. Chemotherapy-based treatments remain the standard of care for TNBC. Unfortunately, chemotherapy resistance is common, and for these patients, outcomes are poor and alternative treatment strategies remain an unmet need. circRNAs are a newly identified class of noncoding RNA molecules with covalently closed circular structures. An increasing number of recent studies including ours have indicated that circRNAs play crucial roles in regulating tumor development and chemoresistance. However, the role of circRNAs in the process of chemotherapy resistance and TNBC progression is not clear. Materials and Methods: As a first step towards identifying circRNAs that participate in the development of chemoresistance in TNBC cells, and to determine if targeting such circRNAs is a novel and efficacious therapeutic strategy, doxorubicin-resistant (Doxo-R), paclitaxel-resistant (PTX-R), and double-resistant (DP-R) cell lines were generated from MDA-MB-231. Human circRNA microarrays were utilized to profile the expression of approximately 14,000 known circRNAs in normal breast tissue, matched patient-derived xenografts (PDX) generated prior to and following neoadjuvant chemotherapy (NAC), and TNBC chemosensitive and chemoresistant cell lines. Top hits were validated using RT-PCR. Results: circRNA microarray profiling identified 429 and 310 transcripts differentially expressed in doxorubicin and paclitaxel resistant cells, respectively, compared to parental chemosensitive cell lines (|FC| ≥ 1.5; p value < 0.05). In comparison to pre-NAC derived xenografts, 1,396 circRNAs were dysregulated among post-NAC PDX models. Further, three circRNAs (hsa_circ_001388, hsa_circ_104652, and hsa_circ_061260) were upregulated in Doxo-R, PTX-R, and post-NAC PDX samples compared to their respective controls. Among these three circRNAs, hsa_circ_001388 (also known as circNSD2) was the only transcript also predicted to be translated into a novel and uncharacterized protein given the presence of a high confidence translation initiation site and IRES sequence. Ongoing studies are aimed at determining the role of circNSD2 protein in breast cancer carcinogenesis, progression, and response to standard of care chemotherapeutics and the mechanistic process by which this protein functions. Additionally, the efficacy of specifically targeting this circRNA as a novel therapeutic approach is being explored. Conclusions: Increasing knowledge of the important functions of circRNAs underlying drug resistance will provide new opportunities for developing efficacious therapeutic strategies and prognostic/predictive biomarkers for TNBC. Citation Format: Xiyin Wang, Michael J. Emch, Xiaojia Tang, Jia Yu, Krishna R. Kalari, Liewei Wang, Matthew P. Goetz, John R. Hawse. The role of circular RNAs in triple negative breast cancer and chemotherapy resistance [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 1689.

Background: Cancer is usually treated with a combination of drugs rather than a single agent. Treating cancer with a drug combination decreases the likelihood of the cancer acquiring resistance to therapy. It also allows for lower dosages of individual drugs, which reduces the occurrence and severity of side effects. Using public datasets, we investigated potential anti-cancer drug combinations using a novel framework, ReCorDE (Recurrent Correlation of Drugs with Enrichment), which identifies correlated drug response patterns for drugs with different primary mechanisms of action. Methods: 250 drugs from CTRPv2, GDSCv2, and PRISM datasets were examined using normalized logIC50 or AUC measurements. ReCorDE consists of a correlation and an enrichment step. For each dataset, we constructed pairwise drug-drug relationships using Spearman’s correlation. Combinations significant (p < 0.05) in less than two datasets after Benjamini-Hochberg (BH) correction were pruned. To turn each drug combination into a drug class combination, all drugs were mapped to their corresponding fourth-level Anatomical Therapeutic Chemical (ATC) codes. Then, using a hypergeometric test, for each class combination in the subset of significant drug combinations (the “inside” set), we tested for enrichment of that class combination in the inside set compared to the class combination universe, the distribution of class combinations for all unique pairwise drug-drug associations (250 unique drugs; 31,125 distinct drug combinations mapped to 1596 distinct class combinations). Multiple testing during enrichment was corrected using the BH method with a significance cutoff of p< 0.05. Results: We found 2764 pairs to be significant at a = 0.05 using the AUC measure. With a correlation coefficient of 0.77, vincristine and YK-4-279, an EWS-FLI1/RNA Helicase A inhibitor, had the lowest computed p-value of of 1.39 × 10−120 Adjusted p-values for >15% of these combinations had p-values < 5 × 10−6. Class combination enrichment on the same set showed 132 class pairs were enriched in the inside set compared to the class combination universe after BH correction. Taxanes/Plk1 inhibitors (OR = 69.5, adjusted p = 5.6 × 10−11); Aurora Kinase inhibitors/histone modifying agents (OR=9.2, adjusted p = 6.13 × 10−11); and pyrimidine analogues/CDK inhibitors (OR=3.43, adjusted p = 0.02) are a few notable enriched class combinations. The enrichment of these class combinations in the inside set suggests that combinations of drugs from these classes may be particularly effective in treating cancer compared to other drug combinations. Conclusions: Our framework, ReCorDE, demonstrates that finding potential drug combinations and characterizing novel, frequently perturbed pathways outside of a drug's primary mechanisms of action can be accomplished by identifying correlated drug-drug pairs from large, publicly accessible databases. Citation Format: Emily T. Ghose, August John, Huanyao Gao, Krishna R. Kalari, Liewei Wang. ReCorDE: A novel computational framework to discover potential combinations of anti-cancer drugs. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4290.

Background: Cyclin dependent 4/6 kinase inhibitors (CDK4/6i) and endocrine therapy (ET) have improved progression-free survival (PFS) and overall survival in hormone-receptor (HR)-positive metastatic breast cancer (MBC), but progression of disease is inevitable. Serum thymidine kinase-1 (TK1) is a secreted marker of proliferation that is prognostic in patients (pts) with HR-positive, HER2-negative MBC and may be predictive of ET and CDK 4/6i response. PROMISE (NCT0281902) is a prospective study enrolling women with HR-positive MBC starting palbociclib (Pb) + letrozole (L) (1st line) or Pb + fulvestrant (2nd line). We undertook a comprehensive “omic” assessment of blood, tumor, urine and the fecal microbiome in order to identify novel genomic variants and pathways associated with an early decline in TK1 (measured after 2 months) and PFS. Additionally, patient derived xenografts/organoids were generated at baseline and progression to test new therapeutic approaches to overcome resistance to CDK4/6i and ET. We report the initial association between the baseline genomic landscape and baseline TK1 levels. Methods: FFPE tumor biopsies were obtained for DNA/RNA sequencing (TempusTM) and blood samples for TK1 (Divitum® assay, Biovica) were collected pretreatment (pre-Pb) and after 2 cycles of Pb + ET (post-Pb2). Both whole-exome (exome capture) sequencing (WES) and RNA-Seq used the Integrated DNA Technologies xGen Exome Research Panel v1.0 capture kit. TK1+ disease was defined as > 200 Du/L and TK1- disease as below limit of detection up to 200 Du/L. We tested the association between genomic and transcriptomic characteristics with baseline TK1 data in pretreatment samples where both WES and RNA-seq and TK1 was available. The data were analyzed using bioinformatics pipelines for somatic and germline mutations and copy number alterations. The current analysis focuses on baseline 1st-line pre-Pb omics data in conjunction with baseline TK1 levels. The database was locked for analysis on 5/29/2020. Results: Thirty-three pts (median age: 59 yrs.) were evaluable, with paired samples for TK1 in 32. Six pts had TK1+ disease pre-Pb and post-Pb2. Twenty-two pts had TK1- disease pre-Pb and post-Pb2. Four pts had a decrease in TK1 after 2 cycles of treatment that altered the classification from TK1+ to TK1-. Both baseline RNA seq and serum TK1 (n=16) were available for 4 TK1+ and 12 TK1- pts. In this group, 476 genes were differentially regulated (398 upregulated; 78 downregulated). Pathway analysis demonstrated enrichment in complement and coagulation cascade pathway, PPAR signaling pathway, and metabolism-related pathways related to up-regulation of CYP and UGT gene families. Further testing for the association of WES data with baseline TK1+ (n=8) and TK1- (n=16) disease demonstrated somatic copy number variations on chromosomes 6, 11, 12 and 15. CDK4 copy number gains were observed in 3/8 TK1+ pts and 0/16 TK1- pts. We also observed that somatic mutations (LOH, copy number and/or SNV/INDELs) were more prevalent in the TK1+ compared to the TK1- pre-Pb group in several cancer-associated genes (FAS [p=0.06] PTEN, PIK3CB, NAB2, SOX9 and FAT1 [p=0.08], TP53, and MAP2K4 [p=0.22]). Conversely, we also noted that 6/7 pts with GATA3 mutations had TK1- disease (p=0.23). Conclusions: Using a comprehensive “omics” approach, our data suggest that a secreted biomarker of proliferation (TK1) obtained prior to initiating CDK4/6i and ET for the first line treatment of HR+ MBC is associated with established and novel genes and pathways associated with prognosis of pts receiving ET and CDK 4/6i. Analysis of on-treatment (after 2 cycles) tumor RNA seq and its association with change in TK1 as well as data from the 2nd-line cohort will be presented at the meeting. Citation Format: Ciara C O Sullivan, Krishna R Kalari, Vera J Suman, Peter T Vedell, Ann Moyer, Abraham D Eyman Casey, Jason Sinnwell, Xiaojia Tang, Kevin Thompson, Alvaro Moreno-Aspitia, Donald W Northfelt, Minetta C Liu, Tufia C Haddad, Saranya Chumsri, Prema Peethambaram, Kathryn J Ruddy, Karthik V Giridhar, Roberto A Leon-Ferre, Adrian Nordmark, Mattias Bergqvist, Brendan P McMenomy, Richard M Weinshilboum, Liewei Wang, Matthew P Goetz. Novel genomic variants and pathways associated with baseline serum thymidine kinase 1 levels in HR-positive HER2-negative MBC patients commencing palbociclib and letrozole [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-24.

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 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.

Chromosome region maintenance protein 1 (CRM1) mediates protein export from the nucleus and is a new target for anticancer therapeutics. Broader application of KPT-330 (selinexor), a first-in-class CRM1 inhibitor recently approved for relapsed multiple myeloma and diffuse large B-cell lymphoma, have been limited by substantial toxicity. We discovered that salicylates markedly enhance the antitumor activity of CRM1 inhibitors by extending the mechanisms of action beyond CRM1 inhibition. Using salicylates in combination enables targeting of a range of blood cancers with a much lower dose of selinexor, thereby potentially mitigating prohibitive clinical adverse effects. Choline salicylate (CS) with low-dose KPT-330 (K+CS) had potent, broad activity across high-risk hematological malignancies and solid-organ cancers ex vivo and in vivo. The K+CS combination was not toxic to nonmalignant cells as compared with malignant cells and was safe without inducing toxicity to normal organs in mice. Mechanistically, compared with KPT-330 alone, K+CS suppresses the expression of CRM1, Rad51, and thymidylate synthase proteins, leading to more efficient inhibition of CRM1-mediated nuclear export, impairment of DNA-damage repair, reduced pyrimidine synthesis, cell-cycle arrest in S-phase, and cell apoptosis. Moreover, the addition of poly (ADP-ribose) polymerase inhibitors further potentiates the K+CS antitumor effect. K+CS represents a new class of therapy for multiple types of blood cancers and will stimulate future investigations to exploit DNA-damage repair and nucleocytoplasmic transport for cancer therapy in general.

Heterogeneity in the clinical presentation of major depressive disorder and response to antidepressants limits clinicians’ ability to accurately predict a specific patient’s eventual response to therapy. Validated depressive symptom profiles may be an important tool for identifying poor outcomes early in the course of treatment. To derive these symptom profiles, we first examined data from 947 depressed subjects treated with selective serotonin reuptake inhibitors (SSRIs) to delineate the heterogeneity of antidepressant response using probabilistic graphical models (PGMs). We then used unsupervised machine learning to identify specific depressive symptoms and thresholds of improvement that were predictive of antidepressant response by 4 weeks for a patient to achieve remission, response, or nonresponse by 8 weeks. Four depressive symptoms (depressed mood, guilt feelings and delusion, work and activities and psychic anxiety) and specific thresholds of change in each at 4 weeks predicted eventual outcome at 8 weeks to SSRI therapy with an average accuracy of 77% (p = 5.5E-08). The same four symptoms and prognostic thresholds derived from patients treated with SSRIs correctly predicted outcomes in 72% (p = 1.25E-05) of 1996 patients treated with other antidepressants in both inpatient and outpatient settings in independent publicly-available datasets. These predictive accuracies were higher than the accuracy of 53% for predicting SSRI response achieved using approaches that (i) incorporated only baseline clinical and sociodemographic factors, or (ii) used 4-week nonresponse status to predict likely outcomes at 8 weeks. The present findings suggest that PGMs providing interpretable predictions have the potential to enhance clinical treatment of depression and reduce the time burden associated with trials of ineffective antidepressants. Prospective trials examining this approach are forthcoming.

Objective To assess the potential impact of Pharmacogenomic (PGx) variation in cytochrome P450 2D6 (CYP2D6) enzyme function, using loss in quality-adjusted life years (QALYs) associated with treatment problems, and the willingness to pay to avoid treatment problems from patients’ and payers’ perspectives. Patients and Methods The study included patients prescribed tramadol or codeine, or both, between January 1, 2005, and December 31, 2017. Demographic information and adverse drug events, including adverse drug events and poor pain control, were collected from the electronic health records using natural language processing techniques and review by trained abstractors. Patients’ willingness to pay and QALY estimates were based on comprehensive literature review. The CYP2D6 phenotypes were divided into 4 groups: ultra-rapid metabolizers, normal metabolizers, intermediate metabolizers, and poor metabolizers. Results Among the 2860 identified patients, 63 (2%) were ultrarapid metabolizers, 1449 (50%) were normal metabolizers, 1155 (40%) were intermediate metabolizers, and 193 (7%) were poor metabolizers. The patients’ average estimated willingness-to-pay value to avoid treatment problems was $23 per month; poor metabolizers developed problems with the highest estimated willingness-to-pay value ($32 per month). The mean QALY loss among all patients was 0.024 QALYs (8.8 healthy days); poor metabolizers had the highest loss (0.027 QALYs, 9.9 healthy days). Conclusion Patients with various phenotypes developed different treatment problem profiles. Poor CYP2D6 metabolizers developed problems with highest willingness to pay, and they might potentially benefit most from PGx-guided treatment and problem prevention.

Triple negative breast cancer (TNBC) patients who fail to achieve a pathological complete response to neoadjuvant chemotherapy (NAC) will likely experience recurrence of the disease within 3-4 years. Prognostic assessment of early recurrence could facilitate clinical decisions and impact disease survival. This study investigated pre- and post-NAC multi-omics data from both an in-house and a public clinical study to identify biomarkers of NAC response. We observed significant transcriptional differences associated with response in the residual disease (post-NAC biopsies), which did not exist in the pre-NAC biopsies. We further refined the post-NAC transcriptional changes to a 17-gene signature and machine learning models were applied to evaluate the signature’s diagnostic potential (area under the curve ≥ 0.8). Interestingly, the signature was enriched with down-regulated immune genes and several of these genes demonstrated prognostic potential in basal TNBC tumors. Our 17-gene signature provides additional insight into TNBC recurrence, which warrants further investigation.

Metastatic castration-resistant prostate cancer (mCRPC) presents very low survival rates due to lack of response or acquired resistance to the available therapies. To date no molecular mechanisms of resistance have been identified, pointing out their complex dynamics. To identify key genes and processes associated with phenotypically-driven regulatory differences, we developed TraRe, a computational method that provides a three-tier analysis: i) at the network level, inferring differentially regulated modules; ii) at the regulon level, identifying regulatory relationships linked to phenotypic differences; and iii) at the single gene level, identifying TFs consistently linked to rewired modules. We applied TraRe (available in Bioconductor with full documentation) to transcriptomic data from 46 mCRPC patients with Abiraterone-response clinical data and uncovered abrogated immune response regulatory modules that showed strong differential regulation in Abi-resistant patients. These modules were replicated in an independent mCRPC study. Further, we experimentally validated key rewiring predictions and their associated transcription factors. Among them, ELK3, MXD1, and MYB were found to have a differential role in cell survival for Abi-response-specific settings. Moreover, we identified the role of ELK3 in cell migration capacity, which could have direct impact on mCRPC. Collectively, these findings shed light on the underlying regulatory mechanisms driving abiraterone response.

The human angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2) proteins play key roles in the cellular internalization of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the coronavirus responsible for the coronavirus disease of 2019 (COVID-19) pandemic. We set out to functionally characterize the ACE2 and TMPRSS2 protein abundance for variant alleles encoding these proteins that contained non-synonymous single-nucleotide polymorphisms (nsSNPs) in their open reading frames (ORFs). Specifically, a high-throughput assay, deep mutational scanning (DMS), was employed to test the functional implications of nsSNPs, which are variants of uncertain significance in these two genes. Specifically, we used a ‘landing pad’ system designed to quantify the protein expression for 433 nsSNPs that have been observed in the ACE2 and TMPRSS2 ORFs and found that 8 of 127 ACE2, 19 of 157 TMPRSS2 isoform 1 and 13 of 149 TMPRSS2 isoform 2 variant proteins displayed less than ~25% of the wild-type protein expression, whereas 4 ACE2 variants displayed 25% or greater increases in protein expression. As a result, we concluded that nsSNPs in genes encoding ACE2 and TMPRSS2 might potentially influence SARS-CoV-2 infectivity. These results can now be applied to DNA sequence data for patients infected with SARS-CoV-2 to determine the possible impact of patient-based DNA sequence variation on the clinical course of SARS-CoV-2 infection.

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype with low overall survival rates and high molecular heterogeneity; therefore, few targeted therapies are available. The luminal androgen receptor (LAR) is the most consistently identified TNBC subtype, but the clinical utility has yet to be established. Here, we constructed a novel genomic classifier, LAR-Sig, that distinguishes the LAR subtype from other TNBC subtypes and provide evidence that it is a clinically distinct disease. A meta-analysis of seven TNBC datasets (n = 1086 samples) from neoadjuvant clinical trials demonstrated that LAR patients have significantly reduced response (pCR) rates than non-LAR TNBC patients (odds ratio = 2.11, 95% CI: 1.33, 2.89). Moreover, deconvolution of the tumor microenvironment confirmed an enrichment of luminal epithelium corresponding with a decrease in basal and myoepithelium in LAR TNBC tumors. Increased immunosuppression in LAR patients may lead to a decreased presence of cycling T-cells and plasma cells. While, an increased presence of myofibroblast-like cancer-associated cells may impede drug delivery and treatment. In summary, the lower levels of tumor infiltrating lymphocytes (TILs), reduced immune activity in the micro-environment, and lower pCR rates after NAC, suggest that new therapeutic strategies for the LAR TNBC subtype need to be developed.

Alternative polyadenylation (APA) is a common genetic regulatory mechanism that generates distinct 3′ ends for RNA transcripts. Changes in APA have been associated with multiple biological processes and disease phenotypes. However, the role of hormones and their drug analogs in APA remains largely unknown. In this study, we investigated transcriptome-wide the impact of glucocorticoids on APA in 30 human B-lymphoblastoid cell lines. We found that glucocorticoids could regulate APA of a subset of genes, possibly by changing the expression of 142 RNA-binding proteins, some with known APA-regulating properties. Interestingly, genes with glucocorticoid-mediated APA were enriched in viral translation-related pathways, while genes with glucocorticoid-mediated expression were enriched in interferon and interleukin pathways, suggesting that glucocorticoid-mediated APA might result in functional consequences distinct from gene expression. Glucocorticoid-mediated APA was also cell-type-specific, suggesting an action of glucocorticoids that may be unique to immune regulation. We also observed evidence for genotype-dependent glucocorticoid-mediated APA, providing potential functional mechanisms for a series of common genetic variants that had previously been associated with immune disorders, but without a clear mechanism. In summary, this study reports a series of novel observations regarding the impact of glucocorticoids on APA, suggesting novel avenues for mechanistic studies of hormonal biology.

Age at depressive onset (AAO) corresponds to unique symptomatology and clinical outcomes. Integration of genome-wide association study (GWAS) results with additional “omic” measures to evaluate AAO has not been reported and may reveal novel markers of susceptibility and/or resistance to major depressive disorder (MDD). To address this gap, we integrated genomics with metabolomics using data-driven network analysis to characterize and differentiate MDD based on AAO. This study first performed two GWAS for AAO as a continuous trait in (a) 486 adults from the Pharmacogenomic Research Network-Antidepressant Medication Pharmacogenomic Study (PGRN-AMPS), and (b) 295 adults from the Combining Medications to Enhance Depression Outcomes (CO-MED) study. Variants from top signals were integrated with 153 p180-assayed metabolites to establish multi-omics network characterizations of early (

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.

Purpose Pharmacogenomics (PGx) studies how inherited genetic variations in individuals affect drug absorption, distribution, and metabolism. PGx panel testing can potentially help improve efficiency and accuracy in individualizing therapy. This study compared the cost-effectiveness between preemptive PGx panel testing, reactive PGx panel testing and usual care (no testing) in cardiovascular disease management. Methods We developed a decision analytic model from the US payer's perspective for a hypothetical cohort of 10,000 patients ≥45 years old, using a short-term decision tree and long-term Markov model. The testing panel included the following gene-drug pairs: CYP2C19-clopidogrel, CYP2C9/VKORC1-warfarin, and SLCO1B1-statins with 30 test-return days. Costs were reported in 2019 US dollars and effectiveness was measured in quality-adjusted life years (QALYs). The primary outcome was incremental cost-effectiveness ratio (ICER = ΔCost/ΔQALY), assuming 3% discount rate for costs and QALYs. Scenario and probabilistic sensitivity analyses were performed to assess the impact of demographics, risk level, and follow-up timeframe. Results Preemptive testing was found to be cost-effective compared with usual care (ICER $86,227/QALY) at the willingness-to-pay threshold of $100,000/QALY while reactive testing was not (ICER $148,726/QALY). Sensitivity analyses suggested that our cost-effectiveness results were sensitive to longer follow-up, and the age group 45-64 years. Conclusion Compared with usual care, preemptive PGx panel testing was cost-effective in cardiovascular disease management.

Next generation antiandrogens such as enzalutamide (Enz) are effective initially for the treatment of castration-resistant prostate cancer (CRPC). However, the disease often relapses and the underlying mechanisms remain elusive. By performing H3-lysine-27 acetylation (H3K27ac) ChIP-seq in Enz-resistant CRPC cells, we identified a group of super enhancers (SEs) that are abnormally activated in Enz-resistant CRPC cells and associated with enhanced transcription of a subset of tumor promoting genes such as CHPT1, which catalyzes phosphatidylcholine (PtdCho) synthesis and regulates choline metabolism. Increased CHPT1 conferred CRPC resistance to Enz in vitro and in mice. While androgen receptor (AR) primarily binds to a putative CHPT1 enhancer and mediates androgen-dependent expression of CHPT1 gene in Enz-sensitive prostate cancer cells, AR binds to a different enhancer within the CHPT1 SE locus and facilities androgen-independent expression of CHPT1 in Enz-resistant cells. We further identified a long-non coding RNA transcribed at CHPT1 enhancer (also known as enhancer RNA) that binds to the H3K27ac reader BRD4 and participates in regulating CHPT1 SE activity and CHPT1 gene expression. Our findings demonstrate that aberrantly activated SE upregulates CHPT1 expression and confers Enz resistance in CRPC, suggesting that SE-mediated expression of downstream effectors such as CHPT1 can be viable targets to overcome Enz resistance in PCa.

Objectives Based on our previous findings that postmenopausal women with estrone (E1) and estradiol (E2) concentrations at or above 1.3 pg/ml and 0.5 pg/ml, respectively, after 6 months of adjuvant anastrozole therapy had a three-fold risk of recurrence, we aimed to identify a single-nucleotide polymorphism (SNP)-based model that would predict elevated E1 and E2 and then validate it in an independent dataset. Patients and methods The test set consisted of 322 women from the M3 study and the validation set consisted of 152 patients from MA.27. All patients were treated with adjuvant anastrozole, had on-anastrozole E1 and E2 concentrations and genome-wide genotyping. Results SNPs were identified from the M3 genome-wide association study. The best model to predict the E1-E2 phenotype with high balanced accuracy was a support vector machine model using clinical factors plus 46 SNPs. We did not have an independent cohort that is similar to the M3 study with clinical, E1-E2 phenotypes and genotype data to test our model. Hence, we chose a nested matched case-control cohort (MA.27 study) for testing. Our E1-E2 model was not validated but we found the MA.27 validation cohort was both clinically and genomically different. Conclusions We identified a SNP-based model that had excellent performance characteristics for predicting the phenotype of elevated E1 and E2 in women treated with anastrozole. This model was not validated in an independent dataset but that dataset was clinically and genomically substantially different. The model will need validation in a prospective study.

Purpose:To determine if the degree of estrogen suppression with aromatase inhibitors (AI: anastrozole, exemestane, letrozole) is associated with efficacy in early-stage breast cancer, and to examine for differences in the mechanism of action between the three AIs.Experimental Design:Matched case–control studies [247 matched sets from MA.27 (anastrozole vs. exemestane) and PreFace (letrozole) trials] were undertaken to assess whether estrone (E1) or estradiol (E2) concentrations after 6 months of adjuvant therapy were associated with risk of an early breast cancer event (EBCE). Preclinical laboratory studies included luciferase activity, cell proliferation, radio-labeled ligand estrogen receptor binding, surface plasmon resonance ligand receptor binding, and nuclear magnetic resonance assays.Results:Women with E1 ≥1.3 pg/mL and E2 ≥0.5 pg/mL after 6 months of AI treatment had a 2.2-fold increase in risk (P = 0.0005) of an EBCE, and in the anastrozole subgroup, the increase in risk of an EBCE was 3.0-fold (P = 0.001). Preclinical laboratory studies examined mechanisms of action in addition to aromatase inhibition and showed that only anastrozole could directly bind to estrogen receptor α (ERα), activate estrogen response element-dependent transcription, and stimulate growth of an aromatase-deficient CYP19A1−/− T47D breast cancer cell line.Conclusions:This matched case–control clinical study revealed that levels of estrone and estradiol above identified thresholds after 6 months of adjuvant anastrozole treatment were associated with increased risk of an EBCE. Preclinical laboratory studies revealed that anastrozole, but not exemestane or letrozole, is a ligand for ERα. These findings represent potential steps towards individualized anastrozole therapy.

Background Sex differences between psychotic depression (PD) and non-psychotic depression (NPD) have received little systematic study. This study was conducted to investigate sex difference in patients with psychotic and non-psychotic major depressive disorder in a Chinese Han population. Methods In this cross-sectional study, a total of 1718 first-episode and drug-naive outpatients with major depressive disorder were recruited. Demographic and clinical characteristics were collected. All subjects were rated on the Hamilton Depression Rating Scale (HAMD), Hamilton Anxiety Rating Scale (HAMA) and the Positive and Negative Syndrome Scale (PANSS). Results The prevalence of PD in female patients (10.97%) was higher than that in male patients (7.99%). Analysis of variance (ANOVA) showed that female patients were older compared with male patients in NPD group, but there were no significant differences in demographic and clinical variables between female and male PD patients. Further, there were no sex differences in the scores of HAMD, HAMA and positive symptom subscale of PANSS in both PD and NPD groups. Two-way ANOVA showed that PD patients had significantly higher scores on the HAMD, HAMA and positive symptom subscale of PANSS than non-PD patients. However, there were no significant effects of sex and sex* subtypes. Limitations The main limitations are cross-sectional design and inability to control selection bias. Conclusions Our findings show significant differences in clinical profiles between PD and NPD patients; however, no sex difference has been observed in the either PD or NPD patients.

Programmed death ligand-1 (PD-L1 or B7-H1) is well known for its role in immune checkpoint regulation, but its function inside the tumor cells has rarely been explored. Here we report that nuclear PD-L1 is important for cancer cell sister chromatid cohesion. We found that depletion of PD-L1 suppresses cancer cell proliferation, colony formation in vitro, and tumor growth in vivo in immune-deficient NSG mice independent of its role in immune checkpoint. Specifically, PD-L1 functions as a subunit of the cohesin complex, and its deficiency leads to formation of multinucleated cells and causes a defect in sister chromatid cohesion. Mechanistically, PD-L1 compensates for the loss of Sororin, whose expression is suppressed in cancer cells overexpressing PD-L1. PD-L1 competes with Wing Apart-Like (WAPL) for binding to PDS5B, and secures proper sister chromatid cohesion and segregation. Our findings suggest an important role for nuclear PD-L1 in cancer cells independent of its function in immune checkpoint.

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.

Purpose To examine the evidence on the cost-effectiveness of implementing pharmacogenomics (PGx) in cardiovascular disease (CVD) care. Methods We conducted a systematic review using multiple databases from inception to 2018. The titles and abstracts of cost-effectiveness studies on PGx-guided treatment in CVD care were screened, and full texts were extracted. Results We screened 909 studies and included 46 to synthesize. Acute coronary syndrome and atrial fibrillation were the predominantly studied conditions (59%). Most studies (78%) examined warfarin–CYP2C9/VKORC1 or clopidogrel–CYP2C19. A payer’s perspective was commonly used (39%) for cost calculations, and most studies (46%) were US-based. The majority (67%) of the studies found PGx testing to be cost-effective in CVD care, but cost-effectiveness varied across drugs and conditions. Two studies examined PGx panel testing, of which one examined pre-emptive testing strategies. Conclusion We found mixed evidence on the cost-effectiveness of PGx in CVD care. Supportive evidence exists for clopidogrel–CYP2C19 and warfarin–CYP2C9/VKORC1, but evidence is limited in other drug–gene combinations. Gaps persist, including unclear explanation of perspective and cost inputs, underreporting of study design elements critical to economic evaluations, and limited examination of PGx panel and pre-emptive testing for their cost-effectiveness. This review identifies the need for further research on economic evaluations of PGx implementation.

Background: The combination of abiraterone acetate and prednisone (AA/P) is used to treat metastatic prostate cancer, but molecular predictors of treatment response are not well elucidated. We evaluated plasma circulating tumor DNA– (ctDNA-) based copy number alterations (CNAs) to determine treatment-related predictive and prognostic biomarkers for metastatic castration-resistant prostate cancer (mCRPC). Methods: Serial plasma specimens were prospectively collected from 88 chemotherapy-naive mCRPC patients before and after 12 weeks of AA/P treatment. Sequencing-based CNA analyses were performed on 174 specimens. We evaluated CNA-associated 12-week responses for primary resistance, time to treatment change (TTTC) for secondary resistance, and overall survival for prognosis (P < .05). Associations with primary resistance were analyzed using the Fisher exact test. Kaplan–Meier survival curves and Cox regression analyses were used to determine the associations of CNAs with acquired resistance and overall survival. Results: ctDNA reduced by 3.89% in responders and increased by 0.94% in nonresponders (P = .0043). Thirty-one prostate cancer–related genes from whole genome CNAs were tested. AR and AR enhancer amplification were associated with primary resistance (P = .0039) and shorter TTTC (P = .0003). ZFHX3 deletion and PIK3CA amplification were associated with primary resistance (P = .026 and P = .017, respectively), shorter TTTC (P = .0008 and P= .0016, respectively), and poor survival (P = .0025 and P = .0022, respectively). CNA-based risk scores combining selected significant associations (AR, NKX3.1, and PIK3CA) at the univariate level with TTTC were predictive of secondary resistance (P = .0002). and established prognoses for survival based on CNAs in ZFHX3, RB1, PIK3CA, and OPHN1 (P = .002). Multigene risk scores were more predictive than individual genes or clinical risk factors (P < .05). Conclusion: Plasma ctDNA CNAs and risk scores can predict mCRPC-state treatment and survival outcomes.

Background: The combination of cyclin dependent 4/6 kinase inhibitors (CDK4/6i) with endocrine therapy (ET) has resulted in clinically significant improvements in progression-free survival (PFS) and overall survival (OS) in hormone-receptor (HR)-positive metastatic breast cancer (MBC). However, most patients’ disease ultimately progresses on CDK4/6i and ET. Therefore, further research is necessary to understand the mechanisms driving primary and secondary resistance. PROMISE is a multicenter prospective cohort study enrolling women with HR-positive MBC commencing treatment with palbociclib + letrozole (1st line) or palbociclib + fulvestrant (2nd Line). The study provides a comprehensive “omic” assessment of blood, tumor, urine and the fecal microbiome to identify molecular or cellular features associated with primary endocrine resistance (e.g. disease progression ≤ 12 months on treatment) and acquired resistance to CDK 4/6i. Additionally, patient derived xenografts and organoids are created to test new drug strategies designed to overcome resistance to CDK 4/6i and ET. Here, we present initial sequencing results from pretreatment biospecimens collected from PROMISE study participants. Methods: On-study tumor biopsies and blood samples were collected for DNA/RNA sequencing (TempusTM). The analyzed biospecimens were all obtained prior to initiation of palbociclib and ET. We correlated patient clinical characteristics (phenotypes) with molecular data and responses to protocol treatment. The data were analyzed using a series of cutting-edge bioinformatics pipelines for somatic and germline mutations in addition to copy number alterations (CNAs). The study database was locked for analysis on 06/20/2019. Results: We analyzed the somatic single nucleotide variants/INDELs (sSNV/INDEL) profiles across the tumor samples to determine the genes that were least likely to occur as a result of background mutation processes. Twenty-six patients had somatic copy number alterations (sCNA) and/or sSNV/INDEL in at least one of 18 genes with the most significant sSNV/INDEL profiles (p < 0.03) which included clinically and biologically relevant genes. The genes with the most statistically significant sSNV/INDEL mutation profiles were GATA3, PIK3CA, CDH1, and ESR1 (p < 0.0009). We observed a high percentage of tumors with somatic alterations in GATA3 (23% sSNV/INDEL, 15% sCNA), PIK3CA (38%, 12%), CDH1 (19%, 50%) and ESR1 (19%, 58%). ESR1 mutations were more frequent in patients receiving 2nd line treatment. Other frequently altered genes included TP53 (15%, 46%), MAP2K4 (8%, 50%), DNAAF1 (8%, 50%), and CDKN1B (8%, 35%). Further, ZNF317 and F3 were altered in 9 and 7 patients, respectively. Twenty-four samples had alterations in at least one of the CDK4/6 pathway genes (RB1, CCNE2, CCND1, CDK6, ESR1, CDKN2A, CCND3, CDK4, CDK2 and CCNE1). Four patients progressed on therapy; three of the four patients had mutations in PIK3CA, and one had a mutation in ESR1. Results of the RNA sequencing data (N=26) will be presented at the SABCS meeting. Conclusions: This is the first report of a prospective study designed to characterize the genomic landscape of ER+/HER2- MBC prior to palbociclib treatment. We observed high frequencies of known targetable alterations in PIK3CA and ESR1, including in patients that progress, which is consistent with previous reports. RNA sequencing data will be presented at the meeting. Citation Format: Ciara C O Sullivan, Krishna R Kalari, Vera J Suman, Peter T Vedell, Ann Moyer, Erin Carlson, Jason Sinnwell, Tejaswi Alaparthi, Xiaojia Tang, Kevin Thompson, Jaeyun Sung, Alvaro Moreno-Aspitia, Donald Northfelt, Minetta C Liu, Tufia C Haddad, Prema Peethambaram, Saranya Chumsri, Kathryn J Ruddy, Karthik V Giridhar, Roberto A Leon-Ferre, Paula Gill, Mohammad Ranginwala, Asad Javed, Sameer Batoo, Brendan P. McMenomy, Richard Weinshilboum, Liewei Wang, Matthew P Goetz. Comprehensive tumor sequencing to identify biomarkers of palbociclib response: First report of the PROMISE study [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P2-11-07.

A piezoelectric valveless pump with flexible actuators (PVPFA) is proposed, which owns advantages of flexible driving, changeable flow direction and high output performance. The PVPFA is consist of three flexible actuators and a pump body. The direction of pumping fluid can be changed by setting different fluid medium volumes in the hydraulic amplification chamber. Each flexible actuator adopts hydraulic amplification structure, which makes PVPFA provide higher output performance. The theoretical model of the factors affecting the flow rate of PVPFA is established and verified by experiments. The experimental results show that PVPFA obtains a maximum flow rate of 98.48 ml/min and a maximum pressure of 10.7 kPa. PVPFA provides a new idea for piezoelectric valveless pump.

PURPOSE: The Mayo-Baylor RIGHT 10K Study enabled preemptive, sequence-based pharmacogenomics (PGx)-driven drug prescribing practices in routine clinical care within a large cohort. We also generated the tools and resources necessary for clinical PGx implementation and identified challenges to be overcome. Furthermore, we measured the frequency of both common genetic variation for which clinical guidelines already exist and that of rare variation that could be detected by DNA sequencing, rather than genotyping. METHODS: Targeted oligonucleotide-capture sequencing of 77 “pharmacogenes” was performed using DNA from 10,077 consented Mayo Clinic Biobank volunteers. The resulting predicted drug response related phenotypes for 13 genes, including CYP2D6 and HLA, affecting 21 drug-gene pairs, were deposited preemptively in the Mayo electronic health record (EHR). RESULTS: For the 13 pharmacogenes of interest, the genomes of 79% of participants carried clinically actionable variants in 3 or more genes and DNA sequencing identified an average of 3.3 additional conservatively predicted deleterious variants that would not have been evident using genotyping. CONCLUSIONS: Implementation of preemptive rather than reactive, sequence-based rather than genotype-based PGx prescribing revealed nearly universal patient applicability and required integrated institution-wide resources to fully realize individualized drug therapy and to demonstrate more efficient use of health care resources.

Understanding the function of non-coding genetic variants represents a formidable challenge for biomedicine. We previously identified genetic variants that influence gene expression only after exposure to a hormone or drug. Using glucocorticoid signaling as a model system, we have now demonstrated, in a genome-wide manner, that exposure to glucocorticoids triggered disease risk variants with previously unclear function to influence the expression of genes involved in autoimmunity, metabolic and mood disorders, osteoporosis and cancer. Integrating a series of pharmacogenomic and pharmacoepigenomic datasets, we identified the cis-regulatory elements and 3-dimensional interactions underlying the ligand-dependent associations between those genetic variants and distant risk genes. These observations increase our understanding of mechanisms of non-coding genetic variant-chemical environment interactions and advance the fine-mapping of disease risk and pharmacogenomic loci.

Combination antidepressant pharmacotherapies are frequently used to treat major depressive disorder (MDD). However, there is no evidence that machine learning approaches combining multi-omics measures (e.g., genomics and plasma metabolomics) can achieve clinically meaningful predictions of outcomes to combination pharmacotherapy. This study examined data from 264 MDD outpatients treated with citalopram or escitalopram in the Mayo Clinic Pharmacogenomics Research Network Antidepressant Medication Pharmacogenomic Study (PGRN-AMPS) and 111 MDD outpatients treated with combination pharmacotherapies in the Combined Medications to Enhance Outcomes of Antidepressant Therapy (CO-MED) study to predict response to combination antidepressant therapies. To assess whether metabolomics with functionally validated single-nucleotide polymorphisms (SNPs) improves predictability over metabolomics alone, models were trained/tested with and without SNPs. Models trained with PGRN-AMPS’ and CO-MED’s escitalopram/citalopram patients predicted response in CO-MED’s combination pharmacotherapy patients with accuracies of 76.6% (p p p p

Clinical pharmacogenomic testing typically uses targeted genotyping, which only detects variants included in the test design and may vary among laboratories. To evaluate the potential patient impact of genotyping compared with sequencing, which can detect common and rare variants, an in silico targeted genotyping panel was developed based on the variants most commonly included in clinical tests and applied to a cohort of 10,030 participants who underwent sequencing for CYP1A2, CYP2C19, CYP2C9, CYP2D6, CYP3A4, CYP3A5, DPYD, SLCO1B1, TPMT, UGT1A1, and VKORC1. The results of in silico targeted genotyping were compared with the clinically reported sequencing results. Of the 10,030 participants, 2780 (28%) had at least one potentially clinically relevant variant/allele identified by sequencing that would not have been detected in a standard targeted genotyping panel. The genes with the largest number of participants with variants only detected by sequencing were SLCO1B1, DPYD, and CYP2D6, which affected 13%, 6.3%, and 3.5% of participants, respectively. DPYD (112 variants) and CYP2D6 (103 variants) had the largest number of unique variants detected only by sequencing. Although targeted genotyping detects most clinically significant pharmacogenomic variants, sequencing-based approaches are necessary to detect rare variants that collectively affect many patients. However, efforts to establish pharmacogenomic variant classification systems and nomenclature to accommodate rare variants will be required to adopt sequencing-based pharmacogenomics.

Triggering receptor expressed on myeloid cell 2 (TREM2) is a surface receptor that, in the central nervous system, is exclusively expressed on microglia. TREM2 variants have been linked to increased risk for neurodegenerative diseases, but the functional effects of microglial TREM2 remain largely unknown. To this end, we investigated TAR-DNA binding protein 43 kDa (TDP-43)-related neurodegenerative disease via viral-mediated expression of human TDP-43 protein (hTDP-43) in neonatal and adult mice or inducible expression of hTDP43 with defective nuclear localization signals in transgenic mice. We found that TREM2 deficiency impaired microglia phagocytic clearance of pathological TDP-43, and enhanced neuronal damage and motor function impairments. Mass cytometry analysis revealed that hTDP-43 induced a TREM2-dependent subpopulation of microglia with high CD11c expression and higher phagocytic ability. Using mass spectrometry and surface plasmon resonance analysis, we further demonstrated an interaction between TDP-43 and TREM2,in vitroandin vivo,in hTDP-43-expressing transgenic mouse brains. We computationally identified the region within hTDP-43 that interacts with TREM2 and observed the potential interaction in ALS patient tissues. Our data reveal the novel interaction between TREM2 and TDP-43, highlighting that TDP-43 is a possible ligand for microglial TREM2 and the interaction mediates neuroprotection of microglial TREM2 in TDP-43-related neurodegeneration.

Androgen receptor (AR)-positive prostate cancers (PCa) and estrogen receptor (ER)-positive luminal breast cancers (BCa) are generally less responsive to immunotherapy compared to certain tumor types such as melanoma. However, the underlying mechanisms are not fully elucidated. Here we found that FOXA1 overexpression inversely correlated with interferon (IFN) signature and antigen presentation gene expression in PCa and BCa patients. FOXA1 bound STAT2 DNA binding domain and suppressed STAT2 DNA binding activity, IFN signaling gene expression and cancer immune response independently of the transactivation activity of FOXA1 and its mutations detected in prostate and breast cancers. Increased FOXA1 expression promoted cancer immuno- and chemotherapy resistance in mice and PCa and BCa patients. These findings were also validated in bladder cancer expressing high level FOXA1. FOXA1 overexpression could be a prognostic factor to predict therapy resistance and a viable target to sensitize luminal prostate, breast and bladder cancer to immuno- and chemotherapy.

BackgroundPrediction of the response of cancer patients to different treatments and identification of biomarkers of drug sensitivity are two major goals of individualized medicine. In this study, we developed a deep learning framework called TINDL, completely trained on preclinical cancer cell lines, to predict the response of cancer patients to different treatments. TINDL utilizes a tissue-informed normalization to account for the tissue and cancer type of the tumours and to reduce the statistical discrepancies between cell lines and patient tumours. In addition, this model identifies a small set of genes whose mRNA expression are predictive of drug response in the trained model, enabling identification of biomarkers of drug sensitivity.ResultsUsing data from two large databases of cancer cell lines and cancer tumours, we showed that this model can distinguish between sensitive and resistant tumours for 10 (out of 14) drugs, outperforming various other machine learning models. In addition, our siRNA knockdown experiments on 10 genes identified by this model for one of the drugs (tamoxifen) confirmed that all of these genes significantly influence the drug sensitivity of the MCF7 cell line to this drug. In addition, genes implicated for multiple drugs pointed to shared mechanism of action among drugs and suggested several important signaling pathways.ConclusionsIn summary, this study provides a powerful deep learning framework for prediction of drug response and for identification of biomarkers of drug sensitivity in cancer.

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.

Triple-negative breast cancer (TNBC) accounts for 15–20% of all invasive breast cancers and tends to have aggressive histological features and poor clinical outcomes. Unlike, estrogen receptor- or HER2-positive diseases, TNBC patients currently lack the US FDA-approved targeted therapies. DNA methylation is a critical mechanism of epigenetic modification. It is well known that aberrant DNA methylation contributes to the malignant transformation of cells by silencing critical tumor suppressor genes. DNA methyltransferase inhibitors reactivate silenced tumor suppressor genes and result in tumor growth arrest, with therapeutic effects observed in patients with hematologic malignancies. The antitumor effect of these DNA methyltransferase inhibitors has also been explored in solid tumors, especially in TNBC that currently lacks targeted therapies.

OBJECTIVE. The purpose of this study is to prospectively compare the size of invasive breast cancer before and after neoadjuvant chemotherapy (NAC) at breast MRI and molecular breast imaging (MBI) and to assess the accuracy of post-NAC MBI and MRI relative to pathologic analysis. SUBJECTS AND METHODS. Women with invasive breast cancer greater than or equal to 1.5 cm were enrolled to compare the longest dimension before and after NAC at MRI and MBI. MBI was performed on a dual-detector cadmium zinc telluride system after administration of 6.5 mCi (240 MBq) 99mTc-sestamibi. The accuracy of MRI and MBI in assessing residual disease (invasive disease or ductal carcinoma in situ) was determined relative to pathologic examination. RESULTS. The longest dimension at MRI was within 1.0 cm of that at MBI in 72.3% of cases before NAC and 70.1% of cases after NAC. The difference between the longest dimension at imaging after NAC and pathologic tumor size was within 1 cm for 58.7% of breast MRI cases and 59.6% of MBI cases. Ninety patients underwent both MRI and MBI after NAC. In the 56 patients with invasive residual disease, 10 (17.9%) cases were negative at MRI and 23 (41.1%) cases were negative at MBI. In the 34 patients with breast pathologic complete response, there was enhancement in 10 cases (29.4%) at MRI and uptake in six cases (17.6%) at MBI. Sensitivity, specificity, positive predictive value, and negative predictive value after NAC were 82.8%, 69.4%, 81.4%, and 71.4%, respectively, for MRI and 58.9%, 82.4%, 84.6%, and 54.9%, respectively, for MBI. CONCLUSION. Breast MRI and MBI showed similar disease extent before NAC. MBI may be an alternative to breast MRI in patients with a contraindication to breast MRI. Neither modality showed sufficient accuracy after NAC in predicting breast pathologic complete response to obviate tissue diagnosis to assess for residual invasive disease. Defining the extent of residual disease compared with pathologic evaluation was also limited after NAC for both breast MRI and MBI.

Greater than 90% of significant genome-wide association study (GWAS) single-nucleotide polymorphisms (SNPs) are in noncoding regions of the genome, but only 25.6% are known expression quantitative trait loci (eQTLs). Therefore, the function of many significant GWAS SNPs remains unclear. We have identified a novel type of eQTL for which SNPs distant from ligand-activated transcription factor (TF) binding sites can alter target gene expression in a SNP genotype-by-ligand-dependent fashion that we refer to as pharmacogenomic eQTLs (PGx-eQTLs)-loci that may have important pharmacotherapeutic implications. In the present study, we integrated chromatin immunoprecipitation-seq with RNA-seq and SNP genotype data for a panel of lymphoblastoid cell lines to identify 10 novel cis PGx-eQTLs dependent on the ligand-activated TF aryl hydrocarbon receptor (AHR)-a critical environmental sensor for xenobiotic (drug) and immune response. Those 10 cis PGx-eQTLs were eQTLs only after AHR ligand treatment, even though the SNPs did not create/destroy an AHR response element-the DNA sequence motif recognized and bound by AHR. Additional functional studies in multiple cell lines demonstrated that some cis PGx-eQTLs are functional in multiple cell types, whereas others displayed SNP-by-ligand-dependent effects in just one cell type. Furthermore, four of those cis PGx-eQTLs had previously been associated with clinical phenotypes, indicating that those loci might have the potential to inform clinical decisions. Therefore, SNPs across the genome that are distant from TF binding sites for ligand-activated TFs might function as PGx-eQTLs and, as a result, might have important clinical implications for interindividual variation in drug response. SIGNIFICANCE STATEMENT: More than 90% of single-nucleotide polymorphisms (SNPs) that are associated with clinical phenotypes are located in noncoding regions of the genome. However, the mechanisms of action of many of those SNPs have not been elucidated, and drugs may unmask functional expression quantitative trail loci (eQTLs). In the current study, we used drugs that bind to the ligand-activated transcription factor aryl hydrocarbon receptor (AHR) and identified SNPs that were associated with interindividual variation in gene expression following drug exposure-termed pharmacogenomic (PGx)-eQTLs. Possibly of greater significance, those PGx-eQTL SNPs were outside of AHR binding sites, indicating that they do not interrupt AHR DNA recognition. PGx-eQTLs such as those described in this work may have crucial implications for interindividual variation in drug.

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 (HRSD17) 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 HRSD17 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 HRSD17. 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 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 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.

Patient-derived xenografts (PDX) are generated in immune deficient mice and demonstrate histologic and molecular features similar to their corresponding human tumors. However, murine tumors (non-human) spontaneously occur in these models. 120 consecutive patients with high-risk primary breast cancer enrolled in the prospective neoadjuvant BEAUTY study had tumor tissue obtained at the time of diagnosis. These tumor cells, including initial tissue and subsequent generations, were injected into either NSG (n = 365) or NOD-SCID (n = 396) female mice. Mice with initial tumor growth sufficient for transfer to the 2nd generation underwent histologic review by pathologists, including Ki67 staining. After passaging the tumors for up to 4 generations, at least one primary mouse tumor was detected from 24 of the 54 PDX-lines, for a frequency of 3.2% (24 mice out of 761 mice), including murine lymphomas (n = 13), mammary tumors (n = 7), osteosarcomas (n = 2), and hemangiosarcomas (n = 2). While true PDX showed scattered strong staining with Ki67, murine tumors were Ki67 negative. No significant differences (p = 0.062) were observed comparing development of murine tumors in NOD-SCID (n = 8) vs NSG mice (n = 16). While PDX are a useful tool in cancer research, there is a potential for spontaneous murine tumors to arise, which could alter results of studies utilizing PDX. Morphologic review by a pathologist, potentially along with Ki67 staining, is necessary to ensure that tumor growth represents the desired PDX prior to use in downstream studies. This study is the first prospective study evaluating the frequency, type, and time frame for development of non-human tumors.

Background Our previous genome-wide association study using the MA.27 aromatase inhibitors adjuvant trial identified SNPs in the long noncoding RNA MIR2052HG associated with breast cancer-free interval. MIR2052HG maintained ERα both by promoting AKT/FOXO3-mediated ESR1 transcription and by limiting ubiquitin-mediated ERα degradation. Our goal was to further elucidate MIR2052HG’s mechanism of action. Methods RNA-binding protein immunoprecipitation assays were performed to demonstrate that the transcription factor, early growth response protein 1 (EGR1), worked together with MIR2052HG to regulate that lemur tyrosine kinase-3 (LMTK3) transcription in MCF7/AC1 and CAMA-1 cells. The location of EGR1 on the LMTK3 gene locus was mapped using chromatin immunoprecipitation assays. The co-localization of MIR2052HG RNA and the LMTK3 gene locus was determined using RNA-DNA dual fluorescent in situ hybridization. Single-nucleotide polymorphisms (SNP) effects were evaluated using a panel of human lymphoblastoid cell lines. Results MIR2052HG depletion in breast cancer cells results in a decrease in LMTK3 expression and cell growth. Mechanistically, MIR2052HG interacts with EGR1 and facilitates its recruitment to the LMTK3 promoter. LMTK3 sustains ERα levels by reducing protein kinase C (PKC) activity, resulting in increased ESR1 transcription mediated through AKT/FOXO3 and reduced ERα degradation mediated by the PKC/MEK/ERK/RSK1 pathway. MIR2052HG regulated LMTK3 in a SNP- and aromatase inhibitor-dependent fashion: the variant SNP increased EGR1 binding to LMTK3 promoter in response to androstenedione, relative to wild-type genotype, a pattern that can be reversed by aromatase inhibitor treatment. Finally, LMTK3 overexpression abolished the effect of MIR2052HG on PKC activity and ERα levels. Conclusions Our findings support a model in which the MIR2052HG regulates LMTK3 via EGR1, and LMTK3 regulates ERα stability via the PKC/MEK/ERK/RSK1 axis. These results reveal a direct role of MIR2052HG in LMTK3 regulation and raise the possibilities of targeting MIR2052HG or LMTK3 in ERα-positive breast cancer. Electronic supplementary material The online version of this article (10.1186/s13058-019-1130-3) contains supplementary material, which is available to authorized users.

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.

CYP2C9 and CYP2C19 are highly polymorphic pharmacogenes; however, clinically actionable genetic variability in drug metabolism due to these genes has been limited to a few common alleles. The identification and functional characterization of less-common open reading frame sequence variation might help to individualize therapy with drugs that are substrates for the enzymes encoded by these genes. The present study identified seven uncharacterized variants each in CYP2C9 and CYP2C19 using next-generation sequence data for 1013 subjects, and functionally characterized the encoded proteins. Constructs were created and transiently expressed in COS-1 cells for the assay of protein concentration and enzyme activities using fluorometric substrates and liquid chromatography- tandem mass spectrometry with tolbutamide (CYP2C9) and (S)-mephenytoin (CYP2C19) as prototypic substrates. The results were compared with the SIFT, Polyphen, and Provean functional prediction software programs. Cytochrome P450 oxidoreductase (CPR) activities were also determined. Positive correlations were observed between protein content and fluorometric enzyme activity for variants of CYP2C9 (P C and CYP2C19 65A>G activities were much lower than predicted based on protein content. Substrate intrinsic clearance values for CYP2C9 218C>T, 343A>C, and CYP2C19 337G>A, 518C>T, 556C>T, and 557G>A were less than 25% of wild-type allozymes. CPR activity levels were similar for all variants. In summary, sequencing of CYP2C9 and CYP2C19 in 1013 subjects identified low-frequency variants that had not previously been functionally characterized. In silico predictions were not always consistent with functional assay results. These observations emphasize the need for high-throughput methods for pharmacogene variant mutagenesis and functional characterization.

Selective serotonin reuptake inhibitors (SSRIs) are first-line antidepressants for the treatment of major depressive disorder (MDD). However, treatment response during an initial therapeutic trial is often poor and is difficult to predict. Heterogeneity of response to SSRIs in depressed patients is partly driven by co-occurring somatic disorders such as coronary artery disease (CAD) and obesity. CAD and obesity may also be associated with metabolic side effects of SSRIs. In this study, we assessed the association of CAD and obesity with treatment response to SSRIs in patients with MDD using a polygenic score (PGS) approach. Additionally, we performed cross-trait meta-analyses to pinpoint genetic variants underpinnings the relationship of CAD and obesity with SSRIs treatment response. First, PGSs were calculated at different p value thresholds (PT) for obesity and CAD. Next, binary logistic regression was applied to evaluate the association of the PGSs to SSRIs treatment response in a discovery sample (ISPC, N = 865), and in a replication cohort (STAR*D, N = 1,878). Finally, a cross-trait GWAS meta-analysis was performed by combining summary statistics. We show that the PGSs for CAD and obesity were inversely associated with SSRIs treatment response. At the most significant thresholds, the PGS for CAD and body mass index accounted 1.3%, and 0.8% of the observed variability in treatment response to SSRIs, respectively. In the cross-trait meta-analyses, we identified (1) 14 genetic loci (including NEGR1, CADM2, PMAIP1, PARK2) that are associated with both obesity and SSRIs treatment response; (2) five genetic loci (LINC01412, PHACTR1, CDKN2B, ATXN2, KCNE2) with effects on CAD and SSRIs treatment response. Our findings implicate that the genetic variants of CAD and obesity are linked to SSRIs treatment response in MDD. A better SSRIs treatment response might be achieved through a stratified allocation of treatment for MDD patients with a genetic risk for obesity or CAD.

An amendment to this paper has been published and can be accessed via the original article.

Aromatase inhibitors (AIs) are the treatment of choice for hormone receptor-positive early breast cancer in postmenopausal women. None of the third-generation AIs are superior to the others in terms of efficacy. We attempted to identify genetic factors that could differentiate between the effectiveness of adjuvant anastrozole and exemestane by examining single-nucleotide polymorphism (SNP)-treatment interaction in 4,465 patients. A group of SNPs were found to be differentially associated between anastrozole and exemestane regarding outcomes. However, they showed no association with outcome in the combined analysis. We followed up common SNPs near LY75 and GPR160 that could differentiate anastrozole from exemestane efficacy. LY75 and GPR160 participate in epithelial-to-mesenchymal transition and growth pathways, in both cases with SNP-dependent variation in regulation. Collectively, these studies identified SNPs that differentiate the efficacy of anastrozole and exemestane and they suggest additional genetic biomarkers for possible use in selecting an AI for a given patient.

Resistance to next-generation anti-androgen enzalutamide (ENZ) constitutes a major challenge for the treatment of castration-resistant prostate cancer (CRPC). By performing genome-wide ChIP-seq profiling in ENZ-resistant CRPC cells we identify a set of androgen receptor (AR) binding sites with increased AR binding intensity (ARBS-gained). While ARBS-gained loci lack the canonical androgen response elements (ARE) and pioneer factor FOXA1 binding motifs, they are highly enriched with CpG islands and the binding sites of unmethylated CpG dinucleotide-binding protein CXXC5 and the partner TET2. RNA-seq analysis reveals that both CXXC5 and its regulated genes including ID1 are upregulated in ENZ-resistant cell lines and these results are further confirmed in patient-derived xenografts (PDXs) and patient specimens. Consistent with the finding that ARBS-gained loci are highly enriched with H3K27ac modification, ENZ-resistant PCa cells, organoids, xenografts and PDXs are hyper-sensitive to NEO2734, a dual inhibitor of BET and CBP/p300 proteins. These results not only reveal a noncanonical AR function in acquisition of ENZ resistance, but also posit a treatment strategy to target this vulnerability in ENZ-resistant CRPC., Resistance to second generation anti-androgen therapies such as enzalutamide (ENZ) can emerge in prostate cancer patients. Here, the authors identify an ENZ-resistant mechanism driven by AR-dependent transcription of noncanonical targets that make resistant cells susceptible to dual inhibition of BET and CBP/p300 signaling.

Purpose: Patient-derived xenografts (PDX) are a research tool for studying cancer biology and drug response phenotypes. While engraftment rates are higher for tumors with more aggressive characteristics, it is uncertain whether engraftment is prognostic for cancer recurrence. Patients and Methods: In a prospective study of patients with breast cancer treated with neoadjuvant chemotherapy (NAC) with taxane ± trastuzumab followed by anthracycline-based chemotherapy, we report the association between breast cancer events and PDX engraftment using tumors derived from treatment naïve (pre-NAC biopsies from 113 patients) and treatment resistant (post-NAC at surgery from 34 patients). Gray test was used to assess whether the cumulative incidence of a breast cancer event differs with respect to either pre-NAC PDX engraftment or post-NAC PDX engraftment. Results: With a median follow-up of 5.7 years, the cumulative incidence of breast cancer relapse did not differ significantly according to pre-NAC PDX engraftment (5-year rate: 13.6% vs. 13.4%; P = 0.89). However, the incidence of a breast event was greater for patients with post-NAC PDX engraftment (5-year rate: 50.0% vs. 19.6%), but this did not achieve significance (P = 0.11). Conclusions: In treatment-naïve breast cancer receiving standard NAC, PDX engraftment was not prognostic for breast cancer recurrence. Further study is needed to establish whether PDX engraftment in the treatment-resistant setting is prognostic for cancer recurrence.