Skip to main content
Download PDF

Epigenetic quantification of circulating immune cells in peripheral blood of triple-negative breast cancer patients

Clinical Epigenetics volume 13, Article number: 207 (2021) Cite this article

Abstract

Background

A shift in the proportions of blood immune cells is a hallmark of cancer development. Here, we investigated whether methylation-derived immune cell type ratios and methylation-derived neutrophil-to-lymphocyte ratios (mdNLRs) are associated with triple-negative breast cancer (TNBC).

Methods

Leukocyte subtype-specific unmethylated/methylated CpG sites were selected, and methylation levels at these sites were used as proxies for immune cell type proportions and mdNLR estimation in 231 TNBC cases and 231 age-matched controls. Data were validated using the Houseman deconvolution method. Additionally, the natural killer (NK) cell ratio was measured in a prospective sample set of 146 TNBC cases and 146 age-matched controls.

Results

The mdNLRs were higher in TNBC cases compared with controls and associated with TNBC (odds ratio (OR) range (2.66–4.29), all Padj. < 1e−04). A higher neutrophil ratio and lower ratios of NK cells, CD4 + T cells, CD8 + T cells, monocytes, and B cells were associated with TNBC. The strongest association was observed with decreased NK cell ratio (OR range (1.28–1.42), all Padj. < 1e−04). The NK cell ratio was also significantly lower in pre-diagnostic samples of TNBC cases compared with controls (P = 0.019).

Conclusion

This immunomethylomic study shows that a shift in the ratios/proportions of leukocyte subtypes is associated with TNBC, with decreased NK cell showing the strongest association. These findings improve our knowledge of the role of the immune system in TNBC and point to the possibility of using NK cell level as a non-invasive molecular marker for TNBC risk assessment, early detection, and prevention.

Introduction

Inflammation plays an important role in almost every stage of cancer development. Many inflammatory markers have been associated with cancer progression and prognosis [1]. Various studies showed that the number and function of blood leukocytes are altered in cancer [2,3,4,5]. A shift in the number of peripheral immune cells is a predictor of cancer patient survival. For instance, an increased neutrophil-to-lymphocyte ratio (NLR) which is indicative of systemic inflammation, could promote cancer cell proliferation, angiogenesis, cellular migration, and metastasis [6]. There is evidence from many studies for a prognostic role of NLR in peripheral blood of various cancer patients [7,8,9,10], including breast cancer [8, 11, 12].

Triple-negative breast cancer (TNBC) accounts for 15–20% of all breast cancers [13]. Due to the lack of targeted therapies, chemotherapy still is the main therapeutic strategy. Therefore, many efforts have been conducted to increase the diagnostic and therapeutic opportunities for the TNBC patients [14]. TNBC is also the most immunogenic subtype [15, 16]. Higher levels of infiltrated T cells are associated with an improved OS and disease-free survival (DFS) of TNBC patients as compared with those affected by other breast cancer subtypes [15, 16]. In another recent study, four spatially distinct tumor immune microenvironment subtypes defined by distinct CD8 + T cell localization patterns and gene expression signatures were identified, which were associated with distinct disease outcomes [17]. Further, NLR is associated with survival in a pre-treatment setting and throughout the treatment course and subsequent follow-up [5, 18, 19]. In addition, higher peripheral lymphocyte counts are associated with a lower mortality from early-stage TNBC suggesting that immune cell functions improve early TNBC treatment [20].

Epigenetic modifications such as DNA methylation play an important role in the cell-specific gene regulation within the hematopoietic system [21, 22]. Since DNA methylation signatures are chemically stable and mitotically heritable, they have been successfully applied to quantify leukocyte subtypes accurately in DNA from peripheral blood [23,24,25].

In the present study, we identified and validated associations of methylation-derived leukocyte subtype ratios and methylation-derived neutrophil-to-lymphocyte ratios (mdNLRs) with TNBC using methylation data of 231 TNBC cases and 231 age-matched controls from a retrospective study. We report associations of various leukocyte subtypes ratios with TNBC, with the natural killer (NK) cell ratio showing the strongest association with disease. Further, we provide evidence for an association of the NK cell ratio with TNBC risk in a prospective sample set of 146 TNBC cases and 146 age-matched controls (for a graphical overview of the work, see Fig. 1).

Fig. 1
figure 1

Graphical representation of the analyses performed in the study

Full size image

Methods

Study populations

In the present study, two TNBC case–control sample sets were analyzed: One set from a retrospective study contains 231 TNBC cases and 231 age-matched controls selected from two different studies: the breast cancer case-only study SKKDKFZS [26] and the breast cancer case–control study GENICA [27]. The prospective sample set includes 146 TNBC cases and 146 controls from the Generations Study (GS) [28]. All study participants were women of European ethnicity. All studies had local ethical approval and all included individuals gave informed written consent. Details on the study populations can be found in Additional file: 9 (available online). The sample sizes and selected characteristics of study participants are presented in Additional file 1, 2: Tables S1 and S2.

DNA methylation analysis

We applied publically available DNA methylation datasets obtained from isolated leukocyte subtypes (NK cells, CD4 + T cells (TCD4 +), CD8 + T cells (TCD8 +), regulatory T cells (Treg), monocytes, neutrophils, and B cells) for selection of cell type-specific differentially methylated sites. Thereafter, DNA methylation levels (beta values) at selected CpG sites were used as proxies for leukocyte subtypes proportions and mdNLRs. The mdNLR was estimated by dividing the beta value at neutrophil-specific CpGs by the beta value of pan-lymphocyte-specific CpGs. Details on the experimental strategies and data analysis are described in Additional file 9 (available online).

In-house genome-wide and locus-specific DNA methylation analysis was performed on the retrospective and prospective sample sets, respectively. Genome-wide DNA methylation profiling was performed on 231 TNBC cases and 231 controls from a retrospective study using the Illumina Infinium HumanMethylation450K BeadChip according to the manufacturer´s instructions. Subsequently, for measurement of the NK cell level in the prospective cohort of 146 TNBC cases and 146 controls, MethyLight droplet digital PCR (ddPCR) assay was carried out on bisulfite converted DNA. More details on the materials and methods are described in Additional file 9.

