Noninvasive prediction of axillary lymph node status in breast cancer using promoter profiling of circulating cell-free DNA
Journal of Translational Medicine volume 20, Article number: 557 (2022)
Lymph node metastasis (LNM) is one of the most important factors affecting the prognosis of breast cancer. The accurate evaluation of lymph node status is useful to predict the outcomes of patients and guide the choice of cancer treatment. However, there is still lack of a low-cost non-invasive method to assess the status of axillary lymph node (ALN). Gene expression signature has been used to assess lymph node metastasis status of breast cancer. In addition, nucleosome footprint of cell-free DNA (cfDNA) carries gene expression information of its original tissues, so it may be used to evaluate the axillary lymph node status in breast cancer.
In this study, we found that the cfDNA nucleosome footprints between the ALN-positive patients and ALN-negative patients showed different patterns by implementing whole-genome sequencing (WGS) to detect 15 ALN-positive and 15 ALN-negative patients. In order to further evaluate its potential for assessing ALN status, we developed a classifier with multiple machine learning models by using 330 WGS data of cfDNA from 162 ALN-positive and 168 ALN-negative samples to distinguish these two types of patients.
We found that the promoter profiling between the ALN-positive patients and ALN-negative patients showed distinct patterns. In addition, we observed 1071 genes with differential promoter coverage and their functions were closely related to tumorigenesis. We found that the predictive classifier based on promoter profiling with a support vector machine model, named PPCNM, produced the largest area under the curve of 0.897 (95% confidence interval 0.86–0.93).
These results indicate that promoter profiling can be used to distinguish ALN-positive patients from ALN-negative patients, which may be helpful to guide the choice of cancer treatment.
Lymph node metastasis (LNM) is one of the most important factors affecting the prognosis of breast cancer . The accurate assessment of lymph node status can predict patients’ outcomes and guide the choice of treatment options . Axillary lymph node dissection (ALND) is the gold standard for evaluating axillary lymph node (ALN) status, but it would bring great harm to patients. Although milder sentinel lymph node biopsy (SLNB) has become routine surgery, it is still risk surgery, which would increase considerable anesthesia time and expense, and cause multiple complications in 3.5–10.9% of patients [3, 4]. Therefore, developing a low-cost non-invasive method to evaluate the status of ALN would be of great benefit to breast cancer patients.
Cell-free DNA (cfDNA) has been an essential biomarker in many cancer applications, such as early detection and outcome prediction of cancer . At present, the most commonly used features are cfDNA level and its sequence information. Previous studies have described the close relationship between abnormal cfDNA levels and ALN metastasis [6, 7], which indicates that cfDNA may be used to assess ALN status. However, the level of cfDNA is influenced levels are affected by many pathological processes, such as infection and inflammation [8,9,10]. In addition, some studies wanted to find ALN metastasis-related ctDNA mutations or ctDNA hypermethylation [1, 11,12,13], however, no relationship was found between them . Thus, novel disease-specific features of cfDNA with high predictive efficacy are needed to be found for predicting LNM.
Recently, cfDNA coverage on gene promoter has found that it carried gene expression information of its original tissues [14, 15]. Plasma cfDNA is mainly released by apoptotic cells after enzymatic processing of chromatin . The DNA bound to the nucleosomes is retained, while the exposed DNA between the nucleosomes is digested. Analysis of cfDNA fragments derived from cancers showed that the promoter regions of active genes exhibited depleted coverage, which meant that nucleosome binding was less in these regions along with increased gene expression . In cancer patients, cfDNA is mainly derived from tumor and hematopoietic cells . More importantly, studies on breast cancer have shown that many gene expression signatures could be used to estimate the risk of distant relapses, and some of which have been commercialized, such as PAM50. In addition, the immune cells have been proved to play an important role in tumor metastasis, and the peripheral blood immunome of breast cancer patients is influenced by the existence and stage of cancer [17, 18]. Therefore, we assume that the cfDNA coverage at the gene promoter has potential to assess the ALN status.
In this study, we first compared the nucleosome footprint around the transcriptional start sites (TSS) of ALN-positive and ALN-negative breast cancer patients to identify genes with differential coverage. In order to further evaluate the potential of promoter profiling for evaluating ALN status, we developed a classifier for distinguishing ALN-positive and ALN-negative patients by using multiple machine learning models. Finally, we incorporated some clinicopathological characteristics in our classifier to test whether its performance would improve.
Participants and study design
From January 2018 to December 2019, before cancer therapy, plasma samples were prospectively collected from 330 breast cancer patients, including 162 ALN-positive and 168 ALN-negative patients. We excluded patients who: (1) were pregnant or lactating, (2) were metastatic breast cancer or had non-infiltrating tumors histologically, (3) had a hematopoietic system or inflammatory breast diseases, and (4) were ALN-negative patients diagnosed with fine needle aspiration biopsy. We reviewed all tumor specimens histopathologically and staged them according to the seventh edition of the American Joint Committee on Cancer (AJCC) staging system for breast cancer. All plasma samples were obtained under institutional review board of The First People's Hospital of Foshan approved protocols with written informed consent from all participants for research use (ID: L-7). Table 1 summarizes the characteristics of patients, including age, T stage, estrogen- (ER) and progesterone-receptor (PR) status, expression of human epidermal growth factor receptor 2 (Her2), proliferative fraction (Ki-67 labeling index), and histological grade.
The ALN status was ascertained clinically by fine needle aspiration biopsy, ALND or SLNB. Because the number of lymph nodes detected by the fine needle aspiration biopsy is limited, some positive lymph nodes may be ignored, which may increase the false positive rate of the evaluation model. Therefore, the patients with ALN-negative detected by fine needle aspiration biopsy were excluded from this study. Indocyanine green with a carbon nanoparticle suspension was used for SLNB and more than three LNs were checked for cancer.
Extracting and sequencing cfDNA
In total, 1 mL peripheral blood was collected using EDTA tubes from each patient and then immediately implemented two-step centrifugation to obtain the plasma. The centrifugation parameters were is 1600g for 10 min, followed by 10 min at 16,000g at 4 °C. Subsequently, the plasma was stored at − 80 °C before use. Each sample yielded at least 1 ng total cfDNA for sequencing. cfDNA was extracted from plasma by QIAamp DNA Blood Mini Kit (Qiagen). A starting amount of approximately 1–5 ng DNA was used for library construction with the Life Sciences Ion Xpress™ Plus Fragment Library Kit. The number of PCR cycles was set to 12. The DNA size distribution of libraries was analyzed on a Bioanalyzer instrument (Agilent Technologies, Singapore). Sequencing was performed with the Ion PI™ Hi-Q™ OT2 200 Kit and the Ion PI™ Hi-Q™ Sequencing 200 Kit on Ion Proton platform (ThermoFisher Scientific, USA) with 520 flow. The mean depth of the sequencing samples was approximately 0.3×.
Sequencing data processing
After sequencing, the raw read was aligned to the human reference genome (hg19) using bwa (ver.0.7.5). Then, SAMtools rmdup function (ver. 0.1.18) was used to remove the polymerase chain reaction duplicates . The GC-bias correction was implemented using the deeptools (ver.3.5.0) with the default setting. The calculation of tumor fraction and copy number-bias correlation were implemented using ichorCNA algorithm .
Promoter profiling calculation
The calculation of promoter profiling was similar to that used in our previous study [15, 21]. In briefly, gene information was downloaded from RefSeq of University of California Santa Cruz . The region ranging from − 1 KB to + 1 KB around the transcriptional start site of each transcript, was defined as the primary transcription start site (pTSS), was first identified. The read counts for each base at the pTSS were calculated using DANPOS with default setting . After read alignment, the read coverage at the pTSS was extracted from the aligned BAM files using bedtools (ver. 2.17.0). Then, the read coverage was normalized by the reads per kilobase per million mapped reads (RPKM)-like method. The normalized value of promoter profiling was calculated by the following formula:
here, the length of each transcript is equal to 2000 because of the pTSS region ranging from − 1 KB to + 1 KB around each transcriptional start site.
Models for evaluating lymph node status
To develop the evaluation classifier, the patients were firstly divided into three cohorts, including discovery, training and validation cohorts. In the discovery cohort, we identified the genes with differential promoter coverage. Then, the plasma samples were then divided into training and validation cohorts in a ratio of 7:3. Based on the training cohort data, we developed classifiers using three models, including support vector machine (SVM), logistic regression (LR), and linear discriminant analysis (LDA) models, to distinguish ALN-positive and ALN-negative tumors. The importance of the features was assessed with the sigFeature package of R. Then we selected top 100 features for further classifier construction. The SVM classifier was constructed with the linear kernel in e1071 package using the default setting. In order to identify the optimal gene combination with the largest area under the curve (AUC), backward method was adopted. To avoid potential bias and over-fitting in the training cohort, the leave-one-out cross validation method was used to evaluate the robustness of these classifiers. Briefly, each subject in the training cohort was withheld in turn, and the rest subjects were submitted to train the model. The trained model was then used to determine the class of the withheld subject. This procedure went on until all subjects in the training cohort were judged. Finally, the efficacy of selected classifiers was evaluated using the validation cohort data.
Wilcoxon rank-sum test or Chi square test were used for analyses that compared the two groups. Benjamini–Hochberg method was used to adjust the raw P-values to the false discovery rate (FDR). Variables with fold change ≥ 1.5 and FDR ≤ 0.05 were considered statistically significant. The genes with differential promoter coverage were used to plot uniform manifold approximation and projection (UMAP) and heat map using uwot package and pheatmap package in R (version 3.0.1), respectively. Receiver operating characteristic (ROC) curves were plotted and differences in the AUC were compared using the pROC package . GO enrichment analysis was implemented by using Metascape with default settings . Housekeeping genes and non-constitutive genes were downloaded from the additional material of a previous study .
cfDNA promoter profiling related to tumor expression profiles
In order to test whether the promoter profiling of cfDNA could be used to predict ALN metastasis, we first studied whether the coverage of gene promoter regions (± 1 KB around TSS) was related to gene expression profiles (Fig. 1). Consistent with previous studies , the promoter coverage of housekeeping genes with high expression levels was significantly reduced compared with those of non-constitutive genes (Fig. 3a). Then, we studied whether the footprint of nucleosomes around the TSS was different between ALN-positive and ALN-negative groups. In ALN-positive breast cancer patients, we observed the loss of related cfDNA signals (Fig. 3b; P = 2.2e−16, Wilcoxon rank sum test).
Genes with differential promoter coverage associated with LNM
The workflow of our study mainly consisted of three stages, including discovery, training and validation stages (Fig. 2). In the discovery cohort, we identified the genes with differential promoter coverage. When comparing the promoter profiling of each gene, we observed 1,071 genes with differential promoter coverage between ALN-positive and ALN-negative patients (Fig. 3e and Additional file 1: Table S1; fold change ≥ 1.5 and FDR ≤ 0.05, Wilcoxon rank sum test). Then, using UMAP, we found that samples from the same groups were clustered together, while the samples from different groups were scattered (Fig. 3c). In addition, the heat map showed distinct patterns of promoter coverage between ALN-positive and ALN-negative breast cancer patients (Fig. 3d). These results indicated that promoter profiling has potential for assessing the ALN status of breast cancer.
By GO enrichment analysis of the genes with differential promoter coverage, we found that most of GO terms were immuno-associated and growth-associated processes (Fig. 3f). Consistent with the existing literature, cfDNA could reflect the expression status of its original tissues. As the expression of tumor and peripheral blood immunome was closely related to cancer stage , the above annotation results may indicate that the genes with differential promoter coverage may be associated with ALN involvement.
Classifiers for evaluating ALN status
To evaluate the potential of promoter profiling for assessing ALN status, we used WGS to characterize the promoter profiling of cfDNA derived from 330 breast cancer patients collected from January 2018 to December 2019, including 162 ALN-positive and 168 ALN-negative patients. The patients were split into training and validation cohorts with a 7:3 ratio and the clinicopathological parameters, such as age, T stage, ER, PR, and Her2 status, were well balanced between the two cohorts of breast patients (Table 1; all P > 0.05).
Then, we used genes with differential promoter coverage in SVM model to develop classifiers to distinguish ALN-positive from ALN-negative patients. ROC analysis was used to evaluate the AUC, sensitivity, specificity and accuracy of the promoter profiling classifiers (Fig. 4a). Among these combinations, a 48-gene combination named PPCNM performed well in the training cohort after LOOCV, with an AUC of 0.936 (95% confidence interval [CI] 0.904–0.967 and an accuracy of 0.848, Fig. 4a and Additional file 1: Table S2). The performance of PPCNM was further evaluated in the validation cohort, and we found that the AUC of PPCNM in the validation cohort was 0.808 (0.730–0.887) (Fig. 4b). These results indicated that a classifier based on promoter profiling can be used to assess ALN status.
Across all cohorts, the average AUC of PPCNM was 0.897 (0.865–0.930), which was used to distinguish ALN-positive and ALN-negative patients, with a sensitivity of 0.914 and a specificity of 0.881 (Fig. 4c). The AUC produced by PPCNM was significantly greater than those of classifiers based on the LR and LDA models (Fig. 4c, LR: 0.829 [0.789–0.870], P = 8.37E−04 and LDA: 0.757 [0.711–0.803], P = 4.03E−10).
PPCNM and tumor DNA fraction
The level of tumor DNA fraction is one of the most important characteristics of tumor. Firstly, we calculated the tumor DNA fraction of ALN-positive and ALN-negative patients, and found that its levels between these two groups were similar (Additional file 1: Fig. S1; P-value = 0.1663). In addition, we found that the efficacy of PPCNM was similar in different concentrations of tumor DNA fraction (all P-value > 0.5; Additional file 1: Table S3). The AUC of tumor DNA fraction used to predict ALN status was 0.544 (0.482–0.606). The efficacy of the combination of PPCNM with tumor DNA fraction has an AUC 0.845 (0.806–0.885), which is significantly lower than that of PPCNM (P = 2.8E−04).
PPCNM combined with clinicopathological characteristics
Previous studies have shown the close relationship between ER, PR, Her2, and Ki67 status with ALN metastasis [26, 27]. Therefore, we first investigated whether the efficacy of our classifiers was different between positive- and negative-status of each feature. The efficacy of the PPCNM model was similar in regards to ER-positive vs. ER-negative, PR-positive vs. PR-negative, Her2-positive vs. Her2-negative, and Ki67-High and Ki67-Low (Fig. 5a–d). We then incorporated these clinical characteristics with PPCNM to see whether its performance would further improve. By evaluating the efficacy of their exhaustive combination with PPCNM, we found that the AUC, accuracy, sensitivity of the PPCNM decreased after being combined with one of the four clinical features (Fig. 5e, f and Additional file 1: Table S4).
We found that there was a significant difference in promoter profiling between ALN-positive and ALN-negative breast cancer patients (Fig. 3). The classifier PPCNM based on promoter profiling using the SVM model, produced the maximum AUC (0.897 [0.865–0.930]) for distinguishing these two groups of patients, and its performance was significantly better than those of classifiers relied on LR and LDA regression models (Fig. 4c; all P < 0.05). In addition, the AUC increased slightly with the incorporation of clinical characteristics. These findings indicate that PPCNM may be a promising non-invasive tool for evaluating ALN status.
There are forty-eight genes in the PPCNM (Additional file 1: Table S2). These genes are closely associated with the metastasis of tumor. For instance, a large number of studies have reported the close relationship between NF-κB signaling pathway and tumor metastasis [28, 29]. NF-κB signaling pathway regulates the expression of its downstream target genes, including MMP9, TNFα, uPA and IL8, thus promoting the invasion and metastasis of breast cancer cells . Besides, BHLHE40 confers a pro-survival and pro-metastatic phenotype to breast cancer cells by modulating HBEGF secretion . And BHLHE40 facilitates the invasion of cancer cell by interacting with SP1 . In addition, USP20 can promote breast cancer metastasis by stabilizing SNAI2 .
ALN status is an essential factor for the prognosis of breast cancer patients and the choice of cancer treatment in breast cancer . Although milder SLNB has become more pervasive, LN surgery for evaluating ALN status still brings various side effects to patients. Therefore, developing a non-invasive method to predict ALN status may be beneficial to breast cancer patients. At present, some studies show that increased cfDNA levels are related to ALN Metastasis [6, 7]. But cfDNA levels were affected by various physiological and pathological processes [8,9,10]. More specific features of cfDNA have to be found for assessing ALN status. Previous studies have reported that cell-free DNA promoter profiling and TF profiling is capable of prediction of tumor subtypes in prostate and detect early-stage colorectal cancer [14, 33]. Therefore, we assume that promoter profiling could be used to evaluate ALN status. In this study, we found the characteristics of specific promoter profile signatures of cfDNA in ALN-positive and ALN-negative patients (Fig. 3e). The classifiers (PPCNM) based on these differential variables achieved high performance with an AUC of 0.897 [0.865–0.930]. We developed a non-invasive method based on plasma cfDNA to assess ALN status, which could dynamically monitor the status of lymph node. More importantly, our method could avoid the heterogeneity of tumor in tissue detection. Nevertheless, there are some limitations in our research. Although the AUC of our classifier achieved 0.897, and 330 WGS data was used in this study, more prospective samples and samples from other external centers were needed to improve the predictive value of efficacy before clinical application.
In summary, our data suggest that PPCNM is a promising tool based on promoter profiling for evaluating ALN status in breast cancer. PPCNM is a non-invasive technique, which only needs low-coverage DNA sequencing and is not affected by cancer heterogeneity. Therefore, the PPCNM classifier may help patients and clinicians to choose appropriate cancer treatment methods, thus improving the curative effects and the quality of life of cancer.
Availability of data and materials
All datasets generated for this study are included in the article/additional material.
Lymph node metastasis
Axillary lymph node dissection
Sentinel lymph node biopsy
Whole genome sequencing
Transcriptional start site
Human epidermal growth factor receptor 2
Primary transcription start site
Support vector machine
Linear discriminant analysis
Leave one out cross-validation
Area under curve
Uniform manifold approximation and projection
Fiegl H, Millinger S, Mueller-Holzner E, Marth C, Ensinger C, Berger A, Klocker H, Goebel G, Widschwendter M. Circulating tumor-specific DNA: a marker for monitoring efficacy of adjuvant therapy in cancer patients. Cancer Res. 2005;65:1141–5.
Lee JH, Jeong H, Choi JW, Oh HE, Kim YS. Liquid biopsy prediction of axillary lymph node metastasis, cancer recurrence, and patient survival in breast cancer: a meta-analysis. Medicine. 2018;97: e12862.
Boughey JC, Moriarty JP, Degnim AC, Gregg MS, Egginton JS, Long KH. Cost modeling of preoperative axillary ultrasound and fine-needle aspiration to guide surgery for invasive breast cancer. Ann Surg Oncol. 2010;17:953–8.
Langer I, Guller U, Berclaz G, Koechli OR, Schaer G, Fehr MK, Hess T, Oertli D, Bronz L, Schnarwyler B, et al. Morbidity of sentinel lymph node biopsy (SLN) alone versus SLN and completion axillary lymph node dissection after breast cancer surgery: a prospective Swiss multicenter study on 659 patients. Ann Surg. 2007;245:452–61.
Cristiano S, Leal A, Phallen J, Fiksel J, Adleff V, Bruhm DC, Jensen SO, Medina JE, Hruban C, White JR, et al. Genome-wide cell-free DNA fragmentation in patients with cancer. Nature. 2019;570:385–9.
Agassi R, Czeiger D, Shaked G, Avriel A, Sheynin J, Lavrenkov K, Ariad S, Douvdevani A. Measurement of circulating cell-free DNA levels by a simple fluorescent test in patients with breast cancer. Am J Clin Pathol. 2015;143:18–24.
Payne RE, Hava NL, Page K, Blighe K, Ward B, Slade M, Brown J, Guttery DS, Zaidi SA, Stebbing J, et al. The presence of disseminated tumour cells in the bone marrow is inversely related to circulating free DNA in plasma in breast cancer dormancy. Br J Cancer. 2012;106:375–82.
Yi J, Zhang Y, Ma Y, Zhang C, Li Q, Liu B, Liu Z, Liu J, Zhang X, Zhuang R, Jin B. Increased plasma cell-free DNA level during HTNV infection: correlation with disease severity and virus load. Viruses. 2014;6:2723–34.
Tovbin D, Novack V, Wiessman MP, Abd Elkadir A, Zlotnik M, Douvdevani A. Circulating cell-free DNA in hemodialysis patients predicts mortality. Nephrol Dial Transplant. 2012;27:3929–35.
Tsai NW, Lin TK, Chen SD, Chang WN, Wang HC, Yang TM, Lin YJ, Jan CR, Huang CR, Liou CW, Lu CH. The value of serial plasma nuclear and mitochondrial DNA levels in patients with acute ischemic stroke. Clin Chim Acta. 2011;412:476–9.
Beaver JA, Jelovac D, Balukrishna S, Cochran R, Croessmann S, Zabransky DJ, Wong HY, Toro PV, Cidado J, Blair BG, et al. Detection of cancer DNA in plasma of patients with early-stage breast cancer. Clin Cancer Res. 2014;20:2643–50.
Mirza S, Sharma G, Parshad R, Srivastava A, Gupta SD, Ralhan R. Clinical significance of promoter hypermethylation of ERbeta and RARbeta2 in tumor and serum DNA in Indian breast cancer patients. Ann Surg Oncol. 2012;19:3107–15.
Gasch C, Oldopp T, Mauermann O, Gorges TM, Andreas A, Coith C, Muller V, Fehm T, Janni W, Pantel K, Riethdorf S. Frequent detection of PIK3CA mutations in single circulating tumor cells of patients suffering from HER2-negative metastatic breast cancer. Mol Oncol. 2016;10:1330–43.
Ulz P, Thallinger GG, Auer M, Graf R, Kashofer K, Jahn SW, Abete L, Pristauz G, Petru E, Geigl JB, et al. Inferring expressed genes by whole-genome sequencing of plasma DNA. Nat Genet. 2016;48:1273–8.
Guo Z, Yang F, Zhang J, Zhang Z, Li K, Tian Q, Hou H, Xu C, Lu Q, Ren Z, et al. Whole-genome promoter profiling of plasma DNA exhibits diagnostic value for placenta-origin pregnancy complications. Adv Sci. 2020;7:1901819.
Snyder MW, Kircher M, Hill AJ, Daza RM, Shendure J. Cell-free DNA comprises an in vivo nucleosome footprint that informs its tissues-of-origin. Cell. 2016;164:57–68.
Ali HR, Chlon L, Pharoah PD, Markowetz F, Caldas C. Patterns of immune infiltration in breast cancer and their clinical implications: a gene-expression-based retrospective study. PLoS Med. 2016;13: e1002194.
Foulds GA, Vadakekolathu J, Abdel-Fatah TMA, Nagarajan D, Reeder S, Johnson C, Hood S, Moseley PM, Chan SYT, Pockley AG, et al. Immune-phenotyping and transcriptomic profiling of peripheral blood mononuclear cells from patients with breast cancer: identification of a 3 gene signature which predicts relapse of triple negative breast cancer. Front Immunol. 2028;2018:9.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Genome Project Data Processing S. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Adalsteinsson VA, Ha G, Freeman SS, Choudhury AD, Stover DG, Parsons HA, Gydush G, Reed SC, Rotem D, Rhoades J, et al. Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors. Nat Commun. 2017;8:1324.
Guo ZW, Xiao WW, Yang XX, Yang X, Cai GX, Wang XJ, Han BW, Li K, Zhai XM, Li FX, et al. Noninvasive prediction of response to cancer therapy using promoter profiling of circulating cell-free DNA. Clin Transl Med. 2020;10: e174.
Casper J, Zweig AS, Villarreal C, Tyner C, Speir ML, Rosenbloom KR, Raney BJ, Lee CM, Lee BT, Karolchik D, et al. The UCSC genome browser database: 2018 update. Nucleic Acids Res. 2018;46:D762–9.
Chen K, Xi Y, Pan X, Li Z, Kaestner K, Tyler J, Dent S, He X, Li W. DANPOS: dynamic analysis of nucleosome position and occupancy by sequencing. Genome Res. 2013;23:341–51.
Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, Muller M. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:77.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, Benner C, Chanda SK. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.
Ahmed AR. HER2 expression is a strong independent predictor of nodal metastasis in breast cancer. J Egypt Natl Canc Inst. 2016;28:219–27.
Van Calster B, Vanden Bempt I, Drijkoningen M, Pochet N, Cheng J, Van Huffel S, Hendrickx W, Decock J, Huang HJ, Leunen K, et al. Axillary lymph node status of operable breast cancers by combined steroid receptor and HER-2 status: triple positive tumours are more likely lymph node positive. Breast Cancer Res Treat. 2009;113:181–7.
Shibata A, Nagaya T, Imai T, Funahashi H, Nakao A, Seo H. Inhibition of NF-kappaB activity decreases the VEGF mRNA expression in MDA-MB-231 breast cancer cells. Breast Cancer Res Treat. 2002;73:237–43.
Wang J, Li S, Li X, Li B, Li Y, Xia K, Yang Y, Aman S, Wang M, Wu H. Circadian protein BMAL1 promotes breast cancer cell invasion and metastasis by up-regulating matrix metalloproteinase9 expression. Cancer Cell Int. 2019;19:182.
Sethuraman A, Brown M, Krutilina R, Wu ZH, Seagroves TN, Pfeffer LM, Fan M. BHLHE40 confers a pro-survival and pro-metastatic phenotype to breast cancer cells by modulating HBEGF secretion. Breast Cancer Res. 2018;20:117.
Zheng Q, Wang C, Wang L, Zhang D, Liu N, Ming X, Zhou H, Guli Q, Liu Y. Interaction with SP1, but not binding to the E-box motifs, is responsible for BHLHE40/DEC1-induced transcriptional suppression of CLDN1 and cell invasion in MCF-7 cells. Mol Carcinog. 2018;57:1116–29.
Li W, Shen M, Jiang YZ, Zhang R, Zheng H, Wei Y, Shao ZM, Kang Y. Deubiquitinase USP20 promotes breast cancer metastasis by stabilizing SNAI2. Genes Dev. 2020;34:1310–5.
Ulz P, Perakis S, Zhou Q, Moser T, Belic J, Lazzeri I, Wolfler A, Zebisch A, Gerger A, Pristauz G, et al. Inference of transcription factor binding from cell-free DNA enables tumor subtype prediction and early detection. Nat Commun. 2019;10:4666.
We would like to thank Shun Xu for his valuable comments.
This work was supported by National Natural Science Foundation of China [81872416, 82173001, 81802435, 81900191]; Medical Scientific Research Foundation of Guangdong Province of China [B2017006].
Ethics approval and consent to participate
This study was approved by the Ethics Committee of The First People’s Hospital of Foshan (ID: L-7).
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The tumor DNA fraction between ALN-positive and ALN-negative patients. Table S1. Genes with differential promoter coverage between ALN-positive and -negative groups. Table S2. The genes in PPCNM. Table S3. The predictive efficacy of PPCNM in different tumor fraction. Table S4. Predictive efficacy of the PPCNM with characteristics.
About this article
Cite this article
Guo, ZW., Liu, Q., Yang, X. et al. Noninvasive prediction of axillary lymph node status in breast cancer using promoter profiling of circulating cell-free DNA. J Transl Med 20, 557 (2022). https://doi.org/10.1186/s12967-022-03724-w