Vedell, Peter T.

ABSTRACT Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has had an unprecedented impact on human health and highlights the need for genomic epidemiology studies to increase our understanding of virus evolution and spread and to inform policy decisions. We sequenced viral genomes from over 22,000 patient samples tested at Mayo Clinic Laboratories between 2020 and 2022 and used Bayesian phylodynamics to describe county and regional spread in Minnesota. The earliest calculated introduction into Minnesota was to Hennepin County from a domestic source around 22 January 2020; 6 weeks before the first confirmed case in the state. This led to the virus spreading to Northern Minnesota and, eventually, the rest of the state. International introductions were most abundant in Hennepin (home to the Minneapolis/St. Paul International Airport) totaling 45 (out of 107) over the 2-year period. Southern Minnesota counties were most common for domestic introductions, with 19 (out of 64), potentially driven by bordering states such as Iowa and Wisconsin as well as Illinois, which is nearby. Hennepin also was, by far, the most dominant source of in-state transmissions to other Minnesota locations ( n = 772) over the 2-year period. We also analyzed the diversity of the location source of SARS-CoV-2 viruses in each county and noted the timing of state-wide policies as well as trends in clinical cases. Neither the number of clinical cases nor the major policy decisions, such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, appeared to have an impact on virus diversity across each individual county. IMPORTANCE We analyzed over 22,000 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes of patient samples tested at Mayo Clinic Laboratories during a 2-year period in the COVID-19 pandemic, which included Alpha, Delta, and Omicron variants of concern to examine the roles and relationships of Minnesota virus transmission. We found that Hennepin County, the most populous county, drove the transmission of SARS-CoV-2 viruses in the state after including the formation of earlier clades including 20A, 20C, and 20G, as well as variants of concern Alpha and Delta. We also found that Hennepin County was the source for most of the county-to-county introductions after an initial predicted introduction with the virus in early 2020 from an international source, while other counties acted as transmission “sinks.” In addition, major policies, such as the end of the lockdown period in 2020 or the end of all restrictions in 2021, did not appear to have an impact on virus diversity across individual counties.

Abstract We identified resistance mechanisms to abiraterone acetate/prednisone (AA/P) in patients with metastatic castration-resistant prostate cancer (mCRPC) in the Prostate Cancer Medically Optimized Genome-Enhanced Therapy (PROMOTE) study. We analyzed whole-exome sequencing (WES) and RNA-sequencing data from 83 patients with metastatic biopsies before (V1) and after 12 weeks of AA/P treatment (V2). Resistance was determined by time to treatment change (TTTC). At V2, 18 and 11 of 58 patients had either short-term (median 3.6 months; range 1.4–4.5) or long-term (median 29 months; range 23.5–41.7) responses, respectively. Nonresponders had low expression of TGFBR3 and increased activation of the Wnt pathway, cell cycle, upregulation of AR variants, both pre- and posttreatment, with further deletion of AR inhibitor CDK11B posttreatment. Deletion of androgen processing genes, HSD17B11, CYP19A1 were observed in nonresponders posttreatment. Genes involved in cell cycle, DNA repair, Wnt-signaling, and Aurora kinase pathways were differentially expressed between the responder and non-responder at V2. Activation of Wnt signaling in nonresponder and deactivation of MYC or its target genes in responders was detected via SCN loss, somatic mutations, and transcriptomics. Upregulation of genes in the AURKA pathway are consistent with the activation of MYC regulated genes in nonresponders. Several genes in the AKT1 axis had increased mutation rate in nonresponders. We also found evidence of resistance via PDCD1 overexpression in responders. Implications: Finally, we identified candidates drugs to reverse AA/P resistance: topoisomerase inhibitors and drugs targeting the cell cycle via the MYC/AURKA/AURKB/TOP2A and/or PI3K_AKT_MTOR pathways.

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