For technical validation of our findings, we applied a reference-based method (Houseman algorithm) to statistically deconvolute the proportions of six immune cell subtypes (NK, TCD4 + , TCD8 + cells, monocytes, neutrophils, and B cells) in the TNBC cases and controls [29]. The mdNLRreference-based method (mdNLRref) was calculated by dividing the estimated proportion of neutrophils by the sum of the lymphoid cell proportions (NK, TCD4 + , TCD8 + , and B cells).

Statistical analyses

Associations of TNBC and leukocyte subtype ratios were analyzed using the methylation beta values from the retrospective study. The ratios of the seven leukocyte subtypes in TNBC cases and controls were calculated and compared using methylation levels (beta value) of 21 selected immune cell-specific unmethylated sites (ISUS). A higher level of methylation at each ISUS corresponds to lower ratio of the corresponding leukocyte subtype.

In primary analysis, the diagnostic performance of individual CpGs and subtype ratios on TNBC status was assessed with univariable conditional logistic regression accounting for age-matched pairs (age matching ± 1 year) and receiver operating characteristic (ROC)/area under the curve (AUC) analysis. In supportive analysis, adjusted effects were estimated based on a multivariable conditional logistic regression model additionally accounting for body mass index (BMI) (continuous), menopausal status (pre/peri, post), ever parous (yes, no), and smoking status (current) (yes, no). Complete data for 221 matched pairs were available for the supportive analysis.

A multivariable logistic regression model based on all pre-selected CpGs was fitted with backward selection at a significance level of 20% for staying in the model. Internal validation of the AUC for the multivariable model including variable selection was done using bootstrapping with 200 repetitions. To assess the diagnostic performance of subtype ratios, differences in logit-transformed beta values were analyzed. In TNBC cases, the association of methylation levels with clinical and epidemiological parameters (age, menopausal status, BMI, smoking status, ever parous and number of children), histopathological tumor parameters (grade, size, node status, stage) and overall survival (OS) was assessed and in controls associations with age, menopausal status, BMI, and smoking status. OS was defined as the time between TNBC diagnosis and death or last follow-up, whichever occurred first. Mann–Whitney test, Jonckheere–Terpstra trend test, and Spearman’s correlation coefficient were used to assess associations between methylation levels and clinical, epidemiological, and histopathological parameters. The impact of beta values on OS was analyzed in a Cox regression model. To account for established prognostic factors, a multivariable Cox regression model including age, tumor grade (G1/G2 vs G3), stage (0–4), tumor size (T1, T2, T3, T4) and N status (N0 vs N1) was fitted. Kaplan–Meier estimates and log-rank test were derived for methylation levels at median cutoff. Individual P values were adjusted for multiple testing using Holm correction to control the family-wise error rate. The obtained “unmethylation “ values from ddPCR experiment were logit-transformed and univariable conditional logistic regression accounting for age-matched and follow-up time matched pairs (age matching ± 5 years) as well as Wilcoxon signed-rank test were used to compare the ratios between cases and controls. All analyses have been done using R 3.6 with add-on packages rms, survival and pROC.

Results

Associations of mdNLRs with TNBC

The methylation levels of the selected neutrophil- and pan-lymphocyte-specific CpG sites are shown in Fig. 2A. The neutrophil-specific sites are methylated in the neutrophils and unmethylated in pan-lymphocytes (NK, TCD4 + , TCD8 + , and B cells) and monocytes. The pan-lymphocyte-specific sites are methylated in the pan-lymphocytes and unmethylated in neutrophils and monocytes. The characteristics of the CpG sites are provided in Table 1.

Fig. 2
figure 2

Heat map visualization of DNA methylation levels at the three pan-lymphocyte- and the three neutrophil-specific CpG sites using GEO datasets (A). Box plots of data on the methylation levels of the pan-lymphocyte- and neutrophil-specific CpG sites in 231 TNBC cases and 231 controls from the retrospective study (B). Outliers appear outside of the whiskers

Full size image
Table 1 Characteristics of selected CpG proxies in neutrophils and pan-lymphocytes
Full size table

Methylation analysis at the selected CpG sites in 231 TNBC cases and 231 controls revealed a significantly higher ratio of neutrophils and lower ratio of pan-lymphocytes in TNBC cases compared with controls (mean ratios: 49% vs. 45% and 36% vs. 40.3%, respectively; conditional logistic regression: all Padj. < 1e−04) (Fig. 2B). The neutrophil ratios were associated with a higher likelihood of being a TNBC case (OR range (2.07–3.02); conditional logistic regression: all Padj. < 1e−04) and the pan-lymphocyte ratios with a lower likelihood (OR range (0.56–0.61); conditional logistic regression: all Padj. < 1e−04) (Additional file 3: Table S3). In cases, there were no associations of neutrophil and pan-lymphocyte ratios with epidemiological parameters, histopathological characteristics, and OS (size/stage: Jonckheere–Terpstra trend test; ever parous: Mann–Whitney test, number of children: Spearman’s correlation coefficient; Grade/Node status: Mann–Whitney test; OS: Cox regression; all P > 0.05). The average means of the mdNLR for cases and controls were 1.5 ± 0.63 and 1.18 ± 0.43, respectively. Logistic regression analysis revealed that all mdNLRs were associated with an increased likelihood of being a TNBC case (OR range (2.66–4.29); all Padj. < 1e−04) (Fig. 3). No associations were observed for any of the nine mdNLRs with epidemiological characteristics, histopathological characteristics, and OS.

Fig. 3
figure 3

