Diet, Sports, and Psychological Stress as Modulators of Breast Cancer Risk: Focus on OPRM1 Methylation

Background: DNA methylation is influenced by environmental factors and contributes to adverse modification of cancer risk and clinicopathological features.

Methods: A case-control study (402 newly diagnosed cases, 470 controls) was conducted to evaluate the effect of environmental factors and OPRM1 methylation in peripheral blood leukocyte (PBL) DNA on the risk of breast cancer. A case-only study (373 cases) was designed to evaluate the effects of environmental factors on OPRM1 methylation in tumor tissue and the relationship of methylation with clinicopathological features.

Results: We found a significant association between hypermethylation of OPRM1 and the risk of breast cancer (OR = 1.914, 95%CI = 1.357–2.777). OPRM1 hypermethylation in PBL DNA combined with low intake of vegetable, garlic, soybean, poultry, and milk; high pork intake; less regular sports and a high psychological stress index significantly increased the risk of breast cancer. Soybean intake (OR = 0.425, 95%CI: 0.231–0.781) and regular sports (OR = 0.624, 95%CI: 0.399–0.976) were associated with OPRM1 hypermethylation in tumor DNA. OPRM1 hypermethylation in tumor tissue was correlated with estrogen receptor (ER) (OR = 1.945, 95%CI: 1.262–2.996) and progesterone receptor (PR) (OR = 1.611, 95%CI: 1.069–2.427) negative status; in addition, OPRM1 hypermethylation in PBL DNA was associated with human epidermal growth factor receptor 2 (HER-2) negative status (OR = 3.673, 95%CI: 1.411–9.564).

Conclusion: A healthy diet, psychosocial adaptability, and regular sports are very beneficial for breast cancer prevention and progress, especially for OPRM1 hypermethylation carriers. Personalized treatment considering the correlation between OPRM1 hypermethylation and ER and PR status may provide a novel benefit for breast cancer patients.

Introduction

Breast cancer is one of the most common public health issues among women (1). According to an International Agency for Research on Cancer (IARC) report, ~2.1 million new female breast cancer cases and 626,800 cancer deaths occurred worldwide in 2018 (2). The Chinese cancer registry estimated 304,000 new breast cancer patients and 70,000 deaths among women in 2015 (2). Despite significant progress being made in treatment, breast cancer is still the main cancer challenge (3).

Breast cancer is a heterogeneous disease with distinct genetic, epigenetic and histopathological characteristics (4). DNA methylation as an important epigenetic process can regulate gene expression and influence the occurrence of cancers by transferring a methyl group to the C5 position of cytosine (5). Several studies have reported a link of peripheral blood leukocyte (PBL) DNA methylation with cancer diagnosis (6). Ji et al. (7) detected 5-methylcytosine (5-mdc) levels in PBL DNA in breast cancer cases and healthy controls and confirmed that leukocyte genome-wide DNA hypomethylation is independently associated with breast cancer risk.

DNA methylation can be altered by environmental factors, which in turn affect breast cancer risk (8). Tobacco exposure is associated with the methylation of several tumor suppressor genes in breast cancer tissues (9). Alcohol ingestion in animals is known to potentially inhibit folate-mediated methionine synthesis, disrupt key methylation processes and increase the risk of breast cancer (10). Moreover, DNA methylation is associated with specific clinical characteristics that affect the progress of tumors, such as estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER-2) (11). ER and PR are poorly expressed in normal breast cells, but ER alpha (ERα) expression increases at the very earliest stages of tumorigenesis, implying that dysregulation of ERα expression contributes to breast tumor formation (12). HER-2 is often overexpressed in a variety of epithelial-derived tumors and is closely related to the development and prognosis of breast cancer (13).

Opioid receptor mu 1 (OPRM1) gene encodes the mu opioid receptor (MOR), one of three opioid receptors in humans (14). MOR is the main target of opioid analgesics such as endogenous opioid peptides, endorphins, and enkephalins (15). Endogenous opioid systems have been implicated in a wide variety of physiological and pathophysiological actions (16). Partington et al. (17) reported that endogenous opioids are involved in neurocognitive processing of social pain and reward, and OPRM1 SNP minor allele carriers are known to be sensitive to the negative effects of social stress (18). Increased methylation of OPRM1 is associated with Alzheimer's disease, opioid exposure, alcohol and drug addiction, etc. (19–21). Zagon et al. (22) showed that endogenous opioids inhibit tumor growth, while opioid antagonists promote cancer. However, the relationship between OPRM1 methylation and breast cancer risk has not been reported yet.

Therefore, this study was carried out to investigate the relationship of OPRM1 methylation in PBL DNA, environmental factors, and their combinations and interactions with breast cancer risk. We also investigated the effects of environmental factors on OPRM1 methylation in breast tumor DNA and the relationship between OPRM1 methylation and clinicopathological features in tumor DNA and PBL DNA.

Materials and Methods

Study Subjects