Abstract 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-positive (HR+) metastatic breast cancer (MBC), but endocrine resistance is a major challenge. PROMISE [NCT0281902; n=63] is a multicenter study that enrolled women with HR+ HER2- MBC commencing palbociclib (Pb) with letrozole (1st line [1L]) or fulvestrant (2nd line [2L]), and was designed to perform a comprehensive “omic” assessment of prospectively collected biospecimens (pre-treatment (M1), at 2 months (M2), and at disease progression). The goal is to identify novel genomic variants and pathways associated with resistance to CDK4/6i and ET and PFS outcomes. Here we report the association between the proteomic, metabolomics, and lipidomics data generated from pre-Pb and 2-month serum samples and PFS. Methods: Untargeted mass spectrometry data was generated from Metabolon, assaying 1308 metabolites and 831 lipids. Additionally, 1436 proteins were assayed on the Olink platform. Cox proportional hazard models were used to evaluate the univariate hazard ratio (HR) for all features with respect to PFS. The analyses were performed on samples from 45 patients (N=33 1Lwith 9 progression events and 12 2L with 8 progression events), obtained from M1 and M2 timepoints on Pb + ET. Enrichment analysis p-values are calculated using Fisher’s exact test. Results: Proteomics: In the M1 timepoint, 93 and 43 proteins were associated with PFS in the 1L and 2L settings, respectively; inflammation genes were enriched among the 1L setting (p= 0.034); 33 proteins presented HRs ranging between 0.026 and 0.56. The FABP9 protein (HR of 1.98, 95% CI 1.02-3.83) was associated with worse PFS. Conversely, inflammation genes were not observed to be enriched in 2L. In the M2 timepoint, we observed 60 and 21 proteins significantly associated with PFS, but no biological function was enriched in 1L and 2L. Metabolites: In the M1 timepoint, metabolism of the sulfur-containing amino acids (methionine, cysteine, SAM and taurine) were enriched in the 1L setting (p= 0.035, HR range 0.15-0.33); and the branched-chain amino acids (leucine, isoleucine, and valine) were significantly associated with PFS in the 2L setting (p= 0.028, HR range 0.013-0.33). At the M2 timepoint, the amino acids were no longer enriched, but fatty acid metabolism was significantly enriched for both 1L and 2L (p= 0.048 and 0.067, respectively). Pathways involving lipids, amino acids, and xenobiotics were enriched in metabolites related to PFS (p &lt;0.05) for both treatment lines at M1 and M2. Lipidomics: In the M1 timepoint, 10 and 19 lipids were associated with PFS for 1L and 2L, respectively. The most notable lipid associated with worse PFS in the 1L was an 18 carbon phosphatidylinositol, PI(18:1/18:2), (HR 7.34 (CI 1.27-42.50); 8 triglycerides were associated with improved PFS (HR range 0.39 and 0.55). In 2L, the 19 lipids associated with PFS included 12 phosphatidylcholines (enrichment p = 5.6X10-8). In the M2 timepoint, 15 and 8 lipids were significantly associated with PFS for 1L and 2L. An enrichment of phosphatidylinositols was observed in 1L (p= 1.2X10-5); none were observed in the 2L.Future Directions: Networks are being constructed using the proximity scores of the proteins, lipids, and metabolites associated with PFS in M1 and M2 for 1L and 2L. Network similarities and analyses will be conducted.Conclusion: Distinct multi-omic changes identified in serum samples obtained from PROMISE participants M1 and M2 on Pb correlate with disease progression in both 1L and 2L settings. Additionally, validation studies will determine the significance of these findings. Citation Format: Krishna R. Kalari, Kevin J. Thompson, Jason Sinnwell, Xiaojia Tang, Vera J. Suman, Jun He, Seul Kee Byeon, Akhilesh Pandey, Abe Eyman Casey, Peter T. Vedell, Ann M. Moyer, 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, Richard M. Weinshilboum, Liewei Wang, Ciara C. O’ Sullivan, Matthew P. Goetz. Multiomics data reveal novel biomarkers for CDK4/6 resistance [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 P4-01-03.

Abstract Introduction. The standard treatment for non-metastatic triple-negative breast cancer (TNBC) is neoadjuvant chemotherapy (NAC) and nearly 50% exhibit pathological complete response (pCR). However, patients with residual disease after NAC are at increased risk for recurrence and death. Prior studies examining the transcriptome of TNBC pre/post-NAC have examined a limited number of genes (&lt;500) in heterogeneous subgroups of TNBC (e.g. LAR and non-LAR). We explored the transcriptome of androgen-receptor (AR) negative (non-LAR) TNBC subtype both pre/post NAC to identify pathways associated with NAC response. Methods. Tumors obtained pre/post NAC from TNBC patients enrolled in the Breast Cancer Genome Guided therapy study (BEAUTY) underwent RNA sequencing and reverse-phase protein array (RPPA). EdgeR was applied for differentially expressed (DE) analysis and regression methods for RPPA. Digital deconvolution method (CIBERSORTx) and TNBC single-cell data were used to obtain cell types. Pathway analysis was carried out using 2972 gene sets and gene set variation analysis (GSVA). Functional enrichment analysis was conducted with significant genes. Results. Of the 44 TNBC patients, 32 patients were excluded from the analysis cohort due to: LAR tumor (6 pts.), non-LAR tumor with pCR (23 pts.), and cell type issues with RNA-seq data (3 pt.). Paired RNA-Seq data were available for 12 TNBC patients (4 with progression &lt;2 years [EP]) and 8 who were progression-free &gt; 4 years [NP]) and paired RPPA data were available for 9 of these 12 patients. Differentially expressed genes, proteins and cell types between EP and NP in post-NAC. We identified 489 genes differentially expressed (DE) between EP and NP (logFC=|2|, FDR &lt; 0.05). Analysis of cytobands from these 489 genes showed an enrichment of genes on chromosome 6p22.1-2 and 17q25.3 regions (enrichment ratio &gt;5; p-value &lt;10E-4). Critical genes identified in the AR- network (p-value &lt; 10E-3) were IL1RN, SLAMF9, KRT81, BHLHE22, B3GALT5, PCP4, TREM1, AQP9, NRTN, and COL2A1.In addition, preliminary results from RPPA data of post-NAC tumors showed astrocytic phosphoprotein (PEA-15), involved in apoptosis, proliferation, glucose metabolism, as well as cell proliferation and Y box binding (YB1) proteins (involved in metastases), were more DE in EP than NP (p &lt; 0.05). CIBERSORTx was applied to estimate the proportions of different cell types in post-NAC tumors. Cancer-associated fibroblasts iCAFs were low and myCAFs are high in EP vs NP. It is known that the cross-talk between CAFs and tumor cells may induce tumor resistance to chemotherapy. Differentially expressed pathways in post and pre-NAC EP tumors. Using genome-wide expression data from the paired 12 tumors and the GSVA method, we obtained individual pathway scores for 2972 pathways. One hundred ninety pathways were downregulated and 61 pathways were upregulated (p-value &lt;= 0.05) in the post-NAC residual disease of EP relative to NP. We further examined these 190 pathways in the paired EPs and found 71% of those pathways were upregulated in the pre-NAC. These 190 downregulated pathways were enriched with FOXO, TGF-beta, PI3k, FGFR1, insulin and others. The 61 upregulated pathways in post-NAC EP tumors were enriched with mismatch repair, purine, tubulin, telomere, polymerase and gap-junction related pathways; 77% of those 61 pathways were downregulated in pre-NAC. Conclusions. Using a comprehensive “omics” approach, we have identified novel cancer and drug response pathways associated with recurrence in AR-TNBC disease. Further work to evaluate these as markers of outcome and potential drug targets is warranted. Citation Format: Krishna R Kalari, Vera J Suman, Xiaojia Tang, Jason P Sinnwell, Kevin J Thompson, Peter T Vedell, Jodi M Carter, Sarah A McLaughlin, Alvaro Moreno Aspitia, Donald W Northfelt, Richard J Gray, Richard Weinshilboum, Liewei Wang, Judy C Boughey, Matthew Goetz. Multi-omics data shows downregulation of mismatch repair, purin and tublin pathways in AR-negative triple-negative chemotherapy-resistant tumors [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 P4-01-05.

Abstract Background: FGFR dysregulation is observed in multiple cancers and targeting FGFR is an emerging therapeutic strategy with FDA approved treatments in bladder and cholangiocarcinoma. Here we examined the prevalence of FGFR mutations, fusions, and high-level amplifications in breast cancer, stratified by receptor subtype and local/metastatic status, in both Foundation Medicine (FM) and institutional Mayo Clinic (MC) cohorts. Methods: For the FM cohort, comprehensive genomic profiling (CGP) examining at least 324 genes for all classes of alterations, including FGFR1-4 was carried out for 32,048 breast cancers during the course of routine clinical care in a Clinical Laboratory Improvement Amendments (CLIA)-certified lab (Foundation Medicine Inc., Cambridge, MA, USA). Tumor mutational burden (TMB) was determined on up to 1.1 Mb, microsatellite instability high (MSI-High) was determined on up to 114 loci and predicted ancestry from &gt;10,000 SNPs. Estrogen receptor (ER) and HER2 status were available for a subset of FM samples. Additionally, 131 patients with metastatic breast cancer from a subset of patients at three Mayo Clinic sites (MC cohort) with clinical characteristics and cancer-panel DNA sequencing data from a CLIA-certified lab (Tempus, Chicago, IL) were included. Results: In the FM cohort, the prevalence of FGFR1-4 high-level amplification (CN≥10) was 10.1%, while mutations (1.5%) and fusions (0.72%) were rare. Most amplifications occurred in FGFR1 (9.2%); most fusions and mutations occurred in FGFR2 (0.46%, 0.77%). FGFR alteration prevalence was highest in ER+/HER2- subtype (14.4%) and lowest in HER2+ disease (7.7%). FGFR alterations were more common in IDC (11.7%) than ILC (7.7%), p&lt;3E-08. FGFR alterations were more prevalent in the metastatic setting relative to breast-biopsied disease (13.6% v 10.1%; OR = 1.4; p=2E-17), especially in the HER2+ (OR =1.9, p=0.004) and ER-/HER2- (OR = 1.9, p = 0.05) disease; no enrichment was seen in the ER+/HER2- metastases (OR =1.0, p = 1). FGFR amplifications were observed at a higher prevalence in patients with predicted East Asian ancestry, relative to patients with European ancestry (12.1% v 10.0%; p = 0.03). Overall, the most common activating mutations in FGFR were FGFR2 N549K (n=85), FGFR1 N546K (n=78), FGFR4 V510M (n=28), FGFR2 K659E (n=28), FGFR4 V510L (n=20), and FGFR2 Y375C (n=15). The most common recurrent fusions were FGFR3:TACC3 (n=36), FGFR2:TACC2 (n=17), FGFR1:TACC1 (n=9), FGFR1:BAG4 (n=6), and FGFR2:ATE1 (n=5). In patients with FGFR amplifications, the most frequently co-occurring alterations were ZNF703 (78.4%), TP53 (51.5%), CCND1 (36.1%), FGF3/4/19 (32.9 - 34.4%), PIK3CA (30.7%), MYC (29.6%), ESR1 (17.2%), EMSY (16.3%), and PTEN (10.6%). Significant co-occurrence was observed for a number of genes including FGF3/4/19, CDK4, and CDK8 (all OR&gt;2, p&lt;1E-07); mutual exclusivity was observed with PIK3R1, BRCA1, and BRCA2 (all OR &lt;0.5, p&lt;4E-13), among other genes. In the 131 metastatic tumors from MC, the prevalence of FGFR1-4 high-level amplifications was 19.8% [FGFR1 (12.4%), FGFR2 (7.4%), and FGFR3 (0.8%)]. The prevalence of high-level FGFR amplifications did not differ by clinical subtypes: HR-/HER2- (7/31), HR+/HER2- (15/79), and HER2+ (2/11), p=0.68. Conclusions: High-level FGFR amplifications are observed in &gt;11% of breast cancers, especially the ER+/HER2- subtype, while mutations/fusions are rare. These data support the ongoing studies evaluating targeted therapies for FGFR amplified ER + breast cancer. Correlations with clinical information (MC cohort) and associations with actionable alterations are ongoing and may inform potential combination strategies. Citation Format: Karthik V Giridhar, Ethan S Sokol, Peter T Vedell, Jason P Sinnwell, Aakash Desai, Tufia C Haddad, Ciara C O’Sullivan, Roberto A Leon-Ferre, Siddhartha Yadav, Kostandinos Sideras, Brenda Ernst, Minetta C Liu, Abe Eyman Casey, Xiaojia Tang, Zoe Fleischmann, Karthikeyan Murugesan, Krishna R Kalari, Matthew P Goetz. The frequency and somatic mutation landscape of Fibroblast growth factor receptor (FGFR) alterations in 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 P3-08-02.

Abstract 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 &gt; 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.

Abstract Background: With few targeted therapies available for triple-negative breast cancer (TNBC), there is a critical need to identify biomarkers that can predict chemotherapy-response. Copy number amplifications (CNAmps) have been implicated in tumorigenesis, progression, and treatment response. Previously, we performed a prospective neoadjuvant chemotherapy study (BEAUTY) and identified a 5 Mb CNAmp in TNBC patients associated with pathological complete response rates to taxane- and anthracycline-based chemotherapy. Here we further interrogated this amplicon and its role in chemotherapy response. Methods: Using siRNA screening of genes identified in the 5 Mb CNAmp, followed by cytotoxicity assays, we identified RDH10 significantly associated with chemo-sensitivity in TNBC cell lines. Functional studies were performed using RDH10 knockdown or overexpression followed by cytotoxicity, cell proliferation, luciferase reporter assays, HPLC or western blot. Results: RDH10 is a rate-limiting step involved in the synthesis of all-trans retinoic acid (ATRA) - the active metabolite of retinoid metabolism that is frequently aberrant in cancer. Previous studies demonstrated ATRA to be a potent inhibitor of PIN1, a key regulator of several oncogenic signaling pathways often amplified in TNBC. We discovered that genetic manipulation of RDH10 affected ATRA intracellular concentrations, contributing to PIN1 degradation and taxane-response. Furthermore, co-treatment of TNBC cells with ATRA significantly increased taxane-sensitivity. Conclusion: RDH10 amplification increases endogenous levels of ATRA and degrades PIN1 in TNBC. ATRA-mediated degradation of PIN1 sensitizes TNBC to chemotherapy, suggesting that RDH10 may be a novel biomarker of taxane response and the addition of ATRA to standard chemotherapy may improve chemo-sensitivity.