Forest plot of the associations of the nine mdNLRs with TNBC and estimated areas under the receiver operating characteristic curves (AUCs) with their corresponding 95% confidence intervals (CIs). OR is given for a 10% increase in methylation level. P values were adjusted (Padj) for multiple testing using the Holm correction. Horizontal lines indicate 95% CIs. mdNLR, methylation-derived neutrophil-to-lymphocyte ratio; OR, odds ratio

Full size image

To validate our mdNLR findings, we used the reference-based method (mdNLRref), which estimates cell proportions. A higher mdNLRref was also observed in TNBC cases compared with controls (2.49 ± 1.53 versus 1.81 ± 0.94) (Additional file 8: Fig. S1), which was associated with a higher likelihood of being a TNBC case (OR = 1.71 [1.38–2.11], conditional logistic regression; Padj. < 1e−04) (Additional file 8: Fig. S2). The mdNLR and mdNLRsref were highly correlated (Additional file 8: Fig. S3), with the mdNLR cg23954655.cg26942829 ratio showing the highest correlation (Spearman’s rank correlation; r = 0.97).

Adjusting additionally for confounding factors (BMI, menopausal status, ever parous, smoking status) in a multivariable conditional logistic regression model, the associations of neutrophil and pan-lymphocyte ratios, mdNLRs, and mdNLRref with TNBC did not change fundamentally (Additional file 3–5: Tables S3–S5).

Associations of leukocyte subtype ratios with TNBC

The methylation levels of the 21 selected ISUS are shown in Fig. 4A. All 21 ISUS were hypomethylated in the target cell types but methylated in the other subtypes. The characteristics of each ISUS are provided in Table 2. In order to investigate the association between leukocyte subtypes ratios with TNBC in the retrospective sample set, the immune cell type ratios of TNBC cases were compared with those in controls using beta values of the specific ISUS. Analysis of these proxies showed that six of seven immune cell type ratios were associated with TNBC. Lower ratios of NK, TCD4 + , TCD8 + cells, monocytes, and B cells in cases compared with controls were associated with TNBC (Fig. 4B), with decreased NK cell ratios showing the strongest association. Further, a higher ratio of neutrophils in TNBC cases compared with controls was associated with TNBC, while no difference between the two groups was found in the Treg cell ratio.

Fig. 4
figure 4

Heat map illustrating the methylation levels at the 21 selected ISUS in target immune cells versus other immune cells (A). Forest plot showing the associations of the leukocyte subtype ratios and TNBC using the ISUS methylation data and estimated areas under the receiver operating characteristic curves (AUCs) with their corresponding 95% confidence intervals (CIs) (B). P values were adjusted (Padj) for multiple testing using the Holm correction. Horizontal lines indicate 95% CIs. ISUS, immune cell-specific unmethylated site; OR, odds ratio

Full size image
Table 2 Characteristics of selected ISUS as proxies and their participation in the diagnostic models
Full size table

To validate our findings obtained with the immune cell methylation proxies, we applied immune cell proportions estimated by reference-based deconvolution method. For each sample, the estimated proportions of the six immune cell types sum to one. Univariable comparison of the various immune cell types between TNBC cases and controls revealed a statistically significant difference in the proportions of neutrophils, NK, TCD4 + , and B cells (Additional file 8: Fig. S4). No difference was observed in the proportions of TCD8 + cells and monocytes. Logistic regression analysis showed associations of neutrophils, NK, TCD4 + , and B cell proportions with TNBC, with a decreased NK cell proportion showing the strongest association signal (Additional file 8: Fig. S2).

After adjustment for confounding factors, the associations of the leukocyte subtype ratios and proportions remained statistically significant except the association with cg07499259, which was no longer statistically significant (Additional file 5, 6: Tables S5 and S6).

Correlation of leukocyte subtype ratios with clinical, epidemiological, and histopathological parameters of the participants of the retrospective study

Correlation analysis of the immune cell type ratios with selected clinical and epidemiological characteristics of the study participants revealed correlations with smoking status and age. The neutrophil and B cell ratios correlated with smoking status in controls. Current smokers had a lower neutrophil and a higher B cell ratio compared with non-smokers (Mann–Whitney test; Padj < 0.05). Two immune cell type ratios correlated with age. The TCD8 + cell ratio showed an inverse correlation with age in both cases and controls and the NK cell ratio a positive correlation in controls (Spearman’s rank correlation; Padj < 0.05). No other correlations were observed. In TNBC cases, there were also no correlations of the immune cell type ratios with ever parous, number of children, and histopathological tumor characteristics (grade, size, node status, stage).

Diagnostic and prognostic performance of ISUS

AUC analysis showed that NK cells and neutrophils had the highest discriminative capability among all immune cell types (Fig. 4B). The estimated NK cell-to-neutrophil ratio was higher in controls compared with cases and slightly improved the discrimination performance between cases and controls with AUC values in the range (0.67–0.71) (Additional file 8: Fig. S5) relative to the values from individual CpGs in the range (0.63–0.67) (Fig. 4B). The NK cell-to-neutrophil ratios were associated with TNBC: a higher NK cell-to-neutrophil ratio was associated with a lower likelihood of being a TNBC cases (OR range (0.52-0-70); conditional logistic regression; all Padj. < 1e−04) (Additional file 7: Table S7). Next, a diagnostic model was developed by fitting a multivariable logistic regression model based on all 21 ISUS and applying backward variable selection. A bootstrap-adjusted AUC was computed to account for overfitting. The final model contained four ISUS that discriminated cases from controls with an AUC of 72%. Of these four ISUS, two were specific for NK cells, one for monocytes, and one for TCD8 + cells (Table 3).

Table 3 Multivariable logistic regression after variable selection
Full size table

The prognostic performance using log-rank test showed that two probes, cg00219921 and cg08326410, which are specific for TCD8 + and NK cells, were associated with survival when using a median split. Higher ratios of TCD8 + and NK cells were associated with a better patient OS (Additional file 8: Fig. S6). Using a multivariable Cox regression model including age, tumor grade, stage, tumor size, and lymph node status, only the association with cg00219921 remained statistically significant (Cox regression; P = 0.04). However, after adjustment for multiple testing, the association lost statistical significance.