A case-control study was designed. A total of 402 breast cancer patients newly diagnosed in the Third Affiliated Clinical Hospital of Harbin Medical University from 2010 to 2014 were enrolled; none of them had received radiation therapy or chemotherapy. In addition, 217 patients from the Department of Orthopedics and Ophthalmology of Harbin Medical University and 253 volunteers from Xiangfang District were assembled as the control group. The control group excluded subjects with a history of breast cancer; subjects with malignant or benign tumors; and pregnant, postpartum, and lactating women. Each participant donated 5 ml of fasting peripheral venous blood before surgery or at registration.

We also conducted a case-only study to analyze the effects of environmental factors on OPRM1 methylation in tumor tissue DNA and the relationship between clinicopathological features and OPRM1 methylation in both PBL and tumor DNA. Tissue specimens of 373 breast cancer patients were collected, frozen in liquid nitrogen as soon as possible, and immediately stored in a −80°C refrigerator.

Ethics Approval and Consent to Participate

All procedures performed in studies involving human participants were granted from the ethical standards of the Human Research and Ethics Committee of Harbin Medical University and the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. All subjects gave informed consent.

Data Collection

All participants were interviewed face-to-face by trained investigators. The content of the questionnaire included demographic characteristics (age, marital status, education level, etc.), behavior (smoking, drinking, regular sports, diet, etc.), menstrual and reproductive history, etc. The clinical and pathological information of cancer patients was extracted from medical records, including the TNM stage, histological records, and other pathological results. We did our best to obtain a complete questionnaire from each participant, but there were still some subjects who did not answer individual questions. Therefore, the n value in the table is not completely corresponding to the sample size.

Genomic DNA Extraction and Sodium Bisulfite Modification

DNA was extracted from blood samples and <25 mg of minced tumor tissues using the QIAamp DNA Blood Mini kit (Qiagen, Hilden, Germany) and PureLinkTM Genomic DNA Kit (Thermo Fisher Scientific, Carlsbad, USA), respectively. All steps were performed in accordance with the manufacturer's protocol, and purified DNA was stored at −80°C as soon as possible. Bisulfite conversion was performed according to the manufacturer's guidelines using 2 μg of DNA and the EpiTect Bisulfite Kit (Qiagen, Hilden, Germany). Before and after the transformation process, the amount of DNA was determined using a Nanodrop 2000 spectrophotometer (Thermo Fisher Scientific, Carlsbad, USA). The transformed DNA was also immediately stored at −80°C.

Plasmid Construction and Standard Curve Dilution

According to the manufacturer's instructions, the TOPO TA Cloning for Sequencing kit (Invitrogen, USA) was used to transfer freshly purified, target-region-specific PCR products to TOP10 competent cells for cloning (23). The extracted plasmid DNA was confirmed by sequencing. Serial dilutions of 10–10⁷ copies of plasmid DNA were used to construct a standard curve for quantification.

Quantification of DNA Methylation With Quantitative Methylation-Specific PCR

DNA methylation was measured by quantitative methylation-specific PCR technology (qMSP). Bisulfite-treated DNA from each specimen served as the template, and qMSP was performed with iTaq Universal SYBR Green Supermix (Bio Rad, America) and ABI 7500 real-time PCR amplifier (Applied Biosystems, Foster City, CA, USA) while strictly following the operating protocols.

Primers were designed for qMSP analysis using Primer Premier 5.0 software as follows: OPRM1: forward primer, 5′-CGGTTATTTATCGTTTGTAGGAGGAAACG-3′ reverse primer; 5′-ATCCAACAACGCTTCTATTCGAATCCG-3′. The internal reference was designed as follows: housekeeping gene (MyOD), forward primer, 5′-CCAACTCCAAATCCCCTCTCTAT-3′; reverse primer, 5′-TGATTAATTTAGATTGGGTTTAGAGAAGGA-3′.

For the amplification, 1 μl bisulfite converted DNA (5 ng) was added to 19 μl amplification Master Mix containing 10 μl SYBR Green Supermix, 0.6 μl of forward and reverse primer, and 7.8 μl RNAse-free water. The reaction protocol was as follows: (1) OPRM1, 6 min at 95°C, followed by 42 cycles of 95°C for 30 s, 58°C for 30 s, and 72°C for 30 s and (2) MyOD, 6 min at 95°C, followed by 42 cycles of 95°C for 30 s, 56°C for 30 s, and 72°C for 30 s. All the procedures were performed in duplicate, and 10–20% of the samples were redetected in our study.

The specific of amplification products were determined with melting curve analysis according to fluorescence data acquired during dissociation steps. Each sample's methylation level was calculated by the following formula: Methylation level = (quantity of target gene/quantity of housekeeping gene) × 100%.

Statistical Analysis

