Wang, Liewei

AbstractBackground:Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype. Patients with TNBC are primarily treated with neoadjuvant chemotherapy (NAC). Response to NAC is prognostic, with reductions in low overall survival and disease-free survival rates in those patients who do not achieve pathological complete response (pCR). Based on this premise, we focused on patients with residual disease and hypothesized that paired analysis of primary and residual TNBC tumors following NAC could identify unique biomarkers associated with post-NAC recurrence.Methods and Results: We investigated 24 TNBC tumors from 12 patients with pre- and post-NAC data, including four patients who recurred shortly after surgery(<24 months) and eight who remained recurrence-free (> 48 months). These tumors were collected from a prospective NAC breast cancer study (BEAUTY) conducted at the Mayo Clinic. Differential expression analysis of pre-NAC biopsies showed minimal gene expression differences between early recurrent and non-recurrent TNBC tumors; however, post-NAC samples demonstrated a significant alteration in expression patterns in response to intervention. Topological-level differences associated with early recurrence were implicated in 251 gene sets, and an independent assessment of microarray gene expression data from the paired NAC I-SPY1 trial confirmed 56 genesets. Underlying these 56 gene sets, 113 genes were observed to be differentially expressed in the I-SPY1 and BEAUTY post-NAC studies. An independent (n=392) breast cancer data set with relapse-free survival (RFS) data was used to refine our gene list to a 17 gene signature. A three-fold cross-validation analysis of the gene signature with the combined BEAUTY and I-SPY1 data yielded an average AUC of 0.88 for six machine learning models. Due to the limited studies with pre-and post-NAC TNBC tumor data, further validation of the signature is needed.Conclusion:Analysis of multi-omics data from post-NAC TNBC chemoresistant tumors showed downregulation of mismatch repair, purin, and tubulin pathways. Additionally, weidentified a 17 gene signature in TNBC associated with post-NAC recurrence enriched with down-regulation of immune genes.

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

Abstract Purpose: To identify molecular predictors of grade 3/4 neutropenic or leukopenic events (NLE) after chemotherapy using a genome-wide association study (GWAS). Experimental Design: A GWAS was performed on patients in the phase III chemotherapy study SUCCESS-A (n = 3,322). Genotyping was done using the Illumina HumanOmniExpress-12v1 array. Findings were functionally validated with cell culture models and the genotypes and gene expression of possible causative genes were correlated with clinical treatment response and prognostic outcomes. Results: One locus on chromosome 16 (rs4784750; NLRC5; P = 1.56E-8) and another locus on chromosome 13 (rs16972207; TNFSF13B; P = 3.42E-8) were identified at a genome-wide significance level. Functional validation revealed that expression of these two genes is altered by genotype-dependent and chemotherapy-dependent activity of two transcription factors. Genotypes also showed an association with disease-free survival in patients with an NLE. Conclusions: Two loci in NLRC5 and TNFSF13B are associated with NLEs. The involvement of the MHC I regulator NLRC5 implies the possible involvement of immuno-oncological pathways.

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

Abstract Background The AMP‐activated protein kinase (AMPK) is a central regulator of energy homeostasis, with deregulation leading to cancer and other diseases. However, how this pathway is dysregulated in cancer has not been well clarified. Methods Using a tandem affinity purification/mass‐spec technique and biochemical analyses, we identified tumor protein D52 (TPD52) as an AMPKα‐interacting molecule. To explore the biological effects of TPD52 in cancers, we conducted biochemical and metabolic assays in vitro and in vivo with cancer cells and TPD52 transgenic mice. Finally, we assessed the clinical significance of TPD52 expression in breast cancer patients using bioinformatics techniques. Results TPD52, initially identified to be overexpressed in many human cancers, was found to form a stable complex with AMPK in cancer cells. TPD52 directly interacts with AMPKα and inhibits AMPKα kinase activity in vitro and in vivo. In TPD52 transgenic mice, overexpression of TPD52 leads to AMPK inhibition and multiple metabolic defects. Clinically, high TPD52 expression predicts poor survival of breast cancer patients. Conclusion The findings revealed that TPD52 is a novel regulator of energy stress‐induced AMPK activation and cell metabolism. These results shed new light on AMPK regulation and understanding of the etiology of cancers with TPD52 overexpression.