Association of the NK cell ratio with TNBC in participants of the prospective case–control study

Since NK cells were the most pronounced immune cell type associated with TNBC in the retrospective sample set, we investigated whether the observed association could be detected in pre-diagnostic DNA samples of TNBC cases compared with controls. In this respect, the NK cell ratio was measured in 146 TNBC cases and 146 controls using a ddPCR TaqMan assay specific for one NK cell-specific unmethylated site (cg23060465). The obtained “unmethylation” levels from ddPCR were Logit-transformed and the Wilcoxon signed-rank test was used to compare the obtained NK cell ratios between cases and controls. A lower NK cell ratio was observed in TNBC cases compared with controls (Wilcoxon signed-rank test, P = 0.019) (Fig. 5). Using conditional logistic regression, a higher NK cell ratio was associated with a reduced TNBC risk at the margin of statistical significance (OR = 0.76, 95% CI [0.58–1.00], conditional logistic regression; P = 0.052). Heterogeneity due to age at blood draw, age of diagnosis, and interval time between blood draw and reference date was tested. No heterogeneity/subgroup effect was observed (Interaction test based on conditional logistic regression; all P > 0.05).

Fig. 5
figure 5

Scatter plot of the DNA “unmethylation” ratio in prediagnosed TNBC cases and controls. The unmethylation level corresponds to the NK cell ratio in sample groups

Full size image

Discussion

In the present large immunomethylomic study, we showed that methylation-derived leukocyte subtype ratios and mdNLRs at subtype-specific CpG sites are associated with TNBC in peripheral blood of TNBC patients and age-matched controls from a retrospective study. We showed that mdNLRs were higher in TNBC cases compared with controls and associated with TNBC. Further, higher ratios of neutrophils and lower ratios of NK, TCD4 + , TCD8 + , monocytes, and B cells were associated with TNBC, with a decreased NK cell ratio showing the strongest association. Associations of mdNLRs and neutrophil, NK, TCD4 + , and B cell ratios were validated in analysis based on immune cell type proportions. Moreover, we confirmed that the NK cell ratio was significantly lower in pre-diagnostic samples of TNBC cases compared with controls.

Overall, a 4% higher ratio of neutrophils and a 4.3% lower ratio of total lymphocytes were observed in TNBC cases compared with controls. It was reported that breast tumors maximize their chance of metastasizing by evoking a systemic inflammatory cascade, which leads to an elevated level of neutrophils. These tumor-induced neutrophils suppress cytotoxic TCD8 + lymphocytes, which ultimately enhancing metastatic seeding in the pre-metastatic lung [30]. Neutrophils also infiltrate many other tumor types, and the tumor microenvironment controls neutrophil recruitment [31]. They may act as potent antitumor effector cells (N1 neutrophils) via direct cytotoxicity or by the activation of different innate and adaptive immune cells. Alternatively, neutrophils also may acquire a protumor activity (N2 neutrophils). This is now supported by a growing number of studies showing that neutrophils correlate with poor cancer prognosis [32]. Unfortunately, due to the lack of N1- and N2-specific methylation datasets of neutrophils, the estimation of these subtypes of neutrophils was not feasible in our sample cohort.

Furthermore, the interactions between neutrophils and lymphocytes play critical roles in carcinogenesis. As a hallmark of cancer, in this study also an elevated mdNLR was observed in cases compared with controls (1.5 vs. 1.18). With slight difference, this result was repeated also by Houseman estimation method (averaged mdNLRref for cases and controls were 2.49 and 1.81, respectively). This finding agrees with that from a previous study, which reported mdNLR values of 2.7 in TNBC cases and 2.4 in controls [3]. A high NLR was associated with adverse survival of patients affected by various solid tumors, including TNBC [5, 7, 9, 18]. However, in contrast to the previous studies, mdNLR was not an independent predictor of OS in our study. The discrepant results obtained in these studies may be explained by differences in the study size, study population, and the methods used for mdNLR estimation.

Among the major immune cell types, higher neutrophil and lower ratios of B cells, TCD4 + , TCD8 + , NK cells, and monocytes were associated with TNBC. NK cell ratio showed the strongest association in the large retrospective sample set of 231 TNBC cases and 231 controls. Since the blood samples were drawn after the diagnosis of TNBC, it cannot be excluded that the shift in the NK cell ratio was induced by the tumor. There is evidence from a previous study that the tumor itself could manipulate and decrease the numbers of immune cells, such as NK cells, in the peripheral blood of TNBC patients by secreting cytokines [33]. However, a lower NK cell ratio was also observed in the prospective sample set of 146 TNBC cases compared with controls, which was associated with a higher TNBC risk at the margin of statistical significance. This finding suggests that the shift in NK cell ratio has occurred before disease manifestation. Interestingly, a lower NK cell level was also reported in head and neck cancer patients compared with controls with individuals in the lowest NK tertile having over fivefold risk of being a case [2].

There is evidence that NK cell levels are linked with survival of cancer patients. One study demonstrated that the NK cell ratio in the peripheral blood is an independent predictor of survival of colorectal cancer patients; those with a higher percentage of NK cells being associated with a better survival than those with a lower percentage [10]. Higher peripheral NK cell counts were also associated with better OS in lymphoma and chronic lymphocytic leukemia patients [34,35,36]. Further, there is some evidence that the number of NK cells in peripheral blood may affect the outcome of B cell non-Hodgkin lymphoma patients receiving immunochemotherapy [37].

NK cells are effector lymphocytes of the innate immune system that control several types of tumors and microbial infections by limiting their spread [38]. They have a crucial role in the control of metastasis by virtue of killing circulating tumor cells and acting as the first line of defense against metastasis from circulating tumor cells [39]. Therefore, the lower numbers of circulating NK cells in patients with TNBC may be an underlying cause of the disease, and individuals with decreased NK cell levels could be predisposed to develop TNBC. One previous study reported that mice with deficiencies in NK cell number and function are more susceptible to transplanted tumors [40]. In addition, it has been shown before that individuals with low cytotoxic activity of peripheral blood lymphocytes, including NK cells, are at higher risk of developing various types of cancer [41].