Categorical and continuous variables were tested by a chi-square test and independent-sample t-test, respectively. Univariate and multivariate logistic regression analysis were used to calculate the crude and adjusted odds ratio (ORs) and 95% confidence interval (95%CIs). Combined and interaction effects were performed by crossover analysis (24) and multivariate logistic regression. Correlations between clinicopathological characteristics and the OPRM1 methylation level in tumor DNA and PBL DNA were evaluated using ORs and 95%CIs derived from unconditional logistic regression. The effect of environment factors on OPRM1 methylation in tumor tissue DNA was calculated using unconditional univariate and multivariate logistic regression. All statistical analyses were performed using SPSS version 24.0, with P-values of <0.05 considered statistically significant.

Results

Quality Evaluation of Quantitative Methylation-Specific PCR for OPRM1

We constructed a standard curve using a 10–10⁷ copy number of plasmid DNA as template DNA in PCR process. Figure 1A shows the standard curve of quantitative methylation-specific PCR for OPRM1. The R² value was 0.998, and the amplification efficiency (EFF %) was 53.218%. The signals of samples were compared to the standard curve. Figure 1B displays the amplification plot of serial dilutions of bisulfite modified, universal methylated DNA and samples. All samples were well-amplified, and the amplification plots reached a plateau. The melt curves of bisulfite modified universal methylated DNA and samples were also compared (Figure 1C).

FIGURE 1

Figure 1. Quantitative methylation-specific PCR (qMSP) curves of OPRM1 methylation. (A) Standard curve. The R² value was 0.998, and the amplificationefficiency (EFF %) was 53.218%. (B) Amplification plot. (C) Melt curve. The melting peaks of serial dilutions of plasmid DNA and samples.

The Relationship Between OPRM1 Methylation in PBL DNA and Breast Cancer Risk

A total of 402 PBL DNA samples from breast cancer patients and 470 control samples were analyzed. The distribution of basic demographic characteristics is shown in Supplementary Table 1. The differences were statistically significant in the categories of marital status (P = 0.040), educational level (P = 0.001), occupation (P = 0.000), and family history of cancer (P = 0.000). Therefore, these characteristics were adjusted as confounding factors in the following analyses. As shown in Table 1, by adjusting the marital status, educational level, occupation, and family history of cancer, we found a significant association between the hypermethylation of OPRM1 in PBL DNA and the risk of breast cancer (OR = 1.914, 95%CI = 1.357–2.777, P = 0.000).

TABLE 1

Table 1. Association of OPRM1 methylation in peripheral blood leukocyte DNA and breast cancer risk.

Relationship of Environmental Exposure and Breast Cancer Risk

Univariate and multivariate logistic regression were used to analyze the relationship between environmental factors and breast cancer risk. As shown in Table 2, according to univariate logistic regression analysis, factors such as pork intake and psychological stress index are risk factors for breast cancer (P < 0.05). Intake of vegetable, garlic, poultry, milk, and soybean and regular sports are protective factors against breast cancer (P < 0.05). Subsequently, those environmental factors identified as significant were examined in the multivariate analyses. The results showed that vegetable intake (≥500 vs. <500 g/d, OR = 0.576, 95%CI: 0.407–0.816), garlic intake (≥1 vs. <1 time/week, OR = 0.468, 95%CI: 0.331–0.662), poultry intake (≥1 vs. <1 time/week, OR = 0.411, 95%CI: 0.281–0.601), milk intake (≥3 vs. <3 times/week, OR = 0.496, 95%CI: 0.335–0.732), regular sports (≥1 vs. <1 time/week OR = 0.604, 95%CI: 0.419–0.871), and soybean intake (≥1 vs. <1 time/week, OR = 0.500, CI: 0.291–0.859) were significantly associated with a decreased risk of breast cancer. A high intake of pork (≥1 vs. <1 time/week, OR = 2.956, 95%CI: 1.792–4.875) and a high psychological stress index (>15 vs. ≤15, OR = 1.938, 95%CI: 1.293–2.903) were significantly associated with an increased risk of breast cancer.

TABLE 2

Table 2. Univariate and multivariate analysis for the association of environmental factors and the risk of breast cancer.

Combined and Interactive Effects of OPRM1 Methylation and Environmental Factors in Breast Cancer

As shown in Table 3, comparing with the reference group, statistically significant combinations of OPRM1 hypermethylation and low vegetable intake (<500 g/day, OR = 2.829, 95%CI: 1.751–4.572, P = 0.000), garlic intake (<1 time/week, OR = 5.762, 95%CI: 3.208–10.350, P = 0.000), soybean intake (<1 time/week, OR = 5.592, 95%CI: 2.905–10.764, P = 0.000), poultry intake (<1 time/week, OR = 4.173, 95%CI: 2.415–7.210, P = 0.002), milk intake (<3 times/week, OR = 5.166, 95%CI: 2.774–9.618, P = 0.000), and less sports (<1 time/week, OR = 4.114, 95%CI: 2.357–7.181, P = 0.000) were observed to be significantly related to increased breast cancer risk. In addition, we also found a significant combination of OPRM1 hypermethylation and a high pork intake (≥1 time/week, OR = 3.628, 95%CI: 1.593–8.262, P = 0.000) and high psychological stress index (>15, OR = 5.340, 95%CI: 2.634–10.829, P = 0.000) for increased breast cancer risk. There was no significant interaction between OPRM1 hypermethylation and any environmental factor in breast cancer.