Introduction: Estrogens have been associated with risk of stroke in both men and women. Estradiol and estrone are synthesized from testosterone via the aromatase enzyme in both men and postmenopausal women. Genetic variation affecting function of the aromatase gene (CYP19A1) has been linked to elevated estradiol levels. We tested the hypothesis that these variants, as well as a modifying gene, are associated with ischemic stroke in three cohorts: the Vitamin Intervention for Stroke Prevention trial (VISP), Women’s Health Initiative (WHI), and Genetic Susceptibility to Early-Onset Stroke study (GEOS). Methods: Genomic data for individuals of European descent with ischemic stroke was compared against an equal number of controls: VISP (1723 cases - 65% men), WHI (436 cases - all postmenopausal women), GEOS (557 cases - all men age <55). Genotyping was funded through the Genomics and Randomized Trials Network (GARNET) and Gene Environment Association Studies (GENEVA) consortia. All cohorts used the Illumina Omni1-Quad chip for genotyping with imputation as needed. Eight candidate single nucleotide polymorphisms (SNPs) in CYP19A1 on chromosome 15 and an additional regulatory SNP in a modifier locus on chromosome 8 were chosen for analysis due to coding status or previous association with estradiol levels. Analyses were adjusted for age, hormone therapy use (women), and the top three principal components. Association results were determined as odds ratios for individual SNPs (effect allele = minor), with significance level of p<0.05. Results: GWAS results are tabulated. There were no significant differences in allele frequency for four SNPs in CYP19A1 among stroke patients and controls. The additional four SNPs were monomorphic in all three populations. The modifier gene SNP (rs1864729) was less common among young men with prior stroke in GEOS, but not men in VISP or in women. Conclusions: Aromatase genetic variation was not associated with risk of stroke in these three cohorts. One SNP that modifies aromatase function was less common among young, but not older men with stroke. Replication in other studies is necessary.

Approximately one‐third of patients with metastatic castration‐resistant prostate cancer (CRPC) exhibited primary abiraterone resistance. To identify alternative treatment for abiraterone nonresponders, we performed drug discovery analyses using the L1000 database using differentially expressed genes identified in tumor biopsies and patient‐derived xenograft (PDX) tumors between abiraterone responders and nonresponders enrolled in PROMOTE trial. This approach identified 3 drugs, including topoisomerase II (TOP2) inhibitor mitoxantrone, CDK4/6 inhibitor palbociclib, and pan‐CDK inhibitor PHA‐793887. These drugs significantly suppressed the growth of abiraterone‐resistant cell lines and PDX models. Moreover, we identified 11 genes targeted by all 3 drugs that were associated with worse outcomes in both the PROMOTE and Stand Up To Cancer cohorts. This 11‐gene panel might also function as biomarkers to select the 3 alternative therapies for this subgroup of patients with CRPC, warranting further clinical investigation.

Background The treatment of depression in children and adolescents is a substantial public health challenge. This study examined artificial intelligence tools for the prediction of early outcomes in depressed children and adolescents treated with fluoxetine, duloxetine, or placebo. Methods The study samples included training datasets ( N = 271) from patients with major depressive disorder (MDD) treated with fluoxetine and testing datasets from patients with MDD treated with duloxetine ( N = 255) or placebo ( N = 265). Treatment trajectories were generated using probabilistic graphical models (PGMs). Unsupervised machine learning identified specific depressive symptom profiles and related thresholds of improvement during acute treatment. Results Variation in six depressive symptoms (difficulty having fun, social withdrawal, excessive fatigue, irritability, low self‐esteem, and depressed feelings) assessed with the Children’s Depression Rating Scale‐Revised at 4–6 weeks predicted treatment outcomes with fluoxetine at 10–12 weeks with an average accuracy of 73% in the training dataset. The same six symptoms predicted 10–12 week outcomes at 4–6 weeks in (a) duloxetine testing datasets with an average accuracy of 76% and (b) placebo‐treated patients with accuracies of 67%. In placebo‐treated patients, the accuracies of predicting response and remission were similar to antidepressants. Accuracies for predicting nonresponse to placebo treatment were significantly lower than antidepressants. Conclusions PGMs provided clinically meaningful predictions in samples of depressed children and adolescents treated with fluoxetine or duloxetine. Future work should augment PGMs with biological data for refined predictions to guide the selection of pharmacological and psychotherapeutic treatment in children and adolescents with depression.

