Joohyuk Sohn, Min Hwan Kim, Jin Mo Ahn, Won-Ji Ryu, Seul-Gi Kim, Jee Hung Kim, Tae Yeong Kim, Hyun Ju Han, Jee Ye Kim, Hyung Seok Park, Seho Park, Byeong Woo Park, Seung Il Kim, Eun Hae Cho, and Gun Min Kim
Background Previous studies proposed low-pass whole genome sequencing (LP-WGS)-based circulating tumor DNA (ctDNA) analysis as a versatile tool for genomic profiling and therapeutic monitoring of cancer patients. Here we demonstrate LP-WGS ctDNA genomic profiles and its clinical significance in metastatic breast cancer patients. Patients and methods This prospective exploratory study enrolled 207 treatment-naïve metastatic breast cancer patients from Feb 2017 to September 2020 in Yonsei Cancer Center. The median follow-up duration of patients was 35 months. The baseline (n=207) and post-progression (n=48) plasma samples were prospectively collected on first-line systemic therapy, and LP-WGS was employed for ctDNA somatic copy number alteration (CNA) analysis. The CNA burden of ctDNA was scored by “I-score” method, which was developed to measure genome-wide chromosomal instabilities, to be matched with therapy response. The unsupervised molecular clustering and homologous recombination deficiency (HRD) estimation by shallowHRD algorithm were performed using locus-level CNA profiles with 1 mega base pair resolution. Results The baseline I-score ctDNA CNA burden was highest in triple-negative breast cancer (TNBC) patients among subtypes, and the patients were dichotomized by median I-score level 5.54 (range 2.55 to 12.98). The high baseline ctDNA I-score was independently associated with poor overall survival (hazard ratio [HR] = 3.98, p < 0.001) with adjustment of tumor subtype, visceral metastasis, and disease status (de novo stage IV versus recurrent). The progression-free survival (PFS) on endocrine plus CDK4/6 inhibitors (HR = 2.75, p = 0.005), anti-HER2 therapy (HR = 2.52, p = 0.032), and cytotoxic chemotherapy (HR = 2.33, p = 0.012) was also shorter in high baseline I-score patients than in low I-score patients. The locus-level CNA profile was analyzed in high I-score patients (n=103), and the patients were classified into five molecular clusters with distinct overall survival by unsupervised k-means clustering of CNA profile: basal-like, EGFR-high basal-like, CCND1-high, luminal, and HER2-enriched clusters. Patients with BCL6 (p = 0.009) and PIK3CA amplification (p < 0.001) on baseline ctDNA showed significantly shorter PFS on CDK4/6 inhibitor treatment. The matched baseline and post-progression ctDNA analysis found emergence of FGFR1 amplification and MYC amplification after CDK4/6 inhibitor treatment (n=1, each). The ctDNA shallowHRD score was highest in TNBC patients among subtypes, and TNBC patients with high shallowHRD score (≥10) showed high response rate on (58.3% versus 28.6%) on platinum-based chemotherapy. Conclusion LP WGS-based ctDNA analysis provides a robust tool for non-invasive genomic clustering, therapy response prediction, and HRD estimation in metastatic breast cancer patients. All patients (n=207)Low I-score (n=104)High I-score (n=103)N (%)N (%)N (%)Age, Median (Interquartile range)54 (46-62)53 (47-60)54(44-62)GenderFemale205 (99)102 (98.1)103Male2 (1)2 (1.9)0SubtypeHR+ HER2-106 (51.2)61 (58.7)45 (43.7)HR- HER2+33 (15.9)14 (13.5)19 (18.4)HR+ HER2+22 (10.6)11 (10.6)11 (10.7)HR- HER2- (TNBC)46 (22.2)18 (17.3)28 (27.2)Disease statusDe novo stage IV74 (35.7)31 (29.8)43 (41.7)Recurrent133 (64.3)73 (70.2)60 (58.3)Primary therapyEndocrine + CDK 4/6 inhibitor97 (46.9)55 (52.9)42 (40.8)Anti-HER2 based therapy54 (26.1)24 (23.1)30 (29.1)Chemotherapy45 (21.7)16 (15.4)29 (28.2)Others11 (5.3)9 (8.7)2 (1.9)Visceral metastasisYes142 (68.6)60 (57.7)82 (79.6)No65 (31.4)44 (42.3)21 (20.4)Metastasis SitesLung89 (43)43 (41.3)46 (44.7)Brain19 (9.2)4 (3.8)15 (14.6)Liver59 (28.5)13 (12.5)46 (44.7)Bone120 (58)47 (45.2)73 (70.9)Lymph node90 (43.7)32 (30.8)58 (56.9)Pleura33 (15.9)17 (16.3)16 (15.5) Citation Format: Joohyuk Sohn, Min Hwan Kim, Jin Mo Ahn, Won-Ji Ryu, Seul-Gi Kim, Jee Hung Kim, Tae Yeong Kim, Hyun Ju Han, Jee Ye Kim, Hyung Seok Park, Seho Park, Byeong Woo Park, Seung Il Kim, Eun Hae Cho, Gun Min Kim. Whole genome sequencing-based circulating tumor DNA profiling of metastatic breast cancer patients for molecular characterization and therapy response prediction [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 PD6-07.