TABLE 3

Table 3. Combined and interactive effects between OPRM1 methylation and environmental factors in breast cancer.

Effects of Exposure to Environmental Factors on OPRM1 Methylation in Tumor Tissue DNA

A total of 373 patients with breast cancer were included in this analysis. The basic demographic characteristics for the hypermethylated and hypomethylated groups are shown in Supplementary Table 2. No significant difference in age, educational level, marital status, occupation, family history of cancer, or BMI distribution (all P > 0.05) between the OPRM1 hypermethylated group and the hypomethylated group was observed.

Based on the results of univariate logistic regression, the history of estrogen therapy (yes vs. no, OR = 2.712, 95%CI: 1.163–6.325), soybean intake (≥1 vs. <1 time/week, OR = 0.420, 95%CI: 0.233–0.758), and regular sports (≥1 vs. <1 time/week, OR = 0.606, 95%CI: 0.392–0.938) are significantly related with OPRM1 hypermethylation in tumor tissue DNA. After introducing these variables into a multivariate model, the soybean intake (≥1 vs. <1 time/week, OR = 0.425, 95%CI: 0.231–0.781) and regular sports (≥1 vs. <1 time/week, OR = 0.624, 95%CI: 0.399–0.976) were still significantly associated with OPRM1 hypermethylation (Table 4).

TABLE 4

Table 4. Relationship between environmental factors exposures and OPRM1 methylation in breast cancer tissue DNA.

Relationship Between Clinicopathological Characteristics and OPRM1 Methylation in PBL and Tumor Tissue DNA

The relationship between clinicopathological characteristics and OPRM1 methylation in breast tumor tissue DNA and PBL DNA is shown in Table 5. We found that OPRM1 hypermethylation in tumor tissue DNA was significantly correlated with ER status (negative vs. positive, OR = 1.945, 95%CI: 1.262–2.996) and PR status (negative vs. positive, OR = 1.611, 95%CI: 1.06–2.427). In addition, OPRM1 hypermethylation in PBL DNA was significantly associated with HER-2 negative status (vs. positive, OR = 3.673, 95%CI: 1.411–9.564). No significant relation between OPRM1 hypermethylation in either tumor tissue DNA or PBL DNA with the TNM stage, pathological type, tumor invasion status, lymphnode involved, metastasis status, histological type or P53 expression was found.

TABLE 5

Table 5. Correlation between clinicopathological characteristics and OPRM1 methylation in breast tumor tissue DNA and peripheral blood leukocyte DNA.

Discussion

Accumulating research has shown the importance of DNA methylation in the impact of environmental factors on breast cancer risk, and this epigenetic mechanism is also associated with specific clinicopathological features (25). In this study, we discovered for the first time a series of valuable associations between OPRM1 methylation and breast cancer. First, we found that OPRM1 hypermethylation in PBL DNA increased the breast cancer risk. Furthermore, OPRM1 hypermethylation combined with low intake of vegetable, garlic, soybean, poultry, and milk; less regular sports; and a high pork intake and high psychological stress index significantly increased the risk of breast cancer. We also found that soybean intake and regular sports reduced OPRM1 methylation in tumor tissue. Additionally, OPRM1 hypermethylation in tumor tissue DNA was significantly correlated with ER and PR negative status, and hypermethylation in PBL DNA was significantly associated with HER-2 negative status.

Some studies have indicated an association between psychological stress and breast cancer risk, which is consistent with our results (26). Our study found that a high psychological stress index (>15) increased the risk of breast cancer (OR = 1.938, 95%CI: 1.293–2.903). Stress hormone signaling may induce DNA damage and promote tumorigenesis by producing reactive oxygen species (ROS)/reactive nitrogen species (RNS) and interfering with DNA repair processes (27). Additionally, we found that high psychological stress combined with OPRM1 hypermethylation in PBL DNA increased the breast cancer risk with a strong correlation (OR = 5.340, 95%CI: 2.634–10.829). Guerrero-Alba et al. (28) reported that psychological stress factors can significantly reduce endogenous opioid levels, block the binding of endogenous opioids to receptors, and thus lead to harmful events. Psychological stress may be correlated with hypermethylation of OPRM1 and low expression of OPRM1-encoded mu-opioid receptors, which promote the occurrence of breast cancer. Thus, we suggest that individuals with OPRM1 hypermethylation in PBL DNA maintain a good psychological state to prevent breast cancer.