Abstract Background: Androgen receptor (AR) is widely expressed in ER+ breast cancers. Accumulated evidence showed that AR is a tumor suppressor in ERα+ breast cancers. Enzalutamide (ENZ), an AR antagonist, blocks AR signalling, and was reported to inhibit estrogen-driven tumor growth as effectively as tamoxifen. Hence, this study was performed to further examine the mechanism of action of ENZ in ERα+ breast cancers. Methods: The AR gene was knocked out using the CRISPR-Cas9 system. Four crRNAs targeting AR and trans-activating crRNA were co-transfected to Cas9-positive MCF7 and T47D cells, two ERα+ breast cancer cell lines with relatively low AR levels based on RNA seq. ESR1 knock down was performed using a lentivirus system. MCF7 and T47D cells were infected with the ESR1-shRNA lentivirus and selected by puromycin for 2 weeks. Single clones of knockdown or knockout cells were tested using western blot. MCF7 and T47D cells were collected for RNA-seq and ChIP-seq after treated with DMSO or 10µM ENZ for 6 hours in complete medium. All sequences were aligned to the Human Reference Genome (hg38, UCSC). Differential gene expression was performed using limma package and upstream analysis was done in Ingenuity Pathway Analysis software. ChIP peaks were called using MACS2 and overlapping was analyzed using ChIPpeakAnno package in R. Tritium labeled ENZ or estradiol were incubated with ERα overnight. The bound radioactivity was measured using a Beckmann LS 6500 liquidscintillation counter. Results: We found the cell growth inhibition effect of ENZ was little impaired by AR gene knock out in MCF7 and T47D cells as was expected. As ENZ has been shown to inhibit estrogen-driven tumor growth, we performed studies to understand the mechanism of ENZ.We performed RNA-seq and identified 193 differential expressed genes (76 up and 117 down) between DMSO and ENZ-treated MCF7 cells (FC&gt;2 or &lt;0.5, FDR&lt;0.05). Upstream analysis showed that ESR1 might be the upstream regulator of those differentially expressed genes and ESR1 signaling was predicted to be inhibited after ENZ treatment (P=1.69×10-10). To examine the alteration of ERα genomic binding, we performed global ERα ChIP-seq and found 39.58% (8,681/21,933) of ERα-bound sites were abolished by ENZ and the binding intensity was significantly decreased upon ENZ treatment in 12,868 shared peaks. Moreover, radiolabled binding assays showed that ENZ displaced estrogen binding to ERα and directly bound to ERα with a KD= 399.60nM. We further observed a reduction of proliferation inhibition of ENZ in ESR1-KD cells compared with control cells. This suggests that ENZ inhibits ER+ cell growth by displacing E2 from ERα binding sites and subsequently diminishes ERα genomic binding and blocks ERα signaling in ER+ cancers. Conclusion: We found that enzalutamide inhibits cell growth not through AR but by antagonizing ERα activity in MCF7 and T47D ER+ cells. These results reveal ER is the main factor affecting enzalutamide efficacy in these two ERα+ breast cancer cell lines. Further clinical trials are needed to investigate AR antagonists in ERα+ breast cancers. Citation Format: Lixuan Wei, Jia Yu, Huanyao Gao, Huan Zhang, Thanh Nguyen, Yayun Gu, Richard M. Weinshilboum, James N. Ingle, Liewei Wang. Estrogen receptor is a target of enzalutamide in ER+ breast cancer [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr PD1-03.

Abstract Background: Cyclin dependent 4/6 kinase inhibitors (CDK4/6i) and endocrine therapy (ET) have improved progression-free survival (PFS) and overall survival in HR+ MBC, but progression of disease ultimately occurs. Apart from HR+ status, there are no clinically available biomarkers that enable oncologists to determine prognosis and predict response to CDK4/6i. An emerging biomarker is serum thymidine kinase 1 (TK1), a secreted marker of proliferation that is prognostic in pts with HR+ HER2- MBC. High levels of TKa are associated with inferior PFS, whereas pts with low TKa levels pretreatment, or TKa levels that decrease on ET and a CDK4/6i, have superior PFS. Notably, TKa levels rebound ≥ 5 days off Pb, with resumption of cell cycling. PROMISE (NCT0281902) is a prospective study that enrolled women with HR+ MBC starting Pb + letrozole (L) in the 1st line [FL] or Pb + fulvestrant in the 2nd line [SL] setting. The trial includes a comprehensive “omic” assessment of blood, tumor, urine and the fecal microbiome to identify novel genomic variants and pathways associated with an early decline in TKa (measured after 2 months or end of cycle [C]2) and PFS. Here, we report the association between i) pre-treatment TKa (pre-TKa) levels and PFS (i.e. from registration to the 1st disease event) and ii) TKa levels at the end of C2 (C2-TKa) and PFS-2 (i.e. from the start of C3 to the 1st disease event).Methods: TKa testing was performed using the DiviTum assay (Biovica). TKa+ disease was defined as ≥ 200 Du/L and TKa- disease as below limit of detection to 200 Du/L. Log-rank test and univariate Cox modeling were used to assess the association between pre-TKa levels and PFS and between end of C2-TKa levels and PFS-2. The database was locked on June 28, 2021. Results: Of 68 pts enrolled, 4 were ineligible and pre-TKa data was unavailable for 4. Of the remaining 60 pts (45 FL, 15 SL), the percentage of pts with pre-TKa+ disease was 33.3% in FL (15/45, 95% CI: 20.0-49.0%), and 46.7% (7/15, 95% CI: 21.4-71.9%) in the SL. The median follow-up time for pts on study was 24 months (range: 2-42 months). There were 22 disease events (13 in FL, 9 in SL). In the FL setting, PFS was significantly shorter for preTKa+ pts compared to preTKa- pts (HR: 4.15, 95% CI:1.35-12.74; p=0.007), but not for SL pts (HR: 1.11, 95% CI: 0.30-4.18, p=0.875). End of C2 TKa data was obtained for pts while on Pb (n=5), or after stopping Pb as follows: 1-4 days (n=9), 5-8 days (n=28) and 9-36 days (n=11). PFS-2 was not associated with C2-TKa in the FL (p=0.834) or SL (p=0.454) settings. An analysis of TKa levels by metastatic site will be presented at the meeting.Conclusions: A secreted biomarker of proliferation (TK1) obtained prior to initiating CDK4/6i and ET for the treatment of HR+ MBC is associated with PFS in pts receiving 1st line Pb + L, but not in those receiving 2nd line Pb + fulvestrant. While the end of C2 TKa levels were not associated with PFS, the interpretability of these data are limited, given treatment delays (0-36 days) prior to starting C3 that may result in TKa rebound. Future studies evaluating the predictive nature of TKa and Pb response should focus on earlier timepoints while on drug. Citation Format: Ciara C O'Sullivan, Jun He, Jason Sinnwell, Vera J Suman, Krishna R Kalari, Peter T Vedell, Ann M Moyer, Xiaojia Tang, Kevin J Thompson, Abe Eyman Casey, Alvaro Moreno-Aspitia, Donald W Northfelt, Minetta C Liu, Tufia C Haddad, Saranya Chumsri, Brendan McMenomy, Prema Peethambaram, Kathryn J Ruddy, Karthik V Giridhar, Roberto A Leon-Ferre, Mattias Bergqvist, Adrian Nordmark, Richard M Weinshilboum, Liewei Wang, Matthew P Goetz. Serum thymidine kinase 1 activity (TKa) levels and progression-free survival (PFS) in patients (pts) with hormone receptor positive (HR+) HER2-negative metastatic breast cancer (MBC) on palbociclib (Pb) and endocrine therapy (ET) [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P5-13-22.