Impairment of NK cell function was also reported to play a role in breast tumorigenesis [42]. Breast tumors modify their environment to evade NK cell antitumor immunity [43] and ex vivo-expanded NK cells showed potent antitumor function against breast cancer cell lines and primary cells isolated from patients [44, 45]. Another study suggested that NK cells are important players in TNBC development and metastasis that could be used as a promising immunotherapeutic against TNBC [46]. In this respect, increasing NK cell antitumor activity and ex vivo expanding NK cell populations is paving the way for a new generation of anticancer immunotherapies [4]. Our findings may highlight the value of NK cell-based immunotherapies for TNBC; given that NK cell-induced lysis was significantly higher in TNBC cell lines compared to estrogen receptor positive breast cancer cell lines [47]. In addition to NK and T cells, peripheral blood also contains NKT cells, which are considered at the interface between innate and adaptive immunity. In contrast to NK cells, which are CD3-negative and CD56-positive, NKT cells express both CD3 and CD56 and were independently associated with cancer survival [48]. To the best of our knowledge, the genome-wide methylation profiling of NKT cells is not available at the present time. Given the lack of NKT cell-specific DNA methylation markers to be used as proxy for this immune cell type, an estimation of the NKT cells in our sample set was not possible. At present, flow cytometry is the most widely applied analytical approach for immune cell quantification [25]. This method, however, is limited to intact cells, but fresh or well-preserved blood samples are not available for many clinical cohorts in which samples were often collected in numerous centers, different geographic regions, and within a time period of several years. Therefore, in various recent studies, including this one, an epigenetic assay has been used [2, 3, 23, 25], which allows cell quantification in samples of limited quality and quantity, is applicable to archival blood or DNA samples, which are available for many clinical cohorts, and is less costly than flow cytometry. On the other hand, applying genome-wide DNA methylation analysis to estimate the immune cell ratios by data deconvolution demands high costs and expertise, which may not be available for some clinical studies. It has been shown that the epigenetic assay performs equivalently to flow cytometry for immune cell quantification [23,24,25].

As is the case for all studies, this work is not without some limitations. One limitation is that epigenetic data could not be validated using classical cell counting by flow cytometry due to the lack of fresh blood samples from TNBC cases and controls. Another limitation is the still relatively limited set of only main immune cell subtypes that were investigated. Another limitation is that although the developed diagnostic model with four selected methylation markers improved the capability to distinguish patients with TNBC from controls to an AUC of 72% (compared with AUC values of individual CpGs in the range of 55% to 67%), this has no clinical utility based on a threshold AUC of > 80%. A higher discriminative capability has to be achieved by integrating other molecular markers of different sources such as cell-free nucleic acids or proteins in future models of non-invasive diagnosis of TNBC. Additionally, single cell expression and DNA methylation profiling of peripheral immune cells, especially NK cells, in TNBC patients and healthy controls may reveal unique NK cell states or signatures in individuals developing TNBC.

Conclusion

In summary, this is the largest study investigating immune cell profiles in TNBC patients and controls using methylation data. We identified and validated associations of mdNLRs and neutrophil, TCD4 + , B, and NK cell subtype ratios/proportions with TNBC, with the latter having the strongest association signal. The NK cell ratio was also significantly lower in pre-diagnostic samples of TNBC cases compared with controls. Future studies on the peripheral NK cell population in TNBC patients may provide NK cell-based cellular/molecular signatures that may be useful as a potential non-invasive blood-based biomarker for TNBC risk assessment, early detection or immunotherapeutic applications.

Availability of data and materials

