Next Article in Journal
Controlling Chronic Diseases and Acute Infections with Vitamin D Sufficiency
Next Article in Special Issue
Prevention and Co-Management of Breast Cancer-Related Osteoporosis Using Resveratrol
Previous Article in Journal
Clinical Studies on the Supplementation of Manufactured Human Milk Oligosaccharides: A Systematic Review
Previous Article in Special Issue
Association of Dietary Intake of Zinc and Selenium with Breast Cancer Risk: A Case-Control Study in Chinese Women
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nutrients
Volume 15
Issue 16
10.3390/nu15163621
nutrients-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Plasma One-Carbon Metabolism-Related Micronutrients and the Risk of Breast Cancer: Involvement of DNA Methylation

by
Fubin Liu
†,
Huijun Zhou
†,
Yu Peng
,
Yating Qiao
,
Peng Wang
,
Changyu Si
,
Xixuan Wang
,
Jianxiao Gong
,
Kexin Chen
* and
Fangfang Song
*
Department of Epidemiology and Biostatistics, Key Laboratory of Molecular Cancer Epidemiology, Key Laboratory of Prevention and Control of Major Diseases in the Population, Ministry of Education, National Clinical Research Center for Cancer, Tianjin Medical University Cancer Institute and Hospital, Tianjin Medical University, Tianjin 300060, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Nutrients 2023, 15(16), 3621; https://doi.org/10.3390/nu15163621
Submission received: 30 June 2023 / Revised: 5 August 2023 / Accepted: 16 August 2023 / Published: 17 August 2023
(This article belongs to the Special Issue Micronutrients and Breast Cancer)
Browse Figures
Versions Notes

Abstract

:
Findings of epidemiologic studies focusing on the association between one-carbon metabolism-related micronutrients and breast cancer risk, along with the involvement of DNA methylation, have been inconsistent and incomprehensive. We conducted a case–control study in China including 107 paired participants and comprehensively detected 12 plasma one-carbon metabolism-related micronutrients. Genomic DNA methylation was measured using an 850 K chip and differential methylation probes (DMPs) were identified. Multivariate logistic regression was performed to estimate the associations between plasma micronutrients and the odds of breast cancer. The mediation of selected DMPs in micronutrient breast cancer associations was examined using mediation analyses. An inverse association of plasma folate, methionine cycling-related micronutrients (methionine, S-adenosylmethionine, and S-adenosylhomocysteine), and all micronutrients in the choline metabolism and enzymatic factor groups, and a positive association of methionine cycling-related cysteine with breast cancer risk were observed. Nine micronutrients (methionine, cysteine, SAM, folate, choline, betaine, P5P, vitamins B2, and B12) were related to global or probe-specific methylation levels (p < 0.05). The selected DMPs mediated the micronutrient breast cancer associations with an average mediation proportion of 36.43%. This study depicted comprehensive associations between circulating one-carbon metabolism-related micronutrients and breast cancer risk mediated by DNA methylation.

1. Introduction

Breast cancer is a serious threat to global women’s health. In 2020, approximately 2.26 million new cases and 680 thousand new deaths from breast cancer occurred worldwide [1]. Due to longer life expectancy, urbanization, and the adoption of Western foods and lifestyles, the incidence of breast cancer has increased significantly in developing countries [2]. As a public health problem, breast cancer will reduce the long-term quality of life of patients and bring heavy economic and social burdens [3]. Therefore, identifying modifiable risk factors of breast cancer to improve strategies for risk prevention is essential and urgent.
Dietary factors and nutritional status play an important role in the occurrence and development of breast cancer. Epidemiological studies have shown that the intake of dietary micronutrients involved in the one-carbon metabolism pathway influenced the onset and progression of breast cancer [4,5,6,7,8,9,10]. One-carbon metabolism is a universal metabolic process in cells in which organic groups containing one-carbon atoms are transferred to participate in nucleotide (e.g., purine and thymidylate) biosynthesis, providing methyl groups for all biological methylation reactions. This process is particularly important in the field of cancer metabolism [11]. Two major components of one-carbon metabolism comprise the folate cycle and methionine cycle involving choline metabolism and the participation of enzymatic factors (certain B vitamins) [12]. Some meta-analyses and observational studies have revealed an inverse association between some individual dietary one-carbon metabolism-related micronutrients (e.g., folate, choline, vitamins B2, and B6) and the risk of breast cancer [4,5,9]. Actually, the circulating level in the blood is a better indicator of the body’s overall metabolic load of nutrients than dietary intake. However, inconsistent findings of single circulating one-carbon metabolism-related micronutrients with breast cancer risk remain inconclusive [4,13,14,15,16,17,18,19,20,21,22]. To date, no study has comprehensively depicted an overview of circulating micronutrients in one-carbon metabolism (folate cycle, methionine cycle, choline metabolism, and enzymatic factors) as a functional set in relation to breast cancer risk.
DNA methylation is the most common and important form of epigenetic modification in breast cancer. The epigenome-wide association (EWAS) studies have shown a hypomethylation pattern of peripheral blood genomic DNA existing in breast cancer patients [23,24]. The micronutrients involved in one-carbon metabolism play a crucial part in maintaining genomic stability and methylation patterns [12]. Deficiency or excess of these micronutrients can affect the one-carbon metabolism process, change the availability of S-adenosylmethionine (SAM) in the methionine cycle, interfere with DNA and histone methylation patterns, and further lead to the occurrence and development of breast cancer [12]. Evidence suggested that some B vitamins appeared to be epigenetically active as they are involved in promoter methylation of breast cancer-related genes [25]. However, human data on whether DNA methylation is influenced by one-carbon metabolism-related micronutrients and mediates the association between these micronutrients and breast cancer are limited and inconsistent.
Herein, we aimed to conduct a hospital-based case–control study among Chinese women to examine the association of multiple plasma one-carbon metabolism-related micronutrients with the risk of breast cancer and DNA methylation, and further explore the mediation of DNA methylation in these biochemically-related micronutrients-breast cancer associations.

2. Materials and Methods

2.1. Study Population

The study population consisted of 107 paired breast cancer patients and controls, simply randomly extracted from a hospital-based case–control study we set up [26]. As described in our previous study [27], the cases were patients with newly diagnosed and histologically confirmed breast cancer from Tianjin Medical University Cancer Institute and Hospital since 1 January 2003. They had no history of other malignancies, no reception of any treatment before admission, and no history of blood transfusion in the last 6 months prior to blood sample collection. The age-matched (±1 year) healthy controls were the population without a history of malignancies who received health examination in the community during the same period. The process and purpose of the study were shown in Figure S1.
All participants enrolled in this study were of Chinese Han ethnicity. The Ethics Committee of Tianjin Medical University Cancer Institute and Hospital approved the study protocol (Ek2022212), and written informed consent was obtained from all patients and controls to participate in this study.

2.2. Data and Sample Collection

Information of each subject about demographics (age, education, household income, job, height, and weight), menstrual and reproductive history (age of menarche and first pregnancy, number of pregnancies and live births, breastfeeding months, menopausal status, and history of abortion, birth control pills use, estrogen replacement therapy, and sterilization surgery), environmental exposures (passive smoking and negative events), lifestyles (smoking and alcohol drinking), and family history of cancer in first-degree relatives was obtained from a structured questionnaire. Regarding negative events, participants were questioned, “Have you ever suffered from a negative life event (e.g., miserable marriage, abortion, bereavement, unemployment, illness in the family, change for adaptation, etc.) within the last year that caused a pronounced and lasting negative emotional experience” in the questionnaire, and asked to answer “yes” or “no” according to their own experiences and feelings. Body mass index (BMI) was calculated as weight (kg) divided by the square of height (m). Information on breast density was obtained from imaging examinations.
Five to 10 mL of fasting peripheral blood samples were collected from all subjects. After quiescence for 1–2 h, the blood samples were centrifuged at 3000 rpm for 15 min. Genomic DNA from leukocytes was isolated by the standard phenol–chloroform extraction. The separated plasma and extracted genomic DNA stock solution were stored in the refrigerator at −80 °C for subsequent use.

2.3. Measurements of Plasma One-Carbon Metabolism-Related Micronutrients