Steensberg et al. (29) reported that increased exercise contributes to the upregulation of IL-6, which enhances the levels of anti-inflammatory cytokines IL-1α and IL-10. In addition, IL-6 inhibits the expression of tumor necrosis factor-α. According to our results, regular sports can reduce the risk of breast cancer significantly (OR = 0.604, 95%CI: 0.419–0.871). In particular, comparing with the individuals of OPRM1 hypomethylation in PBL DNA and regular sports, we found that less frequent sports (<1 time/week) combined with OPRM1 hypermethylation in PBL DNA increased breast cancer risk with a strong association (OR = 4.114, 95%CI: 2.357–7.181). Arida et al. (30) randomly divided rats into acute (7 days) and chronic (30 days) exercise groups and a control group and observed higher opioid receptor binding in the acute-exercise animals; they concluded that physical exercise can stimulate the release of endogenous opioid substances. By inference, less sports and hypermethylation of OPRM1 may have combined action on the inhibition of endogenous opioid substances and increased the risk of breast cancer. Thus, we encourage people with OPRM1 hypermethylation in PBL DNA to increase their physical activity to stay healthy. Strikingly, we also found that regular sport is associated with a decreased OPRM1 methylation level in breast cancer tissues (OR = 0.624, 95%CI: 0.399–0.976). All of these results suggest that regular sports can not only reduce the risk of breast cancer directly but also indirectly change the individual's susceptibility to breast cancer through modification of DNA methylation.

In addition, many dietary factors are also related to breast cancer risk (31). We found that vegetable, garlic, poultry, and milk intake are protection factors and that pork intake is a risk factor for breast cancer (P < 0.05). Vegetables are rich in dietary fiber, which may decrease reabsorption in the gut of estrogen excreted in the biliary system (32). Garlic contains organosulfur compounds, of which diallyl trisulfide inhibits ER-α activity in human breast cancer cells (33). Milk is rich in vitamin D; Xie et al. (34) found that vitamin D analogs suppress IGF-I signaling and promote apoptosis in breast cancer cells. The protective effect of poultry can be explained by its amino acid content; for example, n-3 fatty acids inhibit tumor growth and metastasis by the effect of impaired angiogenesis (35). However, pork intake may increase breast cancer risk due to dietary heme iron, fat and N-glycolylneuraminic acid, which are indicated to possibly increase tumor formation (36). A combination of low vegetable, garlic, poultry, and milk intake and high pork intake and OPRM1 hypermethylation in PBL DNA was observed in our study to be significantly associated with an increased risk of breast cancer (P < 0.05). Therefore, we also encourage that women with OPRM1 hypermethylation should increase their intake of vegetable, garlic, milk and poultry and reduce their intake of pork to prevent breast cancer.

Genistein has been identified as the predominant isoflavone in soybean, which can prevent abundant cell proliferation or abnormal angiogenesis by inhibiting RTK-mediated signaling pathways (37). Zhu et al. (38) found higher consumption of soy protein decreased breast cancer risk significantly (OR = 0.46, 95%CI = 0.24–0.88). And the Shanghai Breast Cancer Survival Study, a large, population-based cohort study of 5,042 female breast cancer survivors in China, found that soy food consumption was significantly associated with decreased risk of death (HR = 0.71, 95%CI, 0.54–0.92) and recurrence (HR = 0.68,95%CI, 0.54–0.87) (39). We found that soybean intake reduced the risk of breast cancer (OR = 0.500, 95%CI: 0.291–0.859) and that low soybean intake combined with OPRM1 hypermethylation increased the risk of breast cancer (OR = 5.592, 95%CI: 2.905–10.764). Moreover, soybean intake reduced OPRM1 methylation in tumor tissue (OR = 0.425, 95%CI: 0.231–0.781). However, Katrin Sak reviewed the existing literature and reported that intake of isoflavones can be associated with a decrease in breast tumorigenesis only in Asian countries where the consumption of soy foods is high but not among Western women with significantly lower ingestion amounts (40). Some other studies revealed that the incidence of hormone-dependent breast cancer is much lower in the Asian population than in the United States and Europe (41). Soybean consumption may be responsible in part for lower levels of hormones and decreased rates of breast cancer in women in Asia compared with Western populations (42). Therefore, as for the inconsistent findings regarding soybeans and breast cancer risk, larger epidemiological studies involving multiracial populations are still needed to confirm the effect of genes and soy on breast cancer risk.

A clinicopathology staging system including ER, PR, and HER-2 was previously validated for the treatment and prognostic value of breast cancer in patients with breast cancer (43). Zagon et al. (44), Gach et al. (45), and Hatzoglou et al. (46) reported that the expression of OPRM1 is related to ER positivity, which leads to ERα activation through MOR activation and mediates its translocation to plasma membrane, synergistically inducing the proliferation of breast cancer cells. OPRM1 hypermethylation leads to a decrease in its gene expression, which also weakens ER activity. Similarly, in our study, we found that OPRM1 hypermethylation was significantly associated with ER negative status and PR negative status in tumor tissue DNA (P < 0.05). As Dunnwald et al. (47) reported, ER- and PR-positive status in tumors not only improves the therapeutic effect but also increases the survival rate, so directional therapy plays a vital role in adjuvant breast cancer treatment. Although the mechanism of the association between OPRM1 hypermethylation and ER/PR-negative status remains to be investigated, as does the association between OPRM1 hypermethylation in PBL DNA and HER-2 negative status, our results may provide new ideas for the treatment of breast cancer patients through further studies on epigenetic regulation. Targeted induction of OPRM1 methylation to regulate gene expression and design of possible therapeutic drugs may provide a novel remedy for breast cancer patients.