The datasets of sorted immune cells analyzed during the current study are available in the Gene Expression Omnibus repository (GEO; https://www.ncbi.nlm.nih.gov/geo/) through GEO accessions: GSE35969, GSE59250, GSE110554, GSE88824, and GSE49667. The in-house generated methylation array data are not publically available. They can be accessed for research purposes upon request.

Abbreviations

TNBC:

Triple-negative breast cancer

NK cell:

Natural killer cell

mdNLR:

Methylation-derived neutrophil-to-lymphocyte ratio

OR:

Odds ratio

ISUS:

Immune cell-specific unmethylated sites

ROC:

Receiver operating characteristic

AUC:

Area under the curve

OS:

Overall survival

ddPCR:

Droplet digital PCR

mdNLRref :

Methylation-derived neutrophil-to-lymphocyte ratio estimated by Housman method

References

  1. Grivennikov SI, Greten FR, Karin M. Immunity, inflammation, and cancer. Cell. 2010;140(6):883–99.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Accomando WP, Wiencke JK, Houseman EA, Butler RA, Zheng S, Nelson HH, et al. Decreased NK cells in patients with head and neck cancer determined in archival DNA. Clin Cancer Res. 2012;18(22):6147–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Koestler DC, Usset J, Christensen BC, Marsit CJ, Karagas MR, Kelsey KT, et al. DNA methylation-derived neutrophil-to-lymphocyte ratio: an epigenetic tool to explore cancer inflammation and outcomes. Cancer Epidemiol Biomark Prev. 2017;26(3):328–38.

    Article  CAS  Google Scholar 

  4. Chiossone L, Dumas PY, Vienne M, Vivier E. Natural killer cells and other innate lymphoid cells in cancer. Nat Rev Immunol. 2018;18(11):671–88.

    Article  CAS  PubMed  Google Scholar 

  5. Patel DA, Xi J, Luo J, Hassan B, Thomas S, Ma CX, et al. Neutrophil-to-lymphocyte ratio as a predictor of survival in patients with triple-negative breast cancer. Breast Cancer Res Treat. 2019;174(2):443–52.

    Article  PubMed  Google Scholar 

  6. Giraldo NA, Becht E, Vano Y, Sautes-Fridman C, Fridman WH. The immune response in cancer: from immunology to pathology to immunotherapy. Virchows Arch. 2015;467(2):127–35.

    Article  CAS  PubMed  Google Scholar 

  7. Ren K, Yin Y, He F, Shao Y, Wang S. Prognostic role of derived neutrophil-to-lymphocyte ratio in surgical triple-negative breast cancer. Cancer Manag Res. 2018;10:4891–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Ethier JL, Desautels D, Templeton A, Shah PS, Amir E. Prognostic role of neutrophil-to-lymphocyte ratio in breast cancer: a systematic review and meta-analysis. Breast Cancer Res. 2017;19(1):2.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Templeton AJ, McNamara MG, Seruga B, Vera-Badillo FE, Aneja P, Ocana A, et al. Prognostic role of neutrophil-to-lymphocyte ratio in solid tumors: a systematic review and meta-analysis. J Natl Cancer Inst. 2014;106(6):dju124.

    Article  PubMed  Google Scholar 

  10. Tang YP, Xie MZ, Li KZ, Li JL, Cai ZM, Hu BL. Prognostic value of peripheral blood natural killer cells in colorectal cancer. BMC Gastroenterol. 2020;20(1):31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Azab B, Bhatt VR, Phookan J, Murukutla S, Kohn N, Terjanian T, et al. Usefulness of the neutrophil-to-lymphocyte ratio in predicting short- and long-term mortality in breast cancer patients. Ann Surg Oncol. 2012;19(1):217–24.

    Article  PubMed  Google Scholar 

  12. de la Pena FA, Soler JCY, Ivars MA, Sanchez AF, Manzano EN, Barrio PD, et al. Neutrophil-lymphocyte ratio (NLR) as a prognostic factor in metastatic breast cancer. Ann Oncol. 2017;28:94.

    Article  Google Scholar 

  13. Metzger-Filho O, Tutt A, de Azambuja E, Saini KS, Viale G, Loi S, et al. Dissecting the heterogeneity of triple-negative breast cancer. J Clin Oncol. 2012;30(15):1879–87.

    Article  CAS  PubMed  Google Scholar 

  14. Garrido-Castro AC, Lin NU, Polyak K. Insights into molecular classifications of triple-negative breast cancer: improving patient selection for treatment. Cancer Discov. 2019;9(2):176–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Liu Z, Li M, Jiang Z, Wang X. A comprehensive immunologic portrait of triple-negative breast cancer. Transl Oncol. 2018;11(2):311–29.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Disis ML, Stanton SE. Triple-negative breast cancer: immune modulation as the new treatment paradigm. Am Soc Clin Oncol Educ Book. 2015;35:25–30.

    Article  Google Scholar 

  17. Gruosso T, Gigoux M, Manem VSK, Bertos N, Zuo D, Perlitch I, et al. Spatially distinct tumor immune microenvironments stratify triple-negative breast cancers. J Clin Investig. 2019;129(4):1785–800.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Pistelli M, De Lisa M, Ballatore Z, Caramanti M, Pagliacci A, Battelli N, et al. Pre-treatment neutrophil to lymphocyte ratio may be a useful tool in predicting survival in early triple negative breast cancer patients. BMC Cancer. 2015;15:195.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Walker PR, Hildebrand JR. Pretreatment neutrophil to lymphocyte ratio (NLR) in patients with triple negative breast cancer (TNBC) treated with neoadjuvant chemotherapy. J Clin Oncol. 2016;34(15):12570.

    Article  Google Scholar 

  20. Afghahi A, Purington N, Han SS, Desai M, Pierson E, Mathur MB, et al. Higher absolute lymphocyte counts predict lower mortality from early-stage triple-negative breast cancer. Clin Cancer Res. 2018;24(12):2851–8.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Bocker MT, Hellwig I, Breiling A, Eckstein V, Ho AD, Lyko F. Genome-wide promoter DNA methylation dynamics of human hematopoietic progenitor cells during differentiation and aging. Blood. 2011;117(19):e182–9.

    Article  CAS  PubMed  Google Scholar 

  22. Calvanese V, Fernandez AF, Urdinguio RG, Suarez-Alvarez B, Mangas C, Perez-Garcia V, et al. A promoter DNA demethylation landscape of human hematopoietic differentiation. Nucleic Acids Res. 2012;40(1):116–31.

    Article  CAS  PubMed  Google Scholar 

  23. Kelsey KT, Wiencke JK. Immunomethylomics: a novel cancer risk prediction tool. Ann Am Thorac Soc. 2018;15(Suppl 2):S76–80.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Accomando WP, Wiencke JK, Houseman EA, Nelson HH, Kelsey KT. Quantitative reconstruction of leukocyte subsets using DNA methylation. Genome Biol. 2014;15(3):R50.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Baron U, Werner J, Schildknecht K, Schulze JJ, Mulu A, Liebert UG, et al. Epigenetic immune cell counting in human blood samples for immunodiagnostics. Sci Transl Med. 2018;10(452).

  26. Stevens KN, Fredericksen Z, Vachon CM, Wang X, Margolin S, Lindblom A, et al. 19p13.1 is a triple-negative-specific breast cancer susceptibility locus. Cancer Res. 2012;72(7):1795–803.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Pesch B, Ko Y, Brauch H, Hamann U, Harth V, Rabstein S, et al. Factors modifying the association between hormone-replacement therapy and breast cancer risk. Eur J Epidemiol. 2005;20(8):699–711.

    Article  PubMed  Google Scholar 

  28. Swerdlow AJ, Jones ME, Schoemaker MJ, Hemming J, Thomas D, Williamson J, et al. The breakthrough generations study: design of a long-term UK cohort study to investigate breast cancer aetiology. Br J Cancer. 2011;105(7):911–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinform. 2012;13:86.

    Article  Google Scholar 

  30. Coffelt SB, Kersten K, Doornebal CW, Weiden J, Vrijland K, Hau CS, et al. IL-17-producing gammadelta T cells and neutrophils conspire to promote breast cancer metastasis. Nature. 2015;522(7556):345–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Uribe-Querol E, Rosales C. Neutrophils in cancer: two sides of the same coin. J Immunol Res. 2015;2015:983698.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Masucci MT, Minopoli M, Carriero MV. Tumor associated neutrophils their role in tumorigenesis, metastasis. Progn Ther Front Oncol. 2019;9:1146.

    Article  Google Scholar 

  33. Lee JC, Lee KM, Kim DW, Heo DS. Elevated TGF-beta1 secretion and down-modulation of NKG2D underlies impaired NK cytotoxicity in cancer patients. J Immunol. 2004;172(12):7335–40.

    Article  CAS  PubMed  Google Scholar 

  34. He L, Zhu HY, Qin SC, Li Y, Miao Y, Liang JH, et al. Low natural killer (NK) cell counts in peripheral blood adversely affect clinical outcome of patients with follicular lymphoma. Blood Cancer J. 2016;6(8):e457.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Wang WT, Zhu HY, Wu YJ, Xia Y, Wu JZ, Wu W, et al. Elevated absolute NK cell counts in peripheral blood predict good prognosis in chronic lymphocytic leukemia. J Cancer Res Clin Oncol. 2018;144(3):449–57.

    Article  PubMed  Google Scholar 

  36. Plonquet A, Haioun C, Jais JP, Debard AL, Salles G, Bene MC, et al. Peripheral blood natural killer cell count is associated with clinical outcome in patients with aaIPI 2–3 diffuse large B-cell lymphoma. Ann Oncol. 2007;18(7):1209–15.

    Article  CAS  PubMed  Google Scholar 

  37. Klanova M, Oestergaard MZ, Trneny M, Hiddemann W, Marcus R, Sehn LH, et al. Prognostic impact of natural killer cell count in follicular lymphoma and diffuse large B-cell lymphoma patients treated with immunochemotherapy. Clin Cancer Res. 2019;25(15):4634–43.

    Article  CAS  PubMed  Google Scholar 

  38. Vivier E, Tomasello E, Baratin M, Walzer T, Ugolini S. Functions of natural killer cells. Nat Immunol. 2008;9(5):503–10.

    Article  CAS  PubMed  Google Scholar 

  39. Larsen SK, Gao Y, Basse PH. NK cells in the tumor microenvironment. Crit Rev Oncog. 2014;19(1–2):91–105.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Guillerey C, Smyth MJ. NK cells and cancer immunoediting. Curr Top Microbiol Immunol. 2016;395:115–45.

    CAS  PubMed  Google Scholar 

  41. Imai K, Matsuyama S, Miyake S, Suga K, Nakachi K. Natural cytotoxic activity of peripheral-blood lymphocytes and cancer incidence: an 11-year follow-up study of a general population. Lancet. 2000;356(9244):1795–9.

    Article  CAS  PubMed  Google Scholar 

  42. Tsavaris N, Kosmas C, Vadiaka M, Kanelopoulos P, Boulamatsis D. Immune changes in patients with advanced breast cancer undergoing chemotherapy with taxanes. Br J Cancer. 2002;87(1):21–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Mamessier E, Sylvain A, Thibult ML, Houvenaeghel G, Jacquemier J, Castellano R, et al. Human breast cancer cells enhance self tolerance by promoting evasion from NK cell antitumor immunity. J Clin Invest. 2011;121(9):3609–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Nham T, Poznanski SM, Fan IY, Vahedi F, Shenouda MM, Lee AJ, et al. Ex vivo-expanded natural killer cells derived from long-term cryopreserved cord blood are cytotoxic against primary breast cancer cells. J Immunother. 2018;41(2):64–72.

    Article  CAS  PubMed  Google Scholar 

  45. Shenouda MM, Gillgrass A, Nham T, Hogg R, Lee AJ, Chew MV, et al. Ex vivo expanded natural killer cells from breast cancer patients and healthy donors are highly cytotoxic against breast cancer cell lines and patient-derived tumours. Breast Cancer Res. 2017;19(1):76.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Dewan MZ, Terunuma H, Takada M, Tanaka Y, Abe H, Sata T, et al. Role of natural killer cells in hormone-independent rapid tumor formation and spontaneous metastasis of breast cancer cells in vivo. Breast Cancer Res Treat. 2007;104(3):267–75.

    Article  CAS  PubMed  Google Scholar 

  47. Engel JB, Honig A, Kapp M, Hahne JC, Meyer SR, Dietl J, et al. Mechanisms of tumor immune escape in triple-negative breast cancers (TNBC) with and without mutated BRCA 1. Arch Gynecol Obstet. 2014;289(1):141–7.

    Article  CAS  PubMed  Google Scholar 

  48. Krijgsman D, de Vries NL, Skovbo A, Andersen MN, Swets M, Bastiaannet E, et al. Characterization of circulating T-, NK-, and NKT cell subsets in patients with colorectal cancer: the peripheral blood immune cell profile. Cancer Immunol Immunother CeII. 2019;68(6):1011–24.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors thank all study participants, clinicians, family doctors, researchers and technicians for their contributions and commitment to this study. We thank J Hoheisel for critical reading of the manuscript. The GENICA Network: Molecular Genetics of Breast Cancer, German Cancer Research Center (DKFZ), Heidelberg, Germany (UH); Dr. Margarete Fischer-Bosch Institute of Clinical Pharmacology, Stuttgart, and University of Tübingen, Tübingen, Germany (HB, W-Y Lo, R Hoppe, S Winter); German Cancer Consortium (DKTK) and German Cancer Research Center (DKFZ), Partner Site Tübingen (HB), funded by the Deutsche Forschungsgemeinschaft (DFG) – EXC 2180-390900677 (HB); Department of Internal Medicine, Evangelische Kliniken Bonn gGmbH, Johanniter Krankenhaus, Bonn, Germany (Y-DK, C Baisch); Institute of Pathology, University of Bonn, Germany (H-P Fischer); Institute for Prevention and Occupational Medicine of the German Social Accident Insurance, Institute of the Ruhr University Bochum (IPA), Bochum, Germany (TB, B Pesch, S Rabstein, A Lotz); and Institute of Occupational Medicine and Maritime Medicine, University Medical Center Hamburg-Eppendorf, Germany (V Harth). O Fletcher’s work is supported by Programme Grants from Breast Cancer Now as part of Programme Funding to the Breast Cancer Now Toby Robins Research Centre. The Generation Study thanks Breast Cancer Now and the Institute of Cancer Research for funding and support. The Institute of Cancer Research acknowledges NHS funding to the Royal Marsden Foundation Trust and the National Institute for the Health Research (NIHR) Biomedical Research Centre.

Funding

Open Access funding enabled and organized by Projekt DEAL. The work was funded by the German Cancer Research Center (DKFZ), Heidelberg.

Author information

Authors and Affiliations

  1. Molecular Genetics of Breast Cancer, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 580, 69120, Heidelberg, Germany

    Mehdi Manoochehri, Nasim Borhani & Ute Hamann

  2. Department of in-Vitro Diagnostics, Fraunhofer Institute for Interfacial Engineering and Biotechnology IGB, Stuttgart, Germany

    Mehdi Manoochehri

  3. Division of Biostatistics, German Cancer Research Center (DKFZ), Im Neuenheimer Feld 280, 69120, Heidelberg, Germany

    Thomas Hielscher

  4. Cancer Epigenomics, German Cancer Research Center (DKFZ), Heidelberg, Germany

    Clarissa Gerhäuser

  5. The Breast Cancer Now Toby Robins Research Centre, The Institute of Cancer Research, London, UK

    Olivia Fletcher

  6. The Institute of Cancer Research, London, UK

    Anthony J. Swerdlow

  7. Division of Genetics and Epidemiology and Division of Breast Cancer Research, The Institute of Cancer Research, London, UK

    Anthony J. Swerdlow

  8. Department of Internal Medicine, Evangelische Kliniken Bonn gGmbH, Johanniter Krankenhaus, 53113, Bonn, Germany

    Yon-Dschun Ko

  9. Dr. Margarete Fischer-Bosch Institute of Clinical Pharmacology, 70376, Stuttgart, Germany

    Hiltrud Brauch

  10. iFIT Cluster of Excellence, University of Tübingen, 72074, Tübingen, Germany

    Hiltrud Brauch

  11. German Cancer Consortium (DKTK) and German Cancer Research Center (DKFZ), Partner Site Tübingen, 72074, Tübingen, Germany

    Hiltrud Brauch

  12. Institute for Prevention and Occupational Medicine of the German Social Accident Insurance, Institute of the Ruhr University Bochum (IPA), 44789, Bochum, Germany

    Thomas Brüning

Authors
  1. Mehdi Manoochehri
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Thomas Hielscher
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Nasim Borhani
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Clarissa Gerhäuser
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Olivia Fletcher
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Anthony J. Swerdlow
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Yon-Dschun Ko
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Hiltrud Brauch
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Thomas Brüning
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Ute Hamann
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

MM and UH participated in the conception and design of the study. OF, AJS, Y-DK, HB, TB provided the study material or patients. CG, TH, NB, MM, and UH carried out the analysis and interpretation of data. MM, TH, CG, and UH wrote the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Mehdi Manoochehri or Ute Hamann.

Ethics declarations

Ethics approval and consent to participate

All procedures performed in studies involving human were in accordance with the ethical standards of the institutional research committee and with the 1964 Helsinki declaration. Informed consent was obtained for all study participants.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1

. Selected characteristics of the TNBC cases and controls from the retrospective study and breast tumor parameters

Additional file 2: Table S2

. Selected characteristics of the TNBC cases and controls from the prospective cohort

Additional file 3: Table S3

. Associations of the neutrophil and pan-lymphocyte ratios with TNBC after adjustment for multiple testing and confounders

Additional file 4: Table S4

. Associations of the nine mdNLRs with TNBC after adjustment for confounders

Additional file 5: Table S5

. Associations of the immune cell type proportions and mdNLRref with TNBC after adjustment for confounders

Additional file 6: Table S6

. Associations of the leukocyte subtype ratios and TNBC after adjustment for confounders

Additional file 7: Table S7

. Associations of NK cell-to-neutrophil ratios with TNBC after adjustment for multiple testing and confounders

Additional file 8: Fig. S1

. Distribution of the mdNLRref between TNBC cases and controls using the reference-based Houseman method. Fig. S2. Forest plot of the associations of the immune cell subtype proportions and mdNLRref with TNBC. Odds ratio (OR) is given for a 10% increase in methylation level. P values were adjusted (Padj) for multiple testing using the Holm correction. Horizontal lines indicate 95% CIs.Fig. S3. Spearman's rank correlation scatter plot of mdNLR and mdNLRref levels in peripheral blood samples of TNBC cases and controls. Fig. S4. Violin plots showing leukocyte subtype proportions in TNBC cases and controls estimated by Houseman method. Fig. S5. Area under the ROC curve of NK cell-to-neutrophil ratio in TNBC cases versus controls. Fig. S6. Kaplan–Meier survival plots of cg08326419 and cg00219921 with prognostic value.

Additional file 9

. Supplementary methods

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Manoochehri, M., Hielscher, T., Borhani, N. et al. Epigenetic quantification of circulating immune cells in peripheral blood of triple-negative breast cancer patients. Clin Epigenet 13, 207 (2021). https://doi.org/10.1186/s13148-021-01196-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-021-01196-1

Keywords

Clinical Epigenetics

ISSN: 1868-7083

Contact us