The one-carbon metabolism included two major pathways, i.e., folate cycle and methionine cycle, involving choline metabolism and some enzymatic factors (Figure S2). Therefore, we considered 12 one-carbon metabolism-related micronutrients and grouped them into four categories: (1) The methionine cycle group, including methionine, cysteine, homocysteine, SAM, and S-adenosylhomocysteine (SAH); (2) the folate cycle group, including folate and 5-methyltetrahydrofolate (5-MTHF); (3) the choline metabolism group, including choline and betaine; and (4) the enzymatic factor group, including pyridoxal 5-phosphate (P5P, the bioactive form of vitamin B6), vitamin B2, and vitamin B12.
Targeted metabolomic analysis of one-carbon metabolism-related micronutrients in plasma samples was performed using validated liquid chromatography–tandem mass spectrometry (LC-MS/MS) methods [28,29]. Sample injection and separation were performed by a SHIMAZDU 20AXR LC System interfaced with a QTRAP5500+® (AB Sciex, Framingham, MA, USA) [29]. The samples were analyzed following injection of 10 μL of extract on a Phexnomonex Biphenyl 2.6 μm, 3.0 × 50 mm maintained at 35 °C and eluted in a gradient with buffer A (100% water with 0.1% formic acid + 0.01% heptafluorobutyric acid) and buffer B (100% methanol with 0.1% formic acid + 0.01% heptafluorobutyric acid). The flow rate was 0.5 mL/min, with a gradient over a total run time of 7 min: 0.0–2.0 min, 10% B; 2.0–3.0 min, 70% B; 3.0–5.0 min, 100% B; and 5.0–7.0 min, 10% B [28]. Mass spectrum conditions: the electric spray ion source and the negative ion mode were used with the ion source parameter atomization voltage of 5500 (+), the auxiliary heater of 70.0, the ion source temperature at 500.0 °C, the curtain gas of 30.0 and the atomization gas of 50.0. Scanning mode was multiple reaction monitoring. All data were collected and analyzed using Analyst software (version 1.7). Stable isotopes were used for each metabolite to account for matrix effects related to ion–suppression, and the quantitative values obtained were expressed as μmol/L or nmol/L. The accuracy of all metabolites was 85–115%.

2.4. DNA Methylation Assay

The pattern and extent of DNA methylation in the samples were determined using an Illumina Infinium Methylation EPIC BeadChip Kit (850K chip). DNA samples were processed by a series of steps according to the manufacturer’s protocols, i.e., DNA denaturation, whole genome amplification, fragmentation, precipitation, resuspension, hybridization, washing, elongation, staining, and scanning processes, to obtain methylation signals. The extent of methylation was determined by the β value, which was calculated as β = M/(M + U + 100), where M and U represented the signal strength of methylation and unmethylation, respectively. The β values ranged from 0–1, indicating completely unmethylated to totally methylated DNA [30]. The global DNA methylation was calculated as the average β values of all CpGs targeted by the 850K chip.

2.5. Data Processing for DNA Methylation Array

After quality control, we retained 98 cases and 102 controls for subsequent methylation analysis. The raw data of DNA methylation intensity were analyzed using the “ChAMP” package (version 2.21.1) [31]. We read 865,918 probes completely and then removed 9492 probes with a detection p-value above 0.01, 6957 probes with a bead count < 3 in more than 10 samples, 2882 probes with non-CpG sites, 94,837 probes with single nucleotide polymorphisms, 11 probes that aligned to multiple locations, and 16,694 probes located on the X and Y chromosomes. The remaining 735,045 probes were used for further analysis. Furthermore, the beta-mixture quantile (BMIQ) normalization method was used to adjust the intensity distribution of the Infinium II probes to fit that of the Infinium I probes [32], and the ComBat method was used for the correction of multiple batch effects [33]. We used Kolmogorov–Smirnov test methods to compare the distribution of array-wide mean methylation signal between case and control groups [34].
Array-wide differential methylation analysis was carried out using the “limma” package implemented by “champ.DMP” function. A CpG site with |Δβ| > 0.2 (case vs. control) and an adjusted p-value less than 0.05 was defined as a differentially methylated probe (DMP). The heatmap with unsupervised clustering was performed to evaluate the differentiating capacity of DMPs.
We subsequently estimated candidate differentially methylated regions (DMRs) using the “Bumphunter” algorithm [35], where clusters were defined by adjacent CpGs sites spaced up to 500 bp, and 1000 bootstrap samples were used to generate a null distribution of regions. We used the 95% quantile as the cutoff value for candidate regions with at least 7 CpGs probes. Considering its conservative nature relative to Benjamini–Hochberg (BH) and false discovery rate (FDR) methods, family-wise error rate (FWER) was used to control false positive rates with a threshold of <0.001.
To identify the potential biological pathways that DMPs and DMRs are involved in, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed with the “RclusterProfiler” package (Version 4.4.4.) [36]. The significance threshold was set at p < 0.05.

2.6. Statistical Analysis