There are still some limitations with the interpretation of the broad results in our research. First, recall bias when collecting environmental factor data is difficult to avoid. In addition, all participants answered the same questionnaire, but different questions were omitted by different participants, so that the n value in the table could not correspond to the sample size. Second, the information about dietary factors in the questionnaire was divided into only two different frequencies, and the analysis of exposure dose was not detailed enough to give quantitative results. The sample size adopted in our subgroup analysis was relatively small, and the results need to be verified by a study with a larger sample size. The relationship between OPRM1 methylation and increased breast cancer risk is based on the results of population epidemiological studies. The mechanism of OPRM1 methylation and tumorigenesis based on in vitro studies still needs to be further developed and verified. It is also necessary to carry out in vivo interventional experiment in OPRM1 hypermethylation population to verify the effect of individualized prevention.

Conclusion

This study identified that OPRM1 hypermethylation in PBL DNA is correlated with increased risk of breast cancer. A healthy diet with sufficient intake of vegetable, garlic, soybean, poultry and milk; a low pork intake; regular sports and healthy psychosocial adaptability is very beneficial for breast cancer prevention, especially for OPRM1 hypermethylation carriers. Personalized promotion on regular sports and soybean intake and precision treatment based on DNA methylation markers (OPRM1) should be encouraged with the consideration of the correlation between OPRM1 hypermethylation and ER/PR, HER-2 negative status in breast cancer.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics Statement

The studies involving human participants were reviewed and approved by Human Research and Ethics Committee of Harbin Medical University. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

FW contributed to the study design. LL and SL contributed to methylation detection, data analysis, and drafting the manuscript. NZ contributed to data analysis. SQ, JD, YC, and CW contributed to data collection. YG, ZH, and YY contributed to data interpretation. YZ contributed to manuscript revision. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from National Nature Science Foundation of China (Nos. 81202262 and 81773503), Postdoctoral Science Special Foundation of China (2013T60390), and Dr. Wu Lien-teh Science Foundation of Harbin Medical University (WLD-QN1106).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2021.747964/full#supplementary-material

References

1. Liu FC, Lin HT, Kuo CF, See LC, Chiou MJ, Yu HP. Epidemiology and survival outcome of breast cancer in a nationwide study. Oncotarget. (2017) 8:16939–50. doi: 10.18632/oncotarget.15207

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Zheng RS, Sun KX, Zhang SW, Zeng HM, Zou XN, Chen R, et al. Report of cancer epidemiology in China, 2015. Zhonghua Zhong Liu Za Zhi. (2019) 41:19–28. doi: 10.3760/cma.j.issn.0253-3766.2019.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Cao M, Li H, Sun D, Chen W. Cancer burden of major cancers in China: a need for sustainable actions. Cancer Commun. (2020) 40:205–10. doi: 10.1002/cac2.12025

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Klajic J, Fleischer T, Dejeux E, Edvardsen H, Warnberg F, Bukholm I, et al. Quantitative DNA methylation analyses reveal stage dependent DNA methylation and association to clinico-pathological factors in breast tumors. BMC Cancer. (2013) 13:456. doi: 10.1186/1471-2407-13-456

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Jjingo D, Conley AB, Yi SV, Lunyak VV, Jordan IK. On the presence and role of human gene-body DNA methylation. Oncotarget. (2012) 3:462–74. doi: 10.18632/oncotarget.497

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Marsit CJ, Koestler DC, Christensen BC, Karagas MR, Houseman EA, Kelsey KT. DNA methylation array analysis identifies profiles of blood-derived DNA methylation associated with bladder cancer. J Clin Oncol. (2011) 29:1133–9. doi: 10.1200/JCO.2010.31.3577

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Choi JY, James SR, Link PA, McCann SE, Hong CC, Davis W, et al. Association between global DNA hypomethylation in leukocytes and risk of breast cancer. Carcinogenesis. (2009) 30:1889–97. doi: 10.1093/carcin/bgp143

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Skinner MK. Environmental epigenetics and a unified theory of the molecular aspects of evolution: a neo-lamarckian concept that facilitates neo-darwinian evolution. Genome Biol Evol. (2015) 7:1296–302. doi: 10.1093/gbe/evv073

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Callahan CL, Bonner MR, Nie J, Han D, Wang Y, Tao MH, et al. Lifetime exposure to ambient air pollution and methylation of tumor suppressor genes in breast tumors. Environ Res. (2018) 161:418–24. doi: 10.1016/j.envres.2017.11.040

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Mason JB, Choi SW. Effects of alcohol on folate metabolism: implications for carcinogenesis. Alcohol. (2005) 35:235–41. doi: 10.1016/j.alcohol.2005.03.012

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Shakeri H, Fakhrjou A, Nikanfar A, Mohaddes-Ardebili SM. Methylation analysis of BRCA1 and APC in breast cancer and it's relationship to clinicopathological features. Clin Lab. (2016) 62:2333–7. doi: 10.7754/Clin.Lab.2016.160418

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Anderson E. The role of oestrogen and progesterone receptors in human mammary development and tumorigenesis. Breast Cancer Res. (2002) 4:197–201. doi: 10.1186/bcr452

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Tsutsui S, Ohno S, Murakami S, Kataoka A, Kinoshita J, Hachitanda Y. Prognostic significance of the coexpression of p53 protein and c-erbB2 in breast cancer. Am J Surg. (2003) 185:165–7. doi: 10.1016/S0002-9610(02)01203-5

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Crist RC, Berrettini WH. Pharmacogenetics of OPRM1. Pharmacol Biochem Behav. (2014) 123:25–33. doi: 10.1016/j.pbb.2013.10.018

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Oladosu FA, Conrad MS, O'Buckley SC, Rashid NU, Slade GD, Nackley AG. Mu opioid splice variant MOR-1K contributes to the development of opioid-induced hyperalgesia. PLoS ONE. (2015) 10:e0135711. doi: 10.1371/journal.pone.0135711

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Beaumont A, Hughes J. Biology of opioid peptides. Annu Rev Pharmacol Toxicol. (1979) 19:245–67. doi: 10.1146/annurev.pa.19.040179.001333

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Partington LC, Borelli JL, Smiley P, Jarvik E, Rasmussen HF, Seaman LC, et al. Parental overcontrol x OPRM1 genotype interaction predicts school-aged children's sympathetic nervous system activation in response to performance challenge. Res Dev Disabil. (2018) 82:39–52. doi: 10.1016/j.ridd.2018.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Swann G, Byck GR, Dick DM, Aliev F, Latendresse SJ, Riley B, et al. Effect of OPRM1 and stressful life events on symptoms of major depression in African American adolescents. J Affect Disord. (2014) 162:12–9. doi: 10.1016/j.jad.2014.03.020

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Xu C, Liu G, Ji H, Chen W, Dai D, Chen Z, et al. Elevated methylation of OPRM1 and OPRL1 genes in Alzheimer's disease. Mol Med Rep. (2018) 18:4297–302. doi: 10.3892/mmr.2018.9424

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Viet CT, Dang D, Aouizerat BE, Miaskowski C, Ye Y, Viet DT, et al. OPRM1 methylation contributes to opioid tolerance in cancer patients. J Pain. (2017) 18:1046–59. doi: 10.1016/j.jpain.2017.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Zhang H, Herman AI, Kranzler HR, Anton RF, Simen AA, Gelernter J. Hypermethylation of OPRM1 promoter region in European Americans with alcohol dependence. J Hum genet. (2012) 57:670–5. doi: 10.1038/jhg.2012.98

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Zagon IS, McLaughlin PJ, Takemori AE, Portoghese PS. beta-Funaltrexamine (beta-FNA) and neural tumor response in mice. Eur J Pharmacol. (1985) 116:165–9. doi: 10.1016/0014-2999(85)90198-0

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Ramalho-Carvalho J, Rui H, Jerónimo C. Methylation-specific PCR. Methods Mol Biol. (2018) 1708:447–72. doi: 10.1007/978-1-4939-7481-8_23

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Wu YZ, Yang H, Zhang L, Zhang YQ, Liu L, Yi D, et al. Application of crossover analysis-logistic regression in the assessment of gene- environmental interactions for colorectal cancer. Asian Pac J Cancer Prev. (2012) 13:2031–7. doi: 10.7314/APJCP.2012.13.5.2031

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Li J, Chen Y, Yu H, Tian J, Yuan F, Fan J, et al. DUSP1 promoter methylation in peripheral blood leukocyte is associated with triple-negative breast cancer risk. Sci Rep. (2017) 7:43011. doi: 10.1038/srep43011

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Chiriac VF, Baban A, Dumitrascu DL. Psychological stress and breast cancer incidence: a systematic review. Clujul Medical. (2018) 91:18–26. doi: 10.15386/cjmed-924

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Flaherty RL, Owen M, Fagan-Murphy A, Intabli H, Healy D, Patel A, et al. Glucocorticoids induce production of reactive oxygen species/reactive nitrogen species and DNA damage through an iNOS mediated pathway in breast cancer. Breast Cancer Res. (2017) 19:35. doi: 10.1186/s13058-017-0823-8

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Guerrero-Alba R, Valdez-Morales EE, Jimenez-Vargas NN, Lopez-Lopez C, Jaramillo-Polanco J, Okamoto T, et al. Stress activates pronociceptive endogenous opioid signalling in DRG neurons during chronic colitis. Gut. (2017) 66:2121–31. doi: 10.1136/gutjnl-2016-311456

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Steensberg A, Fischer CP, Keller C, Møller K, Pedersen BK. IL-6 enhances plasma IL-1ra, IL-10, and cortisol in humans. Am J Physiol Endocrinol Metab. (2003) 285:E433–7. doi: 10.1152/ajpendo.00074.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Arida RM, Gomes da Silva S, de Almeida AA, Cavalheiro EA, Zavala-Tecuapetla C, Brand S, et al. Differential effects of exercise on brain opioid receptor binding and activation in rats. J Neurochem. (2015) 132:206–17. doi: 10.1111/jnc.12976

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Seiler A, Chen MA, Brown RL, Fagundes CP. Obesity, dietary factors, nutrition, and breast cancer risk. Curr Breast Cancer Rep. (2018) 10:14–27. doi: 10.1007/s12609-018-0264-0

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Li Y, Li S, Meng X, Gan RY, Zhang JJ, Li HB. Dietary natural products for prevention and treatment of breast cancer. Nutrients. (2017) 9:728. doi: 10.3390/nu9070728

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Hahm ER, Singh SV. Diallyl trisulfide inhibits estrogen receptor-alpha activity in human breast cancer cells. Breast Cancer Res Treat. (2014) 144:47–57. doi: 10.1007/s10549-014-2841-x

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Xie SP, Pirianov G, Colston KW. Vitamin D analogues suppress IGF-I signalling and promote apoptosis in breast cancer cells. Eur J Cancer. (1999) 35:1717–23. doi: 10.1016/S0959-8049(99)00200-2

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Rose DP, Connolly JM. Regulation of tumor angiogenesis by dietary fatty acids and eicosanoids. Nutr Cancer. (2000) 37:119–27. doi: 10.1207/S15327914NC372_1

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Kim JH, Lee J, Jung SY, Kim J. Dietary factors and female breast cancer risk: a prospective cohort study. Nutrients. (2017) 9:1331. doi: 10.3390/nu9121331

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Varinska L, Gal P, Mojzisova G, Mirossay L, Mojzis J. Soy and breast cancer: focus on angiogenesis. Int J Mol Sci. (2015) 16:11728–49. doi: 10.3390/ijms160511728

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Zhu YY, Zhou L, Jiao SC, Xu LZ. Relationship between soy food intake and breast cancer in China. Asian Pac J Cancer Prev. (2011) 12:2837–40.

PubMed Abstract | Google Scholar

39. Shu XO, Zheng Y, Cai H, Gu K, Chen Z, Zheng W, et al. Soy food intake and breast cancer survival. J Am Med Assoc. (2009) 302:2437–43. doi: 10.1001/jama.2009.1783

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Sak K. Epidemiological evidences on dietary flavonoids and breast cancer risk: a narrative review. Asian Pac J Cancer Prev. (2017) 18:2309–28. doi: 10.22034/APJCP.2017.18.9.2309

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Barnes S. The biochemistry, chemistry and physiology of the isoflavones in soybeans and their food products. Lymphat Res Biol. (2010) 8:89–98. doi: 10.1089/lrb.2009.0030

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Lu LJ, Anderson KE, Grady JJ, Kohen F, Nagamani M. Decreased ovarian hormones during a soya diet: implications for breast cancer prevention. Cancer Res. (2000) 60:4112–21.

PubMed Abstract | Google Scholar

43. Zhang J, Zhao B, Jin F. The assessment of 8th edition AJCC prognostic staging system and a simplified staging system for breast cancer: the analytic results from the SEER database. Breast J. (2019) 25:838–47. doi: 10.1111/tbj.13347

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Zagon IS, McLaughlin PJ, Goodman SR, Rhodes RE. Opioid receptors and endogenous opioids in diverse human and animal cancers. J Natl Cancer Inst. (1987) 79:1059–65.

PubMed Abstract | Google Scholar

45. Gach K, Piestrzeniewicz M, Fichna J, Stefanska B, Szemraj J, Janecka A. Opioid-induced regulation of mu-opioid receptor gene expression in the MCF-7 breast cancer cell line. Biochem. Cell Biol. (2008) 86:217–26. doi: 10.1139/O08-001

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Hatzoglou A, Bakogeorgou E, Castanas E. The antiproliferative effect of opioid receptor agonists on the T47D human breast cancer cell line, is partially mediated through opioid receptors. Eur J Pharmacol. (1996) 296:199–207. doi: 10.1016/0014-2999(95)00703-2

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Dunnwald LK, Rossing MA, Li CI. Hormone receptor status, tumor characteristics, and prognosis: a prospective cohort of breast cancer patients. Breast Cancer Res. (2007) 9:R6. doi: 10.1186/bcr1639

PubMed Abstract | CrossRef Full Text | Google Scholar