Continuous variables were compared using the student’s t-test and categorical variables were compared using the Chi-square test. Twelve one-carbon metabolism-related micronutrients were converted to Z scores and then compared by rank-sum tests. A logistic regression model was applied to estimate the odds ratios (ORs) and 95% confidence intervals (CIs) to assess the association between one-carbon metabolism-related micronutrients and breast cancer risk. Combining the results of the distribution of basic characteristics with those well-known risk factors for breast cancer reported in the literature, the model was adjusted for per capita monthly household income, education, body mass index (BMI), negative events, benign breast disease, family history of cancer, age of menarche, menopausal status, number of pregnancies, number of live births, breastfeeding months, and abortion as covariates. We further examined the association between these micronutrients and breast cancer risk stratified by menopausal status.
We investigated correlations of one-carbon metabolism-related micronutrients with array-wide mean methylation levels and selected DMPs (intersection of DMPs and DMRs analyses), with the threshold of <0.05 for p-values corrected by FDR and Bonferroni methods. To ascertain if DNA methylation mediated the association between one-carbon metabolism-related micronutrients and breast cancer, we carried out a mediation analysis using the “mediator” package with adjustment for the same covariates previously identified to evaluate the direct effects, indirect effects, and proportion of mediation by the selected DMPs.
All analyses were performed using SAS version 9.4 (SAS Institute, Cary, NC, USA) and R software (The R Foundation, http://www.r-project.org, accessed on 1 June 2023, version 4.2.3). A level of <0.05 for two-sided p-values was considered statistically significant.

3. Results

3.1. Characteristics of Study Population

The main characteristics of participants by case and control groups are shown in Table 1. The age of the study population was 53.3 ± 6.0 years. Compared to healthy controls, breast cancer cases tended to be poorly educated, obese, with a low-income level, and have a history of benign breast disease, family history of cancer, and negative events. The differences in the distribution of reproductive factors such as pregnancies, live births, breastfeeding months, and abortion between cases and controls were statistically significant.

3.2. Association between One-Carbon Metabolism-Related Micronutrients and Breast Cancer Risk

As shown in Figure 1, except for homocysteine and 5-MTHF, the concentrations of the other 10 one-carbon metabolism-related micronutrients were significantly different between cases and controls (p < 0.05), with the case group presenting lower levels of methionine, SAM and SAH in the methionine cycle group, folate, and all micronutrients in the enzymatic factor group and choline metabolism group, but a higher level of methionine cycling-related cysteine.
The univariate logistic regression analysis showed that 9 one-carbon metabolism-related micronutrients, i.e., methionine, SAM and SAH in the methionine cycle group, folate, and all micronutrients in the enzymatic factor group and choline metabolism group, were associated with a decreased risk of breast cancer, while a higher level of cysteine contributed to an increased breast cancer risk (Table 2). These significant effects remained in the model adjusted for covariates, and the ORs (95% CIs) was 0.59 (0.39–0.92) for methionine, 0.07 (0.02–0.23) for SAM, 0.60 (0.39–0.92) for SAH, 0.36 (0.18–0.70) for folate, 0.02 (0.004–0.17) for P5P, 0.09 (0.03–0.28) for vitamin B2, 0.01 (0.001–0.06) for vitamin B12, 0.43 (0.25–0.73) for choline, 0.36 (0.17–0.75) for betaine, and 1.79 (1.08–2.97) for cysteine, respectively.
Stratified analyses performed according to menopausal status are shown in Table 3. We found no multiplicative interactions between each of the micronutrients and menopausal status in breast cancer. In the fully-adjusted models, the levels of P5P and vitamin B2 had a similar trend of association with breast cancer risk, regardless of menopausal status. However, methionine was only associated with the risk of premenopausal breast cancer, but SAM, folate, vitamin B12, and choline were only linked to the risk of postmenopausal breast cancer.

3.3. Overall and Differential Analysis of DNA Methylation

We observed the global DNA methylation levels were significantly lower in the case group than that in the control group as tested by the Kolmogorov–Smirnov method (p = 2.2 × 10−16; Figure 2A). To further explore the DNA methylation differences between case and control subjects, we filtered the probes with a cut-off of |Δβ|-value > 0.2 and adjusted p-value of <0.05 by a two-group comparison and identified a total of 524 DMPs (Table S1), with a majority (85.3%, 447/524) exhibiting hypomethylation in the case group (Figure S3). Heatmap visualization with hierarchical clustering using these 524 DMPs showed two main clusters: one containing all cases, and the other containing all controls (Figure 2B).
With regard to the analysis of DMRs, we found 920 significant DMRs between case and control groups at a FWER < 0.001 (Table S2).
The GO pathway enrichment analyses revealed that the genes containing DMPs predominantly participated in one pathway, i.e., secretory granule membrane (Table S3). The 45 GO terms were significantly enriched in DMR analysis with the top 10 pathways representing biological processes such as embryonic development, neuron guidance, and cell adhesion (Table S4). No KEGG terms remained significant for DMPs and DMRs.
After an intersection of the sets of DMPs and DMRs based on genome location, 17 DMPs located in the DMRs were selected for further analysis (detailed in Table 4). Among these, most DMPs (88.2%, 15/17) were hypomethylated. Twelve DMPs were located in the promoter region (TSS200, TSS1500, and 1stExon) and five DMPs were located in CpG island.

3.4. Correlation between One-Carbon Metabolism-Related Micronutrients and DNA Methylation

We observed that methionine, SAM, folate, P5P, vitamin B2, vitamin B12, and betaine were positive, but 5-MTHF was negatively correlated with global methylation levels (Figure 3).
Of note is that among the aforementioned 17 selected DMPs, 3 DMPs (cg05397629, cg02589685, and cg10617037) can almost completely distinguish between cases and controls (Figure S4A–C) with area under curve values ≥ 0.999 (Figure S4D–F) and model accuracy ≥ 98.5% (Figure S4G–I). Therefore, they showed a misfit in multivariate logistic regression models for CpG and breast cancer in mediation analyses performed by the “mediator” R package and no results were received. Ultimately, those 3 DMPs were not included in the subsequent analysis. Thus, we explored the correlations between one-carbon metabolism-related micronutrients and 14 DMPs. As shown in Table S5, one-carbon metabolism-related micronutrients were related to some DMPs at a cut-off of 0.05 of unadjusted, FDR, and Bonferroni p-values. Under the Bonferroni standard, all the micronutrients, except for homocysteine, SAH, and 5-MTHF, were significantly related to DMPs, where SAM in the methionine cycle group, folate in the folate cycle group, betaine in the choline metabolism group, P5P, and vitamins B2 and B12 in the enzymatic factor group seemed more epigenetically active exhibiting the remarkable relationships with more than half of the DMPs.

3.5. Mediation of DNA Methylation in the Micronutrients-Breast Cancer Associations

Based on DMPs that were highly correlated with one-carbon metabolism-related micronutrients, the mediation analysis was conducted to examine the mediation of DMPs in the associations between one-carbon metabolism-related micronutrients and breast cancer risk (Figure 4 and Table S6). We observed an average mediation proportion of 36.43%. Specifically, the mediation proportion of DNA methylation was 6.86–72.88% for methionine cycling-related micronutrients including methionine and SAM, 18.95–79.43% for folate, 27.96–69.01% for choline metabolism-related micronutrients, and 19.30–76.79% for enzymatic factors, respectively. Only cg21697381 mediated the choline–breast cancer association. Moreover, the probe cg04794268 was the most widely used mediator in the associations between all micronutrients except choline and breast cancer with a mediation proportion of 46.10–79.43%.

4. Discussion

This case–control study among Chinese women suggested that higher plasma levels of B vitamins (folate, P5P, vitamins B2 and B12), and micronutrients involved in the methionine cycle (methionine, SAM, and SAH) and the choline metabolism (choline and betaine) were associated with a decreased risk of breast cancer, but cysteine level was positively associated with breast cancer risk. Most micronutrients were related to overall or probe-specific methylation levels. Moreover, DNA methylation to a larger extent mediated the micronutrient-breast cancer associations.
The methionine cycle is a major pathway of one-carbon metabolism, accompanied by the folate cycle and choline metabolism which are involved in homocysteine remethylation. The association between plasma homocysteine and the risk of breast cancer has been reported to be null in Western countries [15,16,17,18] and positive in China [14,19], the latter of which was consistent with our results. This signifies that race may be an important factor for inter-study heterogeneity. Population evidence for cysteine and breast cancer is mixed among American women [15,17,18,20]. We supported the positive association by utilizing a Chinese population. Actually, cysteine levels in our study were much lower (median levels were 5.61 μmol/L in controls) relative to those in other American studies (median levels were >200 μmol/L in controls). Discrepant findings among these studies may be attributable to the heterogeneity of the study populations with different exposure loads, whereas no studies based on human evidence were found on circulating levels of methionine, SAM, and SAH and their relation to breast cancer risk. Herein, we first documented inverse associations of these three micronutrients with breast cancer, where the association was more pronounced in premenopausal women for methionine and in postmenopausal women for SAM. A previous study revealed SAM inhibits the migration of triple-negative breast cancer cells in vitro [37], which may provide a theoretical basis for our findings. Epidemiological evidence from Western populations implied that regardless of menopausal status, the risk of breast cancer was contributed to a circulating level of 5-MTHF [13] but not folate [4,13,15,17]. Among Chinese populations, only an early study showed a negative but insignificant association between plasma folate and breast cancer risk [21]. Our results supported the negative association of plasma folate and no association of plasma 5-MTHF with breast cancer risk. The present study population has a higher concentration of folate and 5-MTHF than those in previous studies. Evidence suggested high intake of folate could saturate the activity of dihydrofolate reductase (DHFR), resulting in the presence of 5-MTHF and unmetabolized folate in the circulation [13]. Given this, our population might consume more dietary folate to overwhelm the activity of DHFR. However, the role of the DHFR gene in folic acid metabolism and breast cancer remains to be further investigated. Moreover, methylenetetrahydrofolate reductase (MTHFR) gene C677T (rs1081133) and A1298C (rs1801131) variants were reported to be factors predisposing breast cancer in Asians [38], but not in Western countries [15]. They may directly affect folate metabolite circulation levels to influence breast cancer onset and progression. For circulating choline and betaine, only a higher level of serum betaine contributed to a lower breast cancer risk in the Chinese population [22]. In the present study, we observed both plasma choline and betaine were protective factors for breast cancer, with the former being more pronounced in postmenopausal women. Taken together, more research is needed to clarify the role of one-carbon metabolism involving various functional groups in breast cancer.
Some other B vitamins, such as vitamins B2, B6, and B12, are also key cofactors for one-carbon metabolism. The latest pooled analysis reported that a higher intake of dietary vitamin B2 might reduce the risk of breast cancer, but vitamins B6 and B12 do not [5]. A small number of studies suggested circulating B vitamins including PLP, vitamins B2 and B12 had no effects on breast cancer [5]. However, there is heterogeneity among studies and a lack of evidence on the Chinese population. We first reported inverse relationships between plasma levels of P5P and vitamins B2 and B12 with breast cancer risk in Chinese women. The effects of P5P and vitamin B2 were among both premenopausal and postmenopausal women and that of vitamin B12 was only among postmenopausal women. Nevertheless, more prospective studies and large-scale populations are warranted to support our results.
Hypermethylation of the promoter in cancer suppressor genes, accompanied by global hypomethylation of the genome, is a common feature of tumor cells [39]. Our results showed a genome-wide hypomethylation pattern in the breast cancer group, which was consistent with the previous EWAS studies of breast cancer using a 450 K chip [23,24]. We subsequently used 524 DMPs to form two independent clusters, suggesting the potential of DNA methylation as a classifier for cases and controls and as an optional biomarker in the future. Micronutrients in the one-carbon metabolism pathway have the ability to influence the epigenome and ultimately an individual’s risk of developing cancer [40]. The dietary methyl donors, such as folate, methionine, choline, and betaine, are able to modulate methylation patterns in animal models and humans [25,41]. Generally, we concluded that most one-carbon metabolism-related micronutrients (involving methyl donors and enzymatic factors) were related to array-wide and 14 probe-specific methylation levels. Additionally, we observed the associations between micronutrients from four groups and breast cancer more or less mediated by DNA methylation. To date, there is little evidence for a mediating role of DNA methylation in the association between these micronutrients and breast cancer. Nevertheless, our findings further clarified the possible involvement of DNA methylation in the process by which micronutrients involved in various parts of one-carbon metabolism influence the onset of breast cancer.
The hypomethylated probe cg04794268, located in the promoter of KIFC2, widely mediated the associations between almost all micronutrients and breast cancer in our study, although no report about the relationship between KIFC2 and breast cancer has been found so far. Another interesting finding was the association of choline with a decreased risk of breast cancer only mediated by the hypomethylation of cg21697381, located in the open sea region of the SLFN12 promoter. The over-expression of SLFN12 can sensitize triple-negative breast cancer cells to chemotherapy drugs and radiotherapy [42]. However, whether the SLFN12 gene is related to a reduced risk of breast cancer conferred by choline needs to be further investigated. Together, our findings make it interesting to systematically evaluate the role of DNA methylation modulation by one-carbon metabolism-related micronutrients in breast cancer in future studies.
To our knowledge, this was the first study to systematically consider one-carbon metabolism-related micronutrients as functional groups (methionine cycle, folate cycle, choline metabolism, and enzymatic factor) and comprehensively evaluate the potential role of plasma levels of these micronutrients in breast cancer and the involvement of DNA methylation in a Chinese population. Moreover, the whole genome methylation levels were first detected by an 850 K chip in a breast cancer study, providing wider coverage, higher data quality, and higher reproducibility. However, several limitations existing in the present study merited special consideration. Firstly, our study had the inherent limitations of a case–control study, such as recall bias and selection bias. Secondly, the main population in this study was citizens of Tianjin (a city in northern China), which may not be representative of the entire Chinese population. Thirdly, even though we had controlled the majority of confounders, the residual confounding from unmeasured or unknown factors might remain. Finally, due to the relatively small sample size of this study, the current analysis may be underpowered and the conclusions need to be verified by large-scale prospective data.

5. Conclusions

Our findings provided the first comprehensive depiction of an inverse association of plasma folate, micronutrients involved in the methionine cycle (methionine, SAM, and SAH) and choline metabolism (choline and betaine) and enzymatic factors in one-carbon metabolism (P5P, vitamins B2 and B12) and a positive association of methionine cycling-related cysteine with breast cancer risk in a Chinese population. Certain one-carbon metabolism-related micronutrients were related to array-wide and probe-specific methylation levels. Furthermore, differentially methylated CpGs partially or completely mediated the micronutrient-breast cancer associations. The above findings enrich the theoretical basis for the relationship between one-carbon metabolism-related micronutrients and individual methylation sites and provide evidence that one-carbon metabolism-related micronutrients may specifically affect breast cancer risk through different methylation sites. Meanwhile, one-carbon metabolism-related micronutrients may serve as a potential epigenetic drug to reverse abnormal methylation patterns in breast cancer, contributing to the development of targeted prevention and control measures for breast cancer. More prospective studies are needed in the future to clarify those associations, as well as to be motivated to uncover the specific mechanisms of action.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nu15163621/s1, Figure S1: flowchart of the study; Figure S2: overview of one-carbon metabolism and its related micronutrients; Figure S3: volcano plot of analysis of differential methylation probes (DMPs); Figure S4: details about distinguishing capacity of three differential methylation probes (DMPs); Table S1: information of 524 differentially methylation probes (DMPs); Table S2: information of 920 differentially methylation regions (DMRs); Table S3: gene ontology terms that were significantly enriched in an analysis of differentially methylation probes (DMPs); Table S4. gene ontology terms that were significantly enriched in an analysis of differentially methylation regions (DMRs); Table S5: correlation analysis between plasma one-carbon metabolism-related micronutrients and 14 differential methylation probes (DMPs); Table S6: mediation analysis of differential methylation probes (DMPs) in the association between plasma one-carbon metabolism-related micronutrients and risk of breast cancer.

Author Contributions

Conceptualization: F.S. and K.C.; Data Curation: F.L.; Formal Analysis: F.L., Y.P. and Y.Q.; Funding Acquisition: F.S.; Project Administration: F.S.; Software: H.Z. and P.W.; Supervision: F.S. and K.C.; Validation: C.S. and X.W.; Visualization: F.L., H.Z. and J.G.; Writing—Original Draft: F.L. and H.Z.; Writing—Review and Editing: F.L., H.Z., Y.P., Y.Q., P.W., C.S., X.W., J.G., F.S. and K.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Program of China (2021YFC2500401), the National Natural Science Foundation of China (81974488), and the Tianjin Key Medical Discipline (Specialty) Construction Project (TJYXZDXK-009A).

Institutional Review Board Statement

The study was approved by the Ethics Committee of Tianjin Medical University Cancer Institute and Hospital (Ek2022212) and conducted in accordance with the Declaration of Helsinki.

Informed Consent Statement

Written informed consent was obtained from all patients and controls to participate in this study.

Data Availability Statement

The data used for this analysis can be made available upon reasonable request to the corresponding authors.

Acknowledgments

We would like to thank the Cancer Biobank of Tianjin Medical University Cancer Institute and Hospital, founded by the National Human Genetic Resources Sharing Service Platform (2005DKA21300).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

5-MTHF, 5-methyltetrahydrofolate; BH, Benjamini–Hochberg; BMI, body mass index; BMIQ, beta-mixture quantile; CI, confidence interval; DHFR, dihydrofolate reductase; DMP, differentially methylated probe; DMR, differentially methylated region; EWAS, epigenome-wide association; FDR, false discovery rate; FWER, family-wise error rate; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; LC-MS/MS, liquid chromatography–tandem mass spectrometry; MTHFR, methylenetetrahydrofolate reductase; OR, odds ratio; P5P, pyridoxal 5-phosphate; SAH, S-adenosylhomocysteine; SAM, S-adenosylmethionine.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  2. De Cicco, P.; Catani, M.V.; Gasperi, V.; Sibilano, M.; Quaglietta, M.; Savini, I. Nutrition and Breast Cancer: A Literature Review on Prevention, Treatment and Recurrence. Nutrients 2019, 11, 1514. [Google Scholar] [CrossRef] [PubMed]
  3. Chen, W.; Zheng, R.; Baade, P.D.; Zhang, S.; Zeng, H.; Bray, F.; Jemal, A.; Yu, X.Q.; He, J. Cancer statistics in China, 2015. CA Cancer J. Clin. 2016, 66, 115–132. [Google Scholar] [CrossRef] [PubMed]
  4. Ren, X.; Xu, P.; Zhang, D.; Liu, K.; Song, D.; Zheng, Y.; Yang, S.; Li, N.; Hao, Q.; Wu, Y.; et al. Association of folate intake and plasma folate level with the risk of breast cancer: A dose-response meta-analysis of observational studies. Aging 2020, 12, 21355–21375. [Google Scholar] [CrossRef]
  5. Zeng, J.; Gu, Y.; Fu, H.; Liu, C.; Zou, Y.; Chang, H. Association Between One-carbon Metabolism-related Vitamins and Risk of Breast Cancer: A Systematic Review and Meta-analysis of Prospective Studies. Clin. Breast Cancer 2020, 20, e469–e480. [Google Scholar] [CrossRef]
  6. Song, H.; Jeong, A.; Tran, T.X.M.; Lee, J.; Kim, M.; Park, B. Association between Micronutrient Intake and Breast Cancer Risk According to Body Mass Index in South Korean Adult Women: A Cohort Study. Nutrients 2022, 14, 2644. [Google Scholar] [CrossRef]
  7. Sun, Y.; Fowke, J.H.; Liang, X.; Mozhui, K.; Sen, S.; Bao, W.; Liu, B.; Snetselaar, L.G.; Wallace, R.B.; Shadyab, A.H.; et al. Changes in Dietary Intake of Methionine, Folate/Folic Acid and Vitamin B12 and Survival in Postmenopausal Women with Breast Cancer: A Prospective Cohort Study. Nutrients 2022, 14, 4747. [Google Scholar] [CrossRef]
  8. Hatami, M.; Vahid, F.; Akbari, M.E.; Sadeghi, M.; Ameri, F.; Eini-Zeinab, H.; Jamshidi-Naeini, Y.; Davoodi, S.H. The Vitamins Involved in One-Carbon Metabolisms are Associated with Reduced Risk of Breast Cancer in Overall and Subtypes. Int. J. Vitam. Nutr. Res. 2020, 90, 131–140. [Google Scholar] [CrossRef]
  9. Van Puyvelde, H.; Dimou, N.; Katsikari, A.; Indave Ruiz, B.I.; Godderis, L.; Huybrechts, I.; De Bacquer, D. The association between dietary intakes of methionine, choline and betaine and breast cancer risk: A systematic review and meta-analysis. Cancer Epidemiol. 2023, 83, 102322. [Google Scholar] [CrossRef]
  10. Wu, W.; Kang, S.; Zhang, D. Association of vitamin B6, vitamin B12 and methionine with risk of breast cancer: A dose-response meta-analysis. Br. J. Cancer 2013, 109, 1926–1944. [Google Scholar] [CrossRef]
  11. Ducker, G.S.; Rabinowitz, J.D. One-Carbon Metabolism in Health and Disease. Cell Metab. 2017, 25, 27–42. [Google Scholar] [CrossRef] [PubMed]
  12. Mentch, S.J.; Locasale, J.W. One-carbon metabolism and epigenetics: Understanding the specificity. Ann. N. Y. Acad. Sci. 2016, 1363, 91–98. [Google Scholar] [CrossRef] [PubMed]
  13. Koenig, K.L.; Scarmo, S.; Afanasyeva, Y.; Clendenen, T.V.; Ueland, P.M.; Zeleniuch-Jacquotte, A. Circulating unmetabolized folic acid and 5-methyltetrahydrofolate and risk of breast cancer: A nested case-control study. Eur. J. Clin. Nutr. 2020, 74, 1306–1315. [Google Scholar] [CrossRef]
  14. Wu, X.; Zou, T.; Cao, N.; Ni, J.; Xu, W.; Zhou, T.; Wang, X. Plasma homocysteine levels and genetic polymorphisms in folate metablism are associated with breast cancer risk in chinese women. Hered. Cancer Clin. Pract. 2014, 12, 2. [Google Scholar] [CrossRef] [PubMed]
  15. Houghton, S.C.; Eliassen, A.H.; Zhang, S.M.; Selhub, J.; Rosner, B.A.; Willett, W.C.; Hankinson, S.E. Plasma B-vitamins and one-carbon metabolites and the risk of breast cancer in younger women. Breast Cancer Res. Treat. 2019, 176, 191–203. [Google Scholar] [CrossRef] [PubMed]
  16. He, Q.; Yang, Z.; Sun, Y.; Qu, Z.; Jia, X.; Li, J.; Lin, Y.; Luo, Y. The Impact of Homocysteine on the Risk of Hormone-Related Cancers: A Mendelian Randomization Study. Front. Nutr. 2021, 8, 645371. [Google Scholar] [CrossRef]
  17. Houghton, S.C.; Eliassen, A.H.; Zhang, S.M.; Selhub, J.; Rosner, B.A.; Willett, W.C.; Hankinson, S.E. Plasma B-vitamin and one-carbon metabolites and risk of breast cancer before and after folic acid fortification in the United States. Int. J. Cancer 2019, 144, 1929–1940. [Google Scholar] [CrossRef]
  18. Lin, J.; Lee, I.M.; Song, Y.; Cook, N.R.; Selhub, J.; Manson, J.E.; Buring, J.E.; Zhang, S.M. Plasma homocysteine and cysteine and risk of breast cancer in women. Cancer Res. 2010, 70, 2397–2405. [Google Scholar] [CrossRef]
  19. Chou, Y.C.; Lee, M.S.; Wu, M.H.; Shih, H.L.; Yang, T.; Yu, C.P.; Yu, J.C.; Sun, C.A. Plasma homocysteine as a metabolic risk factor for breast cancer: Findings from a case-control study in Taiwan. Breast Cancer Res. Treat. 2007, 101, 199–205. [Google Scholar] [CrossRef]
  20. Zhang, S.M.; Willett, W.C.; Selhub, J.; Manson, J.E.; Colditz, G.A.; Hankinson, S.E. A prospective study of plasma total cysteine and risk of breast cancer. Cancer Epidemiol. Biomark. Prev. 2003, 12, 1188–1193. [Google Scholar]
  21. Chou, Y.C.; Wu, M.H.; Yu, J.C.; Lee, M.S.; Yang, T.; Shih, H.L.; Wu, T.Y.; Sun, C.A. Genetic polymorphisms of the methylenetetrahydrofolate reductase gene, plasma folate levels and breast cancer susceptibility: A case-control study in Taiwan. Carcinogenesis 2006, 27, 2295–2300. [Google Scholar] [CrossRef] [PubMed]
  22. Du, Y.F.; Lin, F.Y.; Long, W.Q.; Luo, W.P.; Yan, B.; Xu, M.; Mo, X.F.; Zhang, C.X. Serum betaine but not choline is inversely associated with breast cancer risk: A case-control study in China. Eur. J. Nutr. 2017, 56, 1329–1337. [Google Scholar] [CrossRef] [PubMed]
  23. Severi, G.; Southey, M.C.; English, D.R.; Jung, C.H.; Lonie, A.; McLean, C.; Tsimiklis, H.; Hopper, J.L.; Giles, G.G.; Baglietto, L. Epigenome-wide methylation in DNA from peripheral blood as a marker of risk for breast cancer. Breast Cancer Res. Treat. 2014, 148, 665–673. [Google Scholar] [CrossRef]
  24. van Veldhoven, K.; Polidoro, S.; Baglietto, L.; Severi, G.; Sacerdote, C.; Panico, S.; Mattiello, A.; Palli, D.; Masala, G.; Krogh, V.; et al. Epigenome-wide association study reveals decreased average methylation levels years before breast cancer diagnosis. Clin. Epigenet. 2015, 7, 67. [Google Scholar] [CrossRef]
  25. Xu, X.; Gammon, M.D.; Jefferson, E.; Zhang, Y.; Cho, Y.H.; Wetmur, J.G.; Teitelbaum, S.L.; Bradshaw, P.T.; Terry, M.B.; Garbowski, G.; et al. The influence of one-carbon metabolism on gene promoter methylation in a population-based breast cancer study. Epigenetics 2011, 6, 1276–1283. [Google Scholar] [CrossRef] [PubMed]
  26. Song, F.; Ji, P.; Zheng, H.; Song, F.; Wang, Y.; Hao, X.; Wei, Q.; Zhang, W.; Chen, K. Definition of a functional single nucleotide polymorphism in the cell migration inhibitory gene MIIP that affects the risk of breast cancer. Cancer Res. 2010, 70, 1024–1032. [Google Scholar] [CrossRef] [PubMed]
  27. Peng, Y.; Liu, F.; Qiao, Y.; Wang, P.; Du, H.; Si, C.; Wang, X.; Chen, K.; Song, F. Genetically Modified Circulating Levels of Advanced Glycation End-Products and Their Soluble Receptor (AGEs-RAGE Axis) with Risk and Mortality of Breast Cancer. Cancers 2022, 14, 6124. [Google Scholar] [CrossRef] [PubMed]
  28. Butler, L.M.; Arning, E.; Wang, R.; Bottiglieri, T.; Govindarajan, S.; Gao, Y.T.; Yuan, J.M. Prediagnostic levels of serum one-carbon metabolites and risk of hepatocellular carcinoma. Cancer Epidemiol. Biomark. Prev. 2013, 22, 1884–1893. [Google Scholar] [CrossRef]
  29. Spratlen, M.J.; Grau-Perez, M.; Umans, J.G.; Yracheta, J.; Best, L.G.; Francesconi, K.; Goessler, W.; Bottiglieri, T.; Gamble, M.V.; Cole, S.A.; et al. Targeted metabolomics to understand the association between arsenic metabolism and diabetes-related outcomes: Preliminary evidence from the Strong Heart Family Study. Environ. Res. 2019, 168, 146–157. [Google Scholar] [CrossRef]
  30. Zhang, Y.; Long, H.; Wang, S.; Xiao, W.; Xiong, M.; Liu, J.; Chen, L.; Chen, R.; Wei, X.; Shu, Y.; et al. Genome-Wide DNA Methylation Pattern in Whole Blood Associated With Primary Intracerebral Hemorrhage. Front. Immunol. 2021, 12, 702244. [Google Scholar] [CrossRef]
  31. Morris, T.J.; Butcher, L.M.; Feber, A.; Teschendorff, A.E.; Chakravarthy, A.R.; Wojdacz, T.K.; Beck, S. ChAMP: 450k Chip Analysis Methylation Pipeline. Bioinformatics 2014, 30, 428–430. [Google Scholar] [CrossRef] [PubMed]
  32. Teschendorff, A.E.; Marabita, F.; Lechner, M.; Bartlett, T.; Tegner, J.; Gomez-Cabrero, D.; Beck, S. A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data. Bioinformatics 2013, 29, 189–196. [Google Scholar] [CrossRef]
  33. Johnson, W.E.; Li, C.; Rabinovic, A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 2007, 8, 118–127. [Google Scholar] [CrossRef] [PubMed]
  34. Li, D.; Xie, Z.; Pape, M.L.; Dye, T. An evaluation of statistical methods for DNA methylation microarray data analysis. BMC Bioinform. 2015, 16, 217. [Google Scholar] [CrossRef] [PubMed]
  35. Jaffe, A.E.; Murakami, P.; Lee, H.; Leek, J.T.; Fallin, M.D.; Feinberg, A.P.; Irizarry, R.A. Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies. Int. J. Epidemiol. 2012, 41, 200–209. [Google Scholar] [CrossRef]
  36. Yu, G.; Wang, L.G.; Han, Y.; He, Q.Y. clusterProfiler: An R package for comparing biological themes among gene clusters. Omics 2012, 16, 284–287. [Google Scholar] [CrossRef] [PubMed]
  37. Coppola, A.; Ilisso, C.P.; Stellavato, A.; Schiraldi, C.; Caraglia, M.; Mosca, L.; Cacciapuoti, G.; Porcelli, M. S-Adenosylmethionine Inhibits Cell Growth and Migration of Triple Negative Breast Cancer Cells through Upregulating MiRNA-34c and MiRNA-449a. Int. J. Mol. Sci. 2020, 22, 286. [Google Scholar] [CrossRef]
  38. Rezaee, M.; Akbari, H.; Momeni-Moghaddam, M.A.; Moazzen, F.; Salahi, S.; Jahankhah, R.; Tahmasebi, S. Association of C677T (rs1081133) and A1298C (rs1801131) Methylenetetrahydrofolate Reductase Variants with Breast Cancer Susceptibility Among Asians: A Systematic Review and Meta-Analysis. Biochem. Genet. 2021, 59, 367–397. [Google Scholar] [CrossRef]
  39. Esteller, M. Epigenetics in cancer. N. Engl. J. Med. 2008, 358, 1148–1159. [Google Scholar] [CrossRef]
  40. Chen, J.; Xu, X. Diet, epigenetic, and cancer prevention. Adv. Genet. 2010, 71, 237–255. [Google Scholar] [CrossRef]
  41. Mahmoud, A.M.; Ali, M.M. Methyl Donor Micronutrients that Modify DNA Methylation and Cancer Outcome. Nutrients 2019, 11, 608. [Google Scholar] [CrossRef] [PubMed]
  42. Raafat Elsayed, A.A.; Al-Marsoummi, S.; Vomhof-Dekrey, E.E.; Basson, M.D. SLFN12 Over-expression Sensitizes Triple Negative Breast Cancer Cells to Chemotherapy Drugs and Radiotherapy. Cancer Genom. Proteom. 2022, 19, 328–338. [Google Scholar] [CrossRef] [PubMed]
Nutrients 15 03621 g001
Figure 1. Comparison of plasma one-carbon metabolism-related micronutrients between case and control population. (A) Methionine cycle group, including methionine, cysteine, homocysteine, SAM, and SAH. (B) Folate cycle group, including folate and 5-MTHF. (C) Enzymatic factor group, including P5P and vitamins B2 and B12. (D) Choline metabolism group, including choline and betaine. The p-value represents the significance of rank-sum test. * <0.05; ** <0.01; *** <0.001; **** <0.0001. Abbreviations: 5-MTHF, 5-methyltetrahydrofolate; ns, not significant; P5P, pyridoxal 5-phosphate; SAH, S-adenosylhomocysteine; SAM, S-adenosylmethionine.
Figure 1. Comparison of plasma one-carbon metabolism-related micronutrients between case and control population. (A) Methionine cycle group, including methionine, cysteine, homocysteine, SAM, and SAH. (B) Folate cycle group, including folate and 5-MTHF. (C) Enzymatic factor group, including P5P and vitamins B2 and B12. (D) Choline metabolism group, including choline and betaine. The p-value represents the significance of rank-sum test. * <0.05; ** <0.01; *** <0.001; **** <0.0001. Abbreviations: 5-MTHF, 5-methyltetrahydrofolate; ns, not significant; P5P, pyridoxal 5-phosphate; SAH, S-adenosylhomocysteine; SAM, S-adenosylmethionine.
Nutrients 15 03621 g001
Nutrients 15 03621 g002
Figure 2. Overall and differential analysis of DNA methylation between breast cancer cases and control subjects. (A) Distribution of array-wide mean methylation level between breast cancer cases and control population. The p-value represents the significance of Kolmogorov–Smirnov test. (B) Heatmap of 524 differential methylation probes with unsupervised clustering.
Figure 2. Overall and differential analysis of DNA methylation between breast cancer cases and control subjects. (A) Distribution of array-wide mean methylation level between breast cancer cases and control population. The p-value represents the significance of Kolmogorov–Smirnov test. (B) Heatmap of 524 differential methylation probes with unsupervised clustering.
Nutrients 15 03621 g002
Nutrients 15 03621 g003
Figure 3. Correlation analysis between plasma one-carbon metabolism-related micronutrients and array-wide methylation levels. (A) Methionine cycle group, including methionine, cysteine, homocysteine, SAM, and SAH. (B) Folate cycle group, including folate and 5-MTHF. (C) Enzymatic factor group, including P5P, vitamins B2, and B12. (D) Choline metabolism group, including choline and betaine. The Spearman correlation coefficient and p-value are shown. Abbreviations: 5-MTHF, 5-methyltetrahydrofolate; P5P, pyridoxal 5-phosphate; SAH, S-adenosylhomocysteine; and SAM, S-adenosylmethionine.
Figure 3. Correlation analysis between plasma one-carbon metabolism-related micronutrients and array-wide methylation levels. (A) Methionine cycle group, including methionine, cysteine, homocysteine, SAM, and SAH. (B) Folate cycle group, including folate and 5-MTHF. (C) Enzymatic factor group, including P5P, vitamins B2, and B12. (D) Choline metabolism group, including choline and betaine. The Spearman correlation coefficient and p-value are shown. Abbreviations: 5-MTHF, 5-methyltetrahydrofolate; P5P, pyridoxal 5-phosphate; SAH, S-adenosylhomocysteine; and SAM, S-adenosylmethionine.
Nutrients 15 03621 g003
Nutrients 15 03621 g004
Figure 4. Mediation analysis of differential methylation probes (DMPs) in the association between plasma one-carbon metabolism-related micronutrients and risk of breast cancer. The mediation proportion (%) and significance of 14 DMPs for each one-carbon metabolism-related micronutrient are shown. The blank areas indicate that the micronutrients were not associated with the selected DMPs and thus not included in the mediation analysis. The areas marked with an “X” indicate an insignificant mediation proportion in the mediation analysis. * p < 0.05; *** p < 0.001. Abbreviations: ns, not significant; P5P, pyridoxal 5-phosphate; and SAM, S-adenosylmethionine.
Figure 4. Mediation analysis of differential methylation probes (DMPs) in the association between plasma one-carbon metabolism-related micronutrients and risk of breast cancer. The mediation proportion (%) and significance of 14 DMPs for each one-carbon metabolism-related micronutrient are shown. The blank areas indicate that the micronutrients were not associated with the selected DMPs and thus not included in the mediation analysis. The areas marked with an “X” indicate an insignificant mediation proportion in the mediation analysis. * p < 0.05; *** p < 0.001. Abbreviations: ns, not significant; P5P, pyridoxal 5-phosphate; and SAM, S-adenosylmethionine.
Nutrients 15 03621 g004
Table 1. Baseline characteristics of 107 paired case and control population.
Table 1. Baseline characteristics of 107 paired case and control population.
VariableGroupCase (N = 107)Control (N = 107)p
Age, mean ± SD, years 53.3 ± 6.353.3 ± 5.80.991
Education, N (%)Senior high school or above55 (55.0)66 (61.7)0.043
Junior high school30 (30.0)36 (33.6)
Primary school and below15 (15.0)5 (4.7)
Per capita monthly household income, N (%), CNY<5004 (4.1)2 (1.9)<0.001
500–99915 (15.5)6 (5.7)
1000–199947 (48.5)32 (30.2)
2000–299920 (20.6)33 (31.1)
≥300011 (11.3)33 (31.1)
Job, N (%)No53 (54.6)55 (52.9)0.803
Yes44 (45.4)49 (47.1)
BMI, N (%), kg/m2<23.935 (32.7)55 (51.9)0.001
24.0–27.946 (43.0)43 (40.6)
≥28.026 (24.3)8 (7.5)
Benign breast disease, N (%)No66 (65.3)83 (82.2)0.007
Yes35 (34.7)18 (17.8)
Family history of cancer, N (%)No66 (62.3)79 (75.2)0.042
Yes40 (37.7)26 (24.8)
Family history of breast cancer, N (%)No97 (94.2)104 (99.0)0.087
Yes6 (5.8)1 (1.0)
Smoking, N (%)No101 (97.1)97 (96.0)0.673
Yes3 (2.9)4 (4.0)
Passive smoking, N (%)No33 (41.3)31 (55.4)0.106
Yes47 (58.8)25 (44.6)
Alcohol drinking, N (%)No102 (98.1)100 (98.0)0.984
Yes2 (1.9)2 (2.0)
Negative events, N (%)No74 (75.5)89 (89.9)0.008
Yes24 (24.5)10 (10.1)
Age of menarche, N (%), years≤1327 (26.0)18 (17.1)0.115
14–1541 (39.4)39 (37.1)
≥1636 (34.6)48 (45.7)
Number of pregnancies, N (%)≥353 (51.0)27 (25.7)<0.001
238 (36.5)35 (33.3)
≤113 (12.5)43 (41.0)
Age of first pregnancy, N (%), years<3087 (87.9)86 (85.1)0.572
≥3012 (12.1)15 (14.9)
Number of live births, N (%)≥230 (29.1)13 (12.4)0.003
≤173 (70.9)92 (87.6)
Breastfeeding months, N (%), months>1263 (63.0)43 (43.0)0.005
≤1237 (37.0)57 (57.0)
Abortion, N (%)No18 (17.3)42 (41.2)<0.001
Yes86 (82.7)60 (58.8)
Take birth control pills, N (%)No84 (84.8)84 (82.4)0.633
Yes15 (15.2)18 (17.6)
Estrogen replacement therapy, N (%)No87 (89.7)85 (94.4)0.232
Yes10 (10.3)5 (5.6)
Sterilization surgery, N (%)No93 (91.2)83 (91.2)0.994
Yes9 (8.8)8 (8.8)
Menopausal status, N (%)Premenopausal34 (32.1)33 (32.0)0.995
Postmenopausal72 (67.9)70 (68.0)
Breast density, N (%)Low72 (75.8)51 (63.0)0.064
High23 (24.2)30 (37.0)
Abbreviations: BMI, body mass index; SD, standard deviation.
Table 2. Association of plasma one-carbon metabolism-related micronutrients with risk of breast cancer.
Table 2. Association of plasma one-carbon metabolism-related micronutrients with risk of breast cancer.
GroupMicronutrientsMedian (P25, P75)OR (95% CI)
Case (N = 107)Control (N = 107)CrudeAdjusted
Methionine cycleMethionine, μmol/L18.06 (9.27,37.02)31.33 (15.50,55.15)0.62 (0.46–0.83)0.59 (0.39–0.92)
Cysteine, μmol/L6.95 (4.63,9.57)5.61 (4.53,8.45)1.44 (1.08–1.91)1.79 (1.08–2.97)
Homocysteine, μmol/L4.91 (3.13,9.77)6.47 (3.38,16.25)0.75 (0.49–1.13)0.84 (0.51–1.40)
SAM, μmol/L0.07 (0.04,0.13)0.20 (0.11,0.37)0.12 (0.06–0.26)0.07 (0.02–0.23)
SAH, μmol/L0.04 (0.02,0.06)0.05 (0.04,0.07)0.70 (0.53–0.94)0.60 (0.39–0.92)
Folate cycleFolate, nmol/L8.75 (4.79,15.35)20.63 (8.37,36.57)0.32 (0.21–0.50)0.34 (0.18–0.70)
5-MTHF, μmol/L0.24 (0.13,0.57)0.26 (0.16,0.47)1.09 (0.83–1.42)1.09 (0.72–1.65)
Enzymatic factorP5P, nmol/L9.94 (4.39,16.07)28.55 (15.83,42.29)0.10 (0.05–0.19)0.02 (0.004–0.17)
Vitamin B2, nmol/L2.45 (1.39,3.85)6.87 (3.67,11.34)0.13 (0.07–0.24)0.09 (0.03–0.28)
Vitamin B12, nmol/L1.23 (0.69,1.97)3.76 (2.16,6.84)0.02 (0.01–0.07)0.01 (0.001–0.06)
Choline metabolismCholine, μmol/L6.58 (4.62,10.64)11.84 (6.15,20.17)0.33 (0.20–0.53)0.43 (0.25–0.73)
Betaine, μmol/L15.26 (9.63,22.47)23.83 (13.55,37.73)0.42 (0.28–0.62)0.36 (0.17–0.75)
Adjusted OR was adjusted for per capita monthly household income (CNY < 500, CNY 500–999, CNY 1000–1999, CNY 2000–2999, and CNY ≥ 3000); education (senior high school or above, junior high school, primary school and below); body mass index (<23.9 kg/m2, 24.0–27.9 kg/m2, ≥28.0 kg/m2); benign breast disease (no, yes); family history of cancer (no yes); negative events (no, yes); age of menarche (≤13 years, 14–15 years, and ≥16 years); menopausal status (premenopausal, postmenopausal); number of pregnancies (≥3, 2, ≤1); number of live births (≥2, ≤1); breastfeeding months (>12 months, ≤12 months); and abortion (no, yes). Abbreviations: 5-MTHF, 5-methyltetrahydrofolate; CI, confidence interval; OR, odds ratio; P5P, pyridoxal 5-phosphate; SAM, S-adenosylmethionine; and SAH, S-adenosylhomocysteine.
Table 3. Association of plasma one-carbon metabolism-related micronutrients with risk of breast cancer stratified by menopausal status.
Table 3. Association of plasma one-carbon metabolism-related micronutrients with risk of breast cancer stratified by menopausal status.
GroupMicronutrientsMedian (P25, P75)OR (95% CI)
CaseControlCrudeAdjusted
Premenopausal (N = 67)
Methionine cycleMethionine, μmol/L18.62 (9.43,51.90)41.85 (19.70,59.58)0.62 (0.37–1.03)0.07 (0.004–0.99)
Cysteine, μmol/L7.44 (4.62–9.80)4.96 (4.48–6.91)2.30 (1.20–4.40)0.98 (0.39–2.46)
Homocysteine, μmol/L5.59 (3.09–18.97)11.62 (4.66–18.97)0.71 (0.42–1.18)0.46 (0.08–2.55)
SAM, μmol/L0.09 (0.03–0.16)0.17 (0.12–0.26)0.16 (0.04–0.70)0.19 (0.02–1.52)
SAH, μmol/L0.04 (0.02–0.05)0.05 (0.04–0.07)0.68 (0.41–1.12)0.37 (0.08–1.72)
Folate cycleFolate, nmol/L7.86 (4.21–13.51)11.80 (6.96–34.68)0.36 (0.17–0.78)0.21 (0.01–4.50)
5-MTHF, μmol/L0.20 (0.13–0.62)0.24 (0.15–0.34)1.42 (0.84–2.39)0.32 (0.07–1.41)
Enzymatic factorP5P, nmol/L11.61 (4.80–18.32)32.35 (16.66–50.75)0.17 (0.06–0.43)0.01 (0.00–0.82)
Vitamin B2, nmol/L2.81 (1.19–3.89)6.18 (3.31–8.49)0.18 (0.07–0.47)0.01 (0.00–0.83)
Vitamin B12, nmol/L1.31 (0.73–2.01)3.89 (2.45–8.00)0.01 (0.00–0.10)0.00 (0.00–6.33)
Choline metabolismCholine, μmol/L6.09 (4.03–8.76)11.02 (5.98–23.63)0.18 (0.05–0.63)0.29 (0.05–1.71)
Betaine, μmol/L14.64 (9.74–23.77)23.64 (16.75–37.16)0.39 (0.19–0.81)0.01 (0.00–1.44)
Postmenopausal (N = 143)
Methionine cycleMethionine, μmol/L16.27 (9.14–34.30)30.56 (15.37–51.28)0.60 (0.42–0.86)0.78 (0.38–1.59)
Cysteine, μmol/L6.95 (4.63–9.56)6.06 (4.59–8.68)1.21 (0.87–1.69)1.97 (0.84–4.66)
Homocysteine, μmol/L4.91 (3.13–8.84)4.95 (2.67–13.62)0.86 (0.61–1.22)1.03 (0.53–1.98)
SAM, μmol/L0.06 (0.04–0.11)0.23 (0.11–0.47)0.06 (0.02–0.19)0.02 (0.002–0.21)
SAH, μmol/L0.04 (0.02–0.06)0.05 (0.04–0.07)0.73 (0.52–1.04)0.72 (0.34–1.51)
Folate cycleFolate, nmol/L10.44 (4.98–16.50)20.88 (10.98–40.04)0.30 (0.17–0.53)0.34 (0.13–0.92)
5-MTHF, μmol/L0.28 (0.13–0.46)0.26 (0.16–0.48)1.04 (0.74–1.45)1.42 (0.69–2.91)
Enzymatic factorP5P, nmol/L8.98 (4.26–15.07)25.49 (14.23–40.81)0.08 (0.03–0.20)0.02 (0.001–0.20)
Vitamin B2, nmol/L2.48 (1.42–3.94)7.27 (4.15–12.65)0.12 (0.05–0.26)0.09 (0.02–0.53)
Vitamin B12, nmol/L1.18 (0.61–1.97)3.29 (1.94–5.64)0.09 (0.04–0.22)0.01 (0.00–0.21)
Choline metabolismCholine, μmol/L7.39 (4.62–11.32)11.73 (6.15–18.78)0.54 (0.35–0.83)0.36 (0.14–0.90)
Betaine, μmol/L15.60 (9.54–21.02)23.85 (12.05–36.48)0.53 (0.36–0.79)0.58 (0.20–1.64)
All p-values for multiplicative interactions between micronutrients and menopausal status were >0.05. Adjusted OR was adjusted for per capita monthly household income (CNY < 500, CNY 500–999, CNY 1000–1999, CNY 2000–2999, and CNY ≥ 3000); education (senior high school or above, junior high school, primary school and below); BMI (<23.9 kg/m2, 24.0–27.9 kg/m2, ≥28.0 kg/m2); benign breast disease (no, yes); family history of cancer (no yes); negative events (no, yes); age of menarche (<13 years, 14–15 years, and ≥16 years); number of pregnancies (≥3, 2, ≤1); number of live births (≥2, ≤1); breastfeeding months (>12 months, ≤12 months); and abortion (no, yes). Abbreviations: 5-MTHF, 5-methyltetrahydrofolate; CI, confidence interval; OR, odds ratio; P5P, pyridoxal 5-phosphate; SAM, S-adenosylmethionine; and SAH, S-adenosylhomocysteine.
Table 4. Information on 17 selected differential methylation probes (DMPs).
Table 4. Information on 17 selected differential methylation probes (DMPs).
CpGΔβChrCpG PositionCpG Related TraitsGeneGene PositionGene Molecular FunctionGene Pathways
cg05397629−0.59711opensea/RHOGBodyNucleotide binding; GTPase activityPI5P, PP2A, and IER3 regulate PI3K/AKT Signaling; signaling by Rho GTPases
cg02589685−0.45111shore//IGR//
cg03363289−0.2579islandBreast cancer prognosis; colorectal cancerLHX6BodyDNA binding/
cg04622888−0.2399islandGingivo-buccal oral squamous cell carcinoma; breast cancer prognosis; colorectal cancerLHX6TSS200DNA binding/
cg09491380−0.2344opensea/MAEABodyActin binding; ubiquitin protein transferase activityCiliary landscape
cg10617037−0.2298shore/CYHR11stExonZinc ion binding/
cg15740243−0.2295shore/RNF145TSS1500Zinc ion binding; transferase activity/
cg00901687−0.21817shoreAging; Alzheimer’s disease (AD)MYCBPAPTSS1500Protein binding; phospholipid binding; clathrin binding/
cg24312537−0.2158shoreCrohn’s disease (CD); gestational diabetes mellitus; mortalityHTRA4TSS1500Serine-type endopeptidase activity; protein binding/
cg06536614−0.2155islandBreast cancer riskMIR886TSS200//
cg04657470−0.2132islandTetralogy of FallotHSPE11stExonRNA binding; protein binding; ATP bindingSignaling by Rho GTPases; RAC2 GTPase cycle
cg04794268−0.2088shoreEarly metastasis of uveal melanomaKIFC2TSS1500Nucleotide binding; cytoskeletal motor activityGolgi-to-ER retrograde transport; vesicle-mediated transport
cg21042336−0.20418shore/OSBPL1ATSS1500Protein binding; phospholipid binding; lipid bindingSynthesis of bile acids and bile salts; metabolism
cg21697381−0.20317openseaB acute lymphoblastic leukemiaSLFN12TSS1500RNA nuclease activity; protein binding17q12 copy number variation syndrome
cg01701207−0.20316openseaAllergic asthmaSF3B3BodyNucleic acid binding; protein bindingProcessing of Capped Intron-Containing Pre-mRNA
cg257554280.22919island/MRI1TSS1500Protein binding; isomerase activityMethionine de novo and salvage pathway; sulfur amino acid metabolism
cg138242700.2706shoreBreast cancerPRPF4BTSS1500Nucleotide binding; RNA binding; protein kinase activityProcessing of Capped Intron-Containing Pre-mRNA
The information on CpG-related traits was obtained from EWAS Atlas (ngdc.cncb.ac.cn/ewas/atlas, accessed on 28 July 2023). The information on gene molecular function and gene pathways was obtained from GeneCards (www.genecards.org, accessed on 28 July 2023). Abbreviations: Chr, chromosome; IGR, intergenic region; TSS, transcription start site.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, F.; Zhou, H.; Peng, Y.; Qiao, Y.; Wang, P.; Si, C.; Wang, X.; Gong, J.; Chen, K.; Song, F. Plasma One-Carbon Metabolism-Related Micronutrients and the Risk of Breast Cancer: Involvement of DNA Methylation. Nutrients 2023, 15, 3621. https://doi.org/10.3390/nu15163621

AMA Style

Liu F, Zhou H, Peng Y, Qiao Y, Wang P, Si C, Wang X, Gong J, Chen K, Song F. Plasma One-Carbon Metabolism-Related Micronutrients and the Risk of Breast Cancer: Involvement of DNA Methylation. Nutrients. 2023; 15(16):3621. https://doi.org/10.3390/nu15163621

Chicago/Turabian Style

Liu, Fubin, Huijun Zhou, Yu Peng, Yating Qiao, Peng Wang, Changyu Si, Xixuan Wang, Jianxiao Gong, Kexin Chen, and Fangfang Song. 2023. "Plasma One-Carbon Metabolism-Related Micronutrients and the Risk of Breast Cancer: Involvement of DNA Methylation" Nutrients 15, no. 16: 3621. https://doi.org/10.3390/nu15163621

APA Style

Liu, F., Zhou, H., Peng, Y., Qiao, Y., Wang, P., Si, C., Wang, X., Gong, J., Chen, K., & Song, F. (2023). Plasma One-Carbon Metabolism-Related Micronutrients and the Risk of Breast Cancer: Involvement of DNA Methylation. Nutrients, 15(16), 3621. https://doi.org/10.3390/nu15163621

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop