- Open Access
Dietary protein sources and tumoral overexpression of RhoA, VEGF-A and VEGFR2 genes among breast cancer patients
Genes & Nutrition volume 14, Article number: 22 (2019)
High protein intake may promote angiogenesis giving support to the development of metastasis according to the experimental data. However, nutritional epidemiologic evidence is inconsistent with metastasis. Therefore, we aimed to study the association between dietary intake of protein and tumoral expression levels of Ras homologous gene family member A (RhoA), vascular endothelial growth factor-A (VEGF-A), and VEGF receptor-2 (VEGFR2) in primary breast cancer (BC) patients.
Over this consecutive case series, 177 women primary diagnosed with histopathologically confirmed BC in Tabriz (Iran) were enrolled between May 2011 and November 2016. A validated food frequency questionnaire was completed for eligible participants. Fold change in gene expression was measured using quantitative real-time PCR. Principal component factor analysis (PCA) was used to express dietary groups of proteins.
Total protein intake was associated with the expression level of VEGF-A in progesterone receptor-positive (PR+: β = 0.296, p < 0.01) and VEGFR2 in patients with involvement of axillary lymph node metastasis (ALNM+: β = 0.295, p < 0.01) when covariates were adjusted. High animal protein intake was correlated with overexpression of RhoA in tumors with estrogen receptor-positive (ER+: β = 0.230, p < 0.05), ALNM+ (β = 0.238, p < 0.05), and vascular invasion (VI+: β = 0.313, p < 0.01). Animal protein intake was correlated with the overexpression of VEGFR2 when tumors were positive for hormonal receptors (ER+: β = 0.299, p < 0.01; PR+: β = 0.296, p < 0.01). Based on the PCA outputs, protein provided by whole meat (white and red meat) was associated inversely with RhoA expression in ALNM+ (β = − 0.253, p < 0.05) and premenopausal women (β = − 0.285, p < 0.01) in adjusted models. Whole meat was correlated with VEGFR2 overexpression in VI+ (β = 0.288, p < 0.05) and premenopausal status (β = 0.300, p < 0.05) in adjusted models. A group composed of dairy products and legumes was correlated with the overexpression of RhoA (β = 0.249, p < 0.05) and VEGF-A (β = 0.297, p < 0.05) in VI+.
Based on the multivariate findings, the dietary protein could associate with the overexpression of RhoA and VEGF-VEGFR2 in favor of lymphatic and vascular metastasis in BC patients.
Breast cancer (BC) is the most frequently diagnosed malignancy in females in many countries . In the recent decade, the prevalence rate of BC has been increasing rapidly among Iranian women . Even lately, BC is the main leading cause of female cancer-related mortality . Metastasis is a critical event in malignancy addressed as an indicator of poor prognosis and causes the vast majority of cancer-related deaths .
Sufficient epidemiologic evidence has revealed the association between some lifestyle-related risk factors and the BC risk . However, it is far less understood how dietary factors can take part in cancer progression toward the formation of metastasis . It is widely studied that breast carcinogenesis seems likely attributable to the high dietary intake of total protein, especially from animal sources [5,6,7]. Lately, in a meta-analysis of prospective studies, it has been documented that higher protein intake from red and processed meat is a potential risk factor for BC . Moreover, in a large prospective cohort study, Cho et al.  found out that red meat intake strongly can elevate BC risk in pre-menopause women who had tumors with hormone receptor-positive status. In Nurses’ Health Study II during 20 years follow-up, it has been reported that high red meat intake in early adulthood might increase the risk of BC in later life . Moreover, despite the previous findings suggested an association between protein intake and increased risk of metastasis [9, 10], a few experimental studies showed the effects of high protein content in association with augmented molecular alterations in promoting metastasis .
A review has highlighted that signaling pathways leading to cytoskeletal reprogramming are vital for cancerous cellular motility . Rho (Ras homologous) is a superfamily of small GTPase, involved Ras homologous gene family member A (RhoA) as a key element function in neoplastic invasion and controlling cellular morphogenesis [13, 14]. Many in vitro experimental studies provide insights into the active contribution of upregulated RhoA in neoplastic propagation through holding the rearrangement of cytoskeletons, cellular motility, and subsequently cancer invasion . On the one hand, Rho proteins modulate both F-actin formation and myosin activation, through RhoA-Rho kinase (ROCK) signaling pathway turned out as a molecular switch to catalyze GTP-GDP exchange . The configuration of active GTP bound allows Rho protein to regulate signal transduction . The phosphatidylinositol-3 kinase (PI3K) activates the RhoA-ROCK signaling pathway . On the other hand, hypoxia-inducible factor 1 (HIF-1)-dependent RhoA-ROCK signaling might result in enhanced BC cellular motility thereby increasing the risk of invasion and metastasis . It is well-established that RhoA overexpressed in breast tumors . However, a few clinical trials indicated that dietary protein intake can affect RhoA expression. Hebels et al.  showed that dietary intervention by red meat (7 days) in patients with irritable bowel disease can increase the expression levels of RhoA in colon tissue. High intake of animal proteins could increase the acid load of the blood, whereby low pH might be a stimulus for the formation of stress fiber mediated by RhoA . To our knowledge, epidemiologic study to support the association of protein intake (the type of protein) and the transcription status of RhoA in breast tumors does not exist.
Angiogenesis is the formation of new blood vessels from the former vasculature, which is crucial in the propagation of tumoral cells to be enabled to grow in other anatomic sites and forming distant metastasis [19, 20]. Silent angiogenesis is precisely regulated by pro- and anti-angiogenic factors in healthy tissue, whereas pathologic angiogenesis potentiates tumorigenesis linked to the distorted balance of angiogenic variables grow endothelial cells . Vascular endothelial growth factor (VEGF) is essential for growth and survival of endothelial cells accounted for promoting pathological angiogenesis [21, 22]. The family VEGF proteins are the major regulators of lymphatic and blood vessel formations . While VEGF-C and VEGF-D considered being involved most often in lymphangiogenesis, VEGF-A is the predominant member of the VEGF family that regulates vasculogenesis essential in both physiologic and pathological angiogenesis . VEGF-A could also induce lymphangiogenesis . VEGF receptor-2 (VEGFR2), a tyrosine kinase receptor, is distributed in all endothelial cells as well as cells of lymphatics capillaries and mediates the signaling pathway activated by VEGF-A as a ligand [23, 24]. VEGF-A is secreted by normal and malignant cells, and its overexpression was reported evidently in breast tumors . Overregulated VEGF-A is significantly associated with poor prognosis in BC . Hypoxia plays a major role in the overexpression of VEGF in breast tumors, because of the binding site on the promoter of VEGF gene for HIF-1 . In the absence of HIF-1, PI3K/Rho/ROCK/c-MYC pathway overrides the effect of VEGF by regulating a cis-regulatory element located in the VEGF gene promoter . Much evidence suggested that high animal protein intake can raise the circulating levels of insulin-like growth factor-1 (IGF-1) that plays an important role in breast tumor progression [10, 26]. IGF-1-dependent PI3K/AKT/mammalian target of rapamycin (mTOR) signaling is a critical transcriptional activator pathway conceivably accounted for VEGF upregulation [10, 26]. No published data indicated whether the dietary protein can associate with the alterations in transcription of VEGF-A and VEGFR2 as a nutrigenomic model may interact in lymphangiogenesis.
Given the fact that overexpression of RhoA and VEGF-VEGFR2 is contributed to poor prognosis in malignancies [13, 22, 27] and evidence for linking protein intake and BC risk is often rare in epidemiologic studies , there is a need to investigate the nutrigenomic aspect of angiogenesis in a population of BC patients. Therefore, we aimed to study the association between dietary sources of protein and RhoA, VEGF-A, and VEGFR2 expression levels in primary BC patients.
Materials and methods
We conducted a consecutive case series study, and 177 eligible newly diagnosed BC patients within the age range of 19–73 years old were recruited from May 2011 to November 2016 at Noor-Nejat Private Hospital in Tabriz, northwest of Iran. The inclusion criteria were histopathological confirmation of primary BC, no chronic inflammation, not vegetarian, not pregnant or a breastfed mother, and no longitudinal usage of drugs (e.g., methotrexate, sulfasalazine, anticonvulsants, and contraceptive drugs). The exclusion criteria consisted of a history of confirmed malignancy during life, history of adjuvant therapy, and history of breast benign disease. Also, those who had a body mass index above 45 kg/m2 were excluded from the study. Individuals’ pathological data were recorded from their medical history including histological tumor grade, tumor size, and histopathological subtype (ductal and non-ductal) and immunohistochemistry data such as estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor-2 (HER2-neu).
Face-to-face interview with each participant was conducted by well-trained interviewers. Food frequency questionnaire (FFQ) with 136 food items was used to assess the dietary intake of protein. The accuracy of FFQ-based dietary intake measures of cobalamin and its association with a dietary group of protein had been met in previous studies [28,29,30]. Food intake of participants was questioned on a timely basis of daily, weekly, monthly, and yearly. FFQ was completed for the previous year. The amount of average daily intake of a nutrient was estimated by multiplying the time-dependent frequency of food intake to the amount of consumption in grams. Portion size was defined based on the common household utensil and then converted to gram using a standard reference value . A set of photographs showing the major food groups was on hand to express better the consumption magnitude for the food items. Nutritionist IV software (ver.3.5.2; 1994, N-squared computing, San Bruno, USA) was used to compute the intake levels of total calories, macronutrients, and fiber. Total calorie intake was adjusted in terms of nutrient density as described by Willett .
Extraction of mRNA and quantitative real-time reverse transcriptase polymerase chain reaction
Fresh tissues (tumor and adjacent normal tissues) were collected after dissection and carried in liquid nitrogen to store at − 70 °C. Based on microscopic examination, tumor cells on average constituted > 85% tissue sections. Total mRNA was extracted from malignant tissues by means of QIAzol (Qiagen, USA), and RNA cleanup was carried out using the RNeasy Mini Kit (Qiagen, USA). The quantity of mRNA was measured using NanoDrop 2000 (Thermo Scientific, Germany). For complementary DNA (cDNA) synthesis, samples were synchronized at the desired concentration and subsequently reverse-transcribed by QuantiTect (Qiagen, USA) reverse transcriptase with the integral removal of genomic DNA contamination. The threshold cycle (Ct) and subsequently relative expression levels were measured using qRT-PCR by means of LightCycler 480II (Roche, Germany). Primer sets were designed for human RhoA (F: 5′AAGCAGGTAGAGTTGGCTTTGTG3′; R: 5′-ATCGGTATCTGGGTAGGAGAGG-3′), VEGF-A (F: 5′-CTACCTCCACCATGCCAAGT-3′; R: 5′-CCACTTCGTGATGATTCTGC-3′), and VEGFR2 (F: 5′-CATAGTTGTCGTTGTAGGGTA-3′; R: 5′-CATTTAGTTCAGTTCTTGCT-3′). The reaction master mix contained cDNA (~ 200 ng/μL), 10 × SYBER Green (Nanohelix, South Korea), and primers (~ 200–600 pg/μL) were used to have a reaction in a total volume of 25 μL. Thermal cycling included a first denaturation step at 95 °C for 15 min, followed by 40 cycles consisting of 24 s at 95 °C and 35 s at 62 °C. Duplicate reactions were carried out for a single sample. Quantification of fold changes was computed using 2−ΔΔct formula. The expression of the hypoxanthine-guanine phosphoribosyltransferase gene (HGPRT; F: 5′-TGGACAGGACTGAACGTCTT-3′; R: 5′-CCAGCAGGTCAGCAAAGAAT-3′) was used as an internal normalizing control.
Sample size was calculated based on the formula of mean comparisons and reach 100 subjects after considering the level of significance of α = 0.05 (two-sided), statistical power (1 − β) = 80%, and comparing the mean (SD) relative expression levels of intercellular adhesion molecule-1 (ICAM1) gene (high fiber 1.77 ± 0.47 vs. low fiber 1.94 ± 0.38) provided by Hermsdorff et al. . On the other hand, for conducting principal component factor analysis (PCA), Gorsuch  and Hair et al.  recommended that sample size can be considered at least 100. A ratio of a sample size to observe variables was also highly suggested to be considered at 10 [35, 36]. Taken together, to meet the needed sample size for conducting PCA and to generalize the findings from a sample to a wider primary population, the sample size was estimated to be 177 subjects with BC by considering the necessity of keeping potential covariates controlled in the analysis. The normal distribution of quantitative variables was assessed by using the Kolmogorov-Smirnov test and plotting histogram. Outlier data were assigned by the plotting box plot and removed. Linear regression analysis was carried out to evaluate the correlation magnitude between protein intakes from different dietary sources and fold changes of studied genes. Standardized β coefficients (β) and adjusted β (βadj.) were obtained from crude (unadjusted) and multivariate (adjusted) regression models, respectively. Models detailing the multivariate linear regression analyses consisted of following variables, just in case, when basic covariates were inter-correlated. The covariates were listed as follow: daily intake of energy (kcal/day); fat (g/day); saturated fat (g/day); cholesterol (mg/day); iron (mg/day); folate (μg/day); dietary, crude, soluble, and insoluble fiber (g/day); weight (kg); waist circumference (cm); hip circumference (cm); waist to hip ratio; body mass index (BMI; kg/m2); tumor size (cm); age at diagnosis (years); age at menarche (years); age at first childbirth (years); number of abortion; number of pregnancies; mean duration of breastfeeding (months); and mean duration of oral contraceptive usage (months). The PCA was carried out to derive appropriate dietary protein components (patterns). Accordingly, eight groups from different dietary protein sources were included in PCA. The eigenvalue is a parameter estimated to represent the sum of the variance of all the variables that can be explained by a given principal component and must be greater than 1.00, and the scree plot was applied to simplify making a decision in determining the number of components to be retained . Factor loadings, the correlation between each variable and a certain component, are presented in the component matrix . Since rotation is an important procedure to interpret the retained component, an orthogonal rotation in terms of the varimax procedure with Kaiser’s normalization was carried out . Based on the rotated component matrix, a factor score for each subject on each component is computed to calculate individual’s scores on each variable involved in a component . Indeed, an individual’s score obtained on each variable included in a component was multiplied by the factor loading for the particular variable. The sum of a person’s factor score on a component was then calculated. Subsequently, linear regression analysis was carried out to evaluate the correlation magnitude between the factor scores estimated for each component (pattern of protein sources) and fold changes of studied genes. Scatter plot was used to illustrate the correlations between the dietary patterns of protein and studied genes. Logistic regression analysis was performed to measure the odds ratio (OR) and 95% confidence interval (95% CI). Subgroup analyses based on the dichotomous status of the hormonal receptor of breast tumors (ER, PR, and HER-2) were done to control possibly the related potential effects. Further stratification analyses were relevant to the status of clinical outcomes such as tumor grade, vascular invasion (VI), and involvement of axillary lymph node metastasis (ALNM) to show whether an interested nutrigenomic correlation could be assigned in certain clinicopathological features of BC disease. Data analyses were carried out with the SPSS statistical software, version 16.0 (SPSS Inc., Chicago, IL, USA). All p values were two-tailed, and below 0.05 was considered statistically significant.
General demographic and descriptive characteristics of 177 study participants are shown in Table 1. The mean age at diagnosis of this study population was 46.67 ± 9.03 years. The relative number of premenopausal participants was 65.5% (116/177). The relative frequency of hormone receptor-positive tumors included 86.2% (131/152) of ER+, 84.8% (128/151) of PR+, and 18.1% (27/149) of HER2+ in the whole study population (p < 0.001). Histopathologic outcomes showed 81% (94/116) of patients with VI+ and 64.3% (99/154) with involvement of the ALNM in all available recorded data. Tumor grade II was significantly constituted in 73.3% (107/146) of study subjects.
The correlation of protein intakes and fold changes of RhoA, VEGF-A, and VEGFR2
Findings of univariate and multivariate linear regression analyses to show the association between participant’s protein intake and fold changes in the expression of RhoA, VEGF-A, and VEGFR2 were summarized in Table 2. Moreover, findings of linear regression analyses to show the correlation of dietary proteins with fold changes of genes of interest were presented based on hormonal receptor status (Tables 3 and 4) and clinical-pathological relevant outcomes (Tables 5 and 6).
High total amounts of protein intake (βadj. = 0.178, p = 0.044), animal proteins (βadj. = 0.192, p = 0.043), and red meat protein (βadj. = 0.348, p = 0.001) were observed to associate significantly with the overexpression of RhoA. The ratio of total protein to dietary fiber (TP:DF ratio, βadj. = 0.216, p = 0.014) and animal proteins to plant proteins (AP:PP ratio, βadj. = 0.227, p = 0.010) was positively associated with RhoA expression (Table 2). Moreover, protein intakes from nuts and cereals were inversely associated with RhoA expression in the whole sample population, when adjustments were made for confounders (βadj. = − 0.180, p = 0.034).
In case of ER+ feature, the multivariate-adjusted models showed that animal proteins (βadj. = 0.230, p = 0.045), red meat protein (βadj. = 0.297, 0.002), and the ratio of TP:DF (βadj. = 0.249, p = 0.020) could associate positively with RhoA expression (Table 3). In PR+ breast tumors, greater red meat consumed (βadj. = 0.305, p = 0.002) and high ratio of TP:DF (βadj. = 0.229, p = 0.034) were associated with higher fold changes of RhoA expression (Table 4). Protein obtained from nuts and cereals correlated inversely with fold changes of RhoA expression (ER+: βadj. = − 0.227, p = 0.045; PR+: βadj. = − 0.225, p = 0.049).
High intake levels of animal proteins (βadj. = 0.238, p = 0.041), as well as the ratios of TP:DF (βadj. = 0.242, p = 0.040) and AP:PP (βadj. = 0.273, p = 0.022), displayed significantly positive correlations with RhoA expression, whenever ALNM+ was selected (Table 5). Also, protein provided by red meat (βadj. = 0.230, p = 0.039) and dairy products (βadj. = 0.271, p = 0.030) positively correlated with the RhoA expression in ALNM+ status. High protein intake from nuts and cereals was inversely associated with RhoA expression in ALNM+ patients (βadj. = − 0.239, p = 0.043). Animal proteins (βadj. = 0.313, p = 0.007), red meat protein (βadj. = 0.358, p = 0.001), and protein from dairy products (βadj. = 0.275, p = 0.035) positively correlated with RhoA expression in subjects who had pathologic diagnosis of VI+. High protein intake from nuts and cereals (βadj. = − 0.234, p = 0.045) was inversely associated with RhoA expression in VI+ status (Table 6).
High total protein (βadj. = 0.200, p = 0.029), residual total protein (βadj. = 0.206, p = 0.029), and red meat protein (βadj. = 0.262, p = 0.001) correlated significantly with VEGF-A expression (Table 2). The overexpression of VEGF-A correlated with proteins including processed red meat (βadj. = 0.184, p = 0.043), seafood (βadj. = 0.180, p = 0.038), and dairy products (βadj. = 0.183, p = 0.030) as well.
The correlations of red meat protein (βadj. = 0.283, p = 0.002) and residual total protein intake (βadj. = 0.311, p = 0.004) with fold changes in the expression of VEGF-A were observed in ER+ BC patients (Table 3). Similarly, red meat protein (βadj. = 0.271, p = 0.004), total protein (βadj. = 0.296, p = 0.005), and its residual variable (βadj. = 0.336, p = 0.002) positively correlated with VEGF-A expression in PR+ BC patients (Table 4).
In ALNM+ patients, only red meat intake was in association with VEGF-A expression at crude (β = 0.294, p = 0.005) and adjusted (βadj. = 0.321, p = 0.003) models (Table 5). On the other hand, subpopulation with ALNM− showed a positive correlation between total protein and upregulation of VEGF-A (βadj. = 0.339, p = 0.019). In patients with VI+ (Table 6), there was a positive correlation between red meat protein and VEGF-A expression (βadj. = 0.346, p = 0.002). Similarly, the high dairy protein could associate with the overexpression of VEGF-A (βadj. = 0.263, p = 0.020) in VI+ status.
High intake of total protein (βadj. = 0.217, p = 0.013) and animal proteins (βadj. = 0.237, p = 0.007) as well as protein provided by red meat (βadj. = 0.304, p < 0.001) and processed meats (βadj. = 0.191, p = 0.039) correlated significantly with the overexpression of VEGFR2 gene (Table 2). A greater ratio of AP: PP (βadj. = 0.283, p = 0.001) was associated with the overexpressed VEGFR2.
Animal proteins (ER+: βadj. = 0.299, p = 0.003; PR+: βadj. = 0.296, p = 0.004), red meat protein (ER+: βadj. = 0.315, p = 0.001; PR+: βadj. = 0.315, p = 0.001), and the ratio of AP:PP (ER+: βadj. = 0.331, p < 0.001; PR+: βadj. = 0.331, p = 0.001) positively correlated with fold changes in the expression of VEGFR2 (Tables 3 and 4).
Total protein was positively correlated with VEGFR2 expression (βadj. = 0.295, p = 0.009) in ALNM+ BC patients (Table 5). There was also a positive correlation between protein from processed red meat and VEGFR2 expression (βadj. = 0.265, p = 0.030). The ratios of TP:DF (βadj. = 0.305, p = 0.010) and AP:PP (βadj. = 0.290, p = 0.018) were the variables significantly correlated with fold change in the expression levels of VEGFR2 at subgroup of VI+ (Table 6). In patients with VI+, red meat intake was strongly correlated with VEGFR2 expression (βadj. = 0.467, p < 0.001).
PCA-based dietary patterns of proteins
By performing PCA analyses, three important dietary protein patterns generated as (1) “whole meat,” (2) “legume dairy products,” and (3) “plant proteins” whereby they all three could explain 49.5% of the total variances accounted for (Additional file 1: Table S1). The χ2 for Bartlett’s test of sphericity was statistically significant at p < 0.001, and the Kaiser-Meyer-Olkin measure of sampling adequacy showed a score of 0.521.
The mean intake of nutrients was compared among each component (protein pattern) stratified by tertile and summarized in Additional file 2: Table S2. Individuals at the highest adherence to the “whole meat” pattern consumed higher intake levels of saturated fat (p = 0.015) and cholesterol (p = 0.030). The highest tertile of “legume dairy products” pattern consumed higher intake of energy (p = 0.008), dietary fiber (p < 0.001), insoluble fiber (p = 0.015), crude fiber (p < 0.001), iron (p = 0.001), and folate (p < 0.001) than the intake amounts observed in the lowest tertile. Women in the highest tertile of “plants” pattern had less cholesterol intake than the lowest tertile (p = 0.008).
The first pattern including protein intakes from seafood, poultries, red meats, and processed meats, in terms of “whole meat,” was inversely associated with RhoA expression in ALNM+ (βadj. = − 0.253, p = 0.033) and positively correlated with VEGFR2 expression in VI+ patients (βadj. = 0.288, p = 0.016). Logistic regression analysis showed that the second quartile of “whole meat” pattern appeared to associate inversely with fold changes in the expression of studied genes (ORRhoA = 0.24, 95% CI 0.07–0.83; ORVEGF-A = 0.26, 95% CI 0.07–0.97; ORVEGFR2 = 0.27, 95% CI 0.08–0.96). However, the trend of ORs was not significant even after the adjustments made for covariates.
The second pattern including protein intakes from milk, dairy products, and legumes was positively associated with fold changes in the expression levels of RhoA (βadj. = 0.249, p = 0.031) and VEGF-A (βadj. = 0.297, p = 0.019) in BC patients with VI+. This pattern also correlated with overexpressed VEGF-A of those patients classified as ALNM+ (βadj. = 0.330, p = 0.013). No significant association was observed between the plant protein pattern (protein intakes from fruits, vegetables, soybean, potato, cereals, nuts, and seeds) and fold change of the expressions of interested genes.
The correlation between dietary patterns of protein and studied genes in menopausal status are shown in Fig. 1. In premenopausal status, the “whole meat” as the first pattern was associated inversely with RhoA (βadj. = − 0.285, p = 0.014) and positively with VEGFR2 (βadj. = 0.300, p = 0.009). Protein provided by “legume dairy products” as the second pattern significantly correlated with the overexpression of VEGF-A in premenopausal women (βadj. = 0.356, p = 0.029).
Findings of the present molecular epidemiologic study provided supports in associations between high consumption of protein and upregulation of RhoA and VEGF-VEGFR2. These candidate genes are functionally significant in lymphangiogenesis as a poor determinant of prognosis [23, 24, 27]. Two sets of results were mainly developed. First, proteins from red meat and dairy products were demonstrated to have a significant correlation with the overexpression of RhoA in favor of growing tumor cells to lymph nodes (ALNM) and VI+ patients. Similarly, in both clinical subgroups, red meat correlated with upregulated VEGF-A and VEGFR2 to promote ALNM and VI. Secondly, the upstream regulatory effects of ER and PR signaling seem to be a crucial modifier in specifying what type of dietary protein can modify the transcription levels of study gene.
To our knowledge, this is the first study to investigate the association between dietary protein intake and the expression of RhoA, VEGF-A, and VEGFR2. A considerable number of prospective cohort-based studies showed that dietary protein can associate with substantially increased BC risk [7, 8, 37]. Cho and co-workers , in a large prospective cohort study, reported that red meat intake strongly elevates BC risk. Lately, Wu et al.  conducted a dose-response meta-analysis of prospective studies and revealed that protein obtained by red and processed meat may enhance BC risk. In another cohort study, total red meat intake was positively associated with increased risk of metastasis of breast tumor . Furthermore, in Sweden, Larsson et al.  suggested that fried meat intake may enhance the risk of ER+/PR− breast carcinogenesis. By contrast, no population-based study exists showing how dietary protein can affect the molecular responses in metastasis of BC tumors .
The overexpression of RhoA is the most common feature of metastasis which is unraveled in association with modifiable dietary risk factors [14, 16]. Present findings showed that animal proteins, especially red meat, can increase strongly the expression levels of RhoA. This result was also re-emphasized particularly in ER+ and PR+ BCs, showing the possible interfering impacts of active ER and PR signaling over this nutrigenomic model. In detail, plant sources of protein specifically from nuts and cereals correlated inversely with RhoA expression varied dramatically by ER+ and PR+ statuses. Several different mechanisms have been proposed to figure out why red meat consumption causes susceptibility to cancer development. Carcinogenic heterocyclic amines (HCAs) formed in meat are contingent with the type of meat (red, white) and the factors associated with the cooking process such as temperature and duration [39, 40]. In vitro experiments showed that specific HCAs, involved 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP), are a byproduct component with estrogenic effects that could induce a mitogenic response through ER signaling [39, 40]. Moreover, PhIP could also induce upregulation of PR in MCF-7 cells . High intake of animal proteins could increase the acid load in the circulation , thereby promoting the formation of stress fiber mediated by RhoA in enhancing focal adhesion mechanism . Oster and co-workers  undertook microarray analyses on skeletal muscle tissue of 253 offspring of German gilts fed with iso-caloric diets showed that the transcriptional pathway of Rho GTPase decreased in “low-protein/high-carbohydrate” diet after 188 days follow-up at postnatal subjects. Garcia et al.  indicated that treatment with arachidonic acid in MDA-MB-435 human melanoma cells can activate RhoA promoting cell adhesion via p38 MAPK-RhoA signaling pathway. Arachidonic acid is rich in animal protein, especially red meat. Thromboxane A2 is a pro-inflammatory factor produced from arachidonic acid and potently participates in platelet aggregation and vascular contraction . Thromboxane A2 is a potent enhancer of RhoA transcription . Red meat is rich in N-glycolylneuraminic acid (Neu5Gc). Samraj et al.  indicated induced hepatocellular tumorigenesis by Neu5Gc acid intervention in human-like Neu5Gc-deficient mice giving support to epidemiologic findings showing the correlation of red meat consumption and risk of cancer incidence. Over an interventional clinical trial-manipulated high intake of red meat in the diet of patients with irritable bowel disease, the expression of RhoA increased in colon tissue . Although different aspects exist to support the carcinogenic effects of high protein consumption particularly from animal red meat, our findings can provide new insight indicating significant alterations of RhoA expression levels in association with independent dietary factors such as red meat and dairy products.
Animal proteins especially red meat and dairy products were correlated with RhoA overexpression dependent on the involvement of ALNM and VI positivity, supported our hypothesis expressed the triple axis of diet, RhoA transcription and lymphatic anomalies in BC participants. Although, based on PCA outputs, “whole meat” was associated inversely with RhoA expression in ALNM+ and premenopausal status. It is noteworthy to mention that at the present study, white meat (poultries and seafood) was consumed 70% more than red and processed meat. This would explain why the extent of correlation by every individual food item could be hindered or even masked by other prominent dietary counterparts in the model defining the component. Given that very limited information exist to discuss how dietary factors can contribute in the infiltration of tumor cells to the lymph nodes and vessels, epidemiologic studies are widely warranted to study the correlation of protein sources and other metastatic transcripts in future researches.
At the present study, protein from the group of foods including nuts and cereals were inversely associated with RhoA expression levels, particularly in ER+ and PR+ subjects and involvement of ALNM and VI positivity, suggesting the possible preventive correlation of nuts and cereals on RhoA overexpression. Similarly, in a large population-based case control study, Liu et al.  reported that nuts and vegetable proteins in adolescence may associate with reduced risk of BC later in life. However, a meta-analysis of prospective cohort studies has found no linear correlation between nut intake and BC risk . Our results are not consistent with those of a prospective cohort study conducted by Farvid et al.  showing that replacing legumes and nuts instead of red meat in early adulthood could reduce BC risk later in life. The nutrigenomic aspect of consuming proteins could be better represented when the planning of models is accounted for controlling the hormonal receptor status of tumors.
It is well established that VEGF-A takes part as a potent angiogenic growth factor in the nurturing malignant solid tumors mainly through binding to VEGFR2 . Molecular evidence indicated that co-expression of VEGF-A and VEGFR2 can associate with poor prognosis and worse clinical outcomes in BC patients . Present results revealed that protein intake from red meat can increase transcript levels of both VEGF-A and VEGFR2 whereby ER and PR were expressed positively. The Western diet including in particular high red meat increases the serum concentrations of free estradiol thereby promoting BC risk . However, little is known about the effects of protein quantity and quality in regulating the molecular pathways which control carcinogenesis. The IGF-1/mTOR signaling pathway is principally regulated by protein content [10, 26]. Studies suggested that dietary protein restriction would be more effective rather than putting a restriction on calories or fat intake in order to decrease IGF-1 levels by inhibiting the PI3K/AKT/mTOR pathway [9, 10, 47]. In an experimental study, mice fed a low-protein diet showed 45% smaller tumor size and 30% less serum concentration of IGF-1 than high protein consumers . In a xenograft model of prostate cancer, mice fed low-protein diet (7% of total calorie) represented a significant decrease in expression of enhancer of zeste homolog 2  which could enhance VEGF-A expression . Isocaloric diet in animal consumed plant proteins significantly inhibit BC growth in human xenograft models of tumorigenesis seems to be mediated by the reduction in serum IGF-1 levels and downregulation of intratumor mTOR activity . Porcine hepatocytes and HepG2 cell line exposed to 4× amino acid concentration showed increased expressions of IGF-1, peroxisome proliferator-activated receptor γ (PPARγ), and activated protein-2 (ap-2) are also significantly overexpressed . Downregulation of ap-2 can lead to the inhibition of VEGF expression in human H1299 cell line . High transcript levels of PPARγ could stimulate angiogenesis in various carcinoma through increasing VEGF expression . Our findings showed that protein obtained from red meat correlated with upregulated VEGF-A and processed meat to higher levels of VEGFR2 in association with the feature of ALNM+. The overexpression of VEGF-A and VEGFR2 were attributed to high protein intake from red meat when patients’ tumor accounted for the involvement of VI. Moreover, based on PCA data, the protein obtained from “whole meat” was positively associated with VEGFR2 expression in VI+ patients and premenopausal status. No previous published data is available to compromise how protein does associate with lymphangiogenesis dependent on alteration in genomic profile (such as VEGFR2). Given the association of protein intake with overexpression of VEGF and VEGFR2 suggest that this nutrigenomic model can correlate with determining the predispose population to spread tumor cells into lymph node and vessels as important clinicopathologic variables.
The results of this study showed that proteins obtained from legumes (beans, peas, and lentil) and dairy products correlated positively with the increased fold change in the expression of VEGF-A at premenopausal status or who characterized by ALNM+ and VI+. While nutraceutical effects of legume consumption raise the issue of possible anti-carcinogenic effects by active ingredients assigned to legumes , some epidemiologic studies were unsuccessful to reveal associations between legume intake and cancers of the breast , prostate , and colon . Our findings showed the positive correlation between “legume dairy products” pattern and VEGF-A expressions while it is noteworthy to highlight that the association of legumes might be hampered by the significant correlation of dairy products to the variances displayed by component 2 (legume dairy products). The PCA-independent data was also reassured that, despite legumes, just the group of dairy products was correlated with the upregulation of VEGF-A. The nature of milk proteins may explain the positive association between dairy products and the overexpression of VEGF-A. Cow’s milk contains two important glycoproteins, lactadherin, and angiogenin-2, which they can modulate angiogenesis process . Lactadherin which is secreted into milk promotes VEGF-dependent Akt phosphorylation consequently induces neovascularization . Moreover, hormone-containing cow’s milk may predispose milk consumers to increased IGF-1 and estrogen in time-dependent condition linked as a potent risk factor for BC . The incidence of mammary tumor in rats which exposed to methylnitrosourea and fed a casein-based diet was 80% higher than rats fed with a soy protein diet (42% incidence rate) . By contrast, a recent meta-analysis indicates that increased consumption of total dairy food, but not milk, may associate with a reduced risk of BC . Higher intake of legumes and dairy products may associate with high estradiol levels and increased risk of BC in postmenopausal women . Our findings can provide new insight specifying significant alterations of VEGF-A and VEGFR2 expression levels in association with legumes and dairy products.
This study had some limitations. Aside from the possibility of recall bias which cannot be completely excluded, the sample size was small. A larger population study would be desirable in order to perform the molecular epidemiologic study.
In conclusion, the findings suggested that high intake of animal proteins especially red meat may associate with the overexpression of RhoA and VEGF-VEGFR2 in patients characterized by the involvement of ALNM and VI. Wherein the combination of legume dairy products correlated with RhoA and VEGF-A, either a clinicopathologic feature of lymphatic or vascular metastasis was remarkable in BC patients. Less intake of “whole meat” was associated with less fold change in the expression of interested genes and may suggest the prevention of metastasis in BC patients. Thus, for future studies, it is highly recommended to study the association between different dietary sources of protein and a larger genomic profile including various metastatic and angiogenic genes.
Availability of data and materials
The data that support the findings of this study are available from TBZMED, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. However, data are available from the authors upon reasonable request and with permission from TBZMED.
Torre LA, Bray F, Siegel RL, Ferlay J. Global cancer statistics, 2012. Cancer J Clin. 2015;65:87–108.
Mousavi SM, Montazeri A, Mohagheghi MA, Jarrahi AM, Harirchi I, Najafi M, Ebrahimi M. Breast cancer in Iran: an epidemiological review. Breast J. 2007;13(4):383–91.
Redig AJ, McAllister SS. Breast cancer as a systemic disease: a view of metastasis. J Intern Med. 2013;274(2):113–26.
Ponder BA. Cancer genetics. Nature. 2001;411(6835):336–41.
Meadows GG. Diet, nutrients, phytochemicals, and cancer metastasis suppressor genes. Cancer Metastasis Rev. 2012;31(3-4):441–54.
Wu J, Zeng R, Huang J, Li X, Zhang J, Ho JC-M, Zheng Y. Dietary protein sources and incidence of breast cancer: a dose-response meta-analysis of prospective studies. Nutrients. 2016;8(11):730.
Cho E, Chen WY, Hunter DJ, Stampfer MJ, Colditz GA, Hankinson SE, Willett WC. Red meat intake and risk of breast cancer among premenopausal women. Arch Intern Med. 2006;166(20):2253–9.
Farvid MS, Cho E, Chen WY, Eliassen AH, Willett WC. Dietary protein sources in early adulthood and breast cancer incidence: prospective cohort study. BMJ. 2014;348:g3437.
Fontana L, Adelaiye RM, Rastelli AL, Miles KM, Ciamporcero E, Longo VD, Nguyen H, Vessella R, Pili R. Dietary protein restriction inhibits tumor growth in human xenograft models of prostate and breast cancer. Oncotarget. 2013;4(12):2451–61.
Levine ME, Suarez JA, Brandhorst S, Balasubramanian P, Cheng C-W, Madia F, Fontana L, Mirisola MG, Guevara-Aguirre J, Wan J. Low protein intake is associated with a major reduction in IGF-1, cancer, and overall mortality in the 65 and younger but not older population. Cell Metabol. 2014;19(3):407–17.
Z’graggen K, Warshaw AL, Werner J, Graeme-Cook F, Jimenez RE, Fernández-del Castillo C. Promoting effect of a high-fat/high-protein diet in DMBA-induced ductal pancreatic cancer in rats. Ann Surg. 2001;233(5):688.
Mizukami Y, Kohgo Y, Chung DC. Hypoxia inducible factor-1–independent pathways in tumor angiogenesis. Clin Cancer Res. 2007;13(19):5670–4.
Rosenthal D, Brenner J, Merajver S. Rho proteins in cancer. In: The Rho GTPases in cancer. Van Golen KL: Springer; 2010. p. 29–42.
Gilkes DM, Xiang L, Lee SJ, Chaturvedi P, Hubbi ME, Wirtz D, Semenza GL. Hypoxia-inducible factors mediate coordinated RhoA-ROCK1 expression and signaling in breast cancer cells. Proc Natl Acad Sci. 2014;111(3):E384–93.
Hall A. Rho GTPases and the actin cytoskeleton. Science. 1998;279(5350):509–14.
Alho I, Costa L, Bicho M, Coelho C. Low molecular weight protein tyrosine phosphatase isoforms regulate breast cancer cells migration through a RhoA dependent mechanism. PloS one. 2013;8(9):e76307.
Hebels DG, Sveje KM, de Kok MC, van Herwijnen MH, Kuhnle GG, Engels LG, Vleugels-Simon CB, Mares WG, Pierik M, Masclee AA. Red meat intake-induced increases in fecal water genotoxicity correlate with pro-carcinogenic gene expression changes in the human colon. Food Chem Toxicol. 2012;50(2):95–103.
Yang X, Huang H-C, Yin H, Alpern RJ, Preisig PA. RhoA required for acid-induced stress fiber formation and trafficking and activation of NHE3. Am J Physiol Renal Physiol. 2007;293(4):F1054–64.
Miller KD, Dul CL. Breast cancer: the role of angiogenesis and antiangiogenic therapy. Hematol Oncol Clin North Am. 2004;18(5):1071–86.
Fox SB, Generali DG, Harris AL. Breast tumour angiogenesis. Breast Cancer Res. 2007;9(6):216.
Ferrara N, Gerber H-P, LeCouter J. The biology of VEGF and its receptors. Nat Med. 2003;9(6):669–78.
Rydén L, Linderholm B, Nielsen NH, Emdin S, Jönsson P-E, Landberg G. Tumor specific VEGF-A and VEGFR2/KDR protein are co-expressed in breast cancer. Breast Cancer Res Treat. 2003;82(3):147–54.
Nagy JA, Vasile E, Feng D, Sundberg C, Brown LF, Detmar MJ, Lawitts JA, Benjamin L, Tan X, Manseau EJ. Vascular permeability factor/vascular endothelial growth factor induces lymphangiogenesis as well as angiogenesis. J Exp Med. 2002;196(11):1497–506.
Rapisarda A, Melillo G. Role of the VEGF/VEGFR axis in cancer biology and therapy. In: Advances in cancer research. Volume 114. Daar IO: Elsevier; 2012. p. 237–67.
Bos R, Zhong H, Hanrahan CF, Mommers EC, Semenza GL, Pinedo HM, Abeloff MD, Simons JW, van Diest PJ, van der Wall E. Levels of hypoxia-inducible factor-1α during breast carcinogenesis. J Natl Cancer Inst. 2001;93(4):309–14.
Klement R, Fink M. Dietary and pharmacological modification of the insulin/IGF-1 system: exploiting the full repertoire against cancer. Oncogenesis. 2016;5(2):e193.
Mohammed RAA, Green A, El-Shikh S, Paish EC, Ellis IO, Martin SG. Prognostic significance of vascular endothelial cell growth factors -A, -C and -D in breast cancer and their relationship with angio- and lymphangiogenesis. Br J Cancer. 2007;96:1092–100.
Pirouzpanah S, Taleban F, Mehdipour P, Atri M, Hooshyareh-Rad A, Sabour S. The biomarker-based validity of a food frequency questionnaire to assess the intake status of folate, pyridoxine and cobalamin among Iranian primary breast cancer patients. Eur J Clin Nutr. 2014;68(3):316–23.
Pirouzpanah S, Taleban F-A, Sabour S, Mehdipour P, Atri M, Farrin N, Houshyar-Rad A, Jalali-Farahani S, Karimian Khosroshahi N. Validation of food frequency questionnaire to assess folate intake status in breast cancer patients. Razi J Med Sci. 2012;18(92):31–41.
Pirouzpanah S, Taleban F-A, Mehdipour P, Atri M, Foroutan-Ghaznavi M. Plasma total homocysteine level in association with folate, pyridoxine, and cobalamin status among iranian primary breast cancer patients. Nutr Cancer. 2014;66(7):1097–108.
Ghaffarpour M, Houshiar-Rad A, Kianfar H. The manual for household measures, cooking yields factors and edible portion of foods, vol. 7. Tehran: Nashre Olume Keshavarzy; 1999.
Willet W. Nutritional epidemiology. New York: Oxford University Press; 1998.
Hermsdorff HHM, Zulet MÁ, Puchau B, Martínez JA. Fruit and vegetable consumption and proinflammatory gene expression from peripheral blood mononuclear cells in young adults: a translational study. Nutr Metab. 2010;7(1):42–53.
Gorsuch RL. Factor analysis Lawrence Erlbaum Associates. Hillsdale; 1983.
Hair JF, Black WC, Babin BJ, Anderson RE. Exploratory factor analysis. In: Hair JF, Black WC, Babin BJ, Anderson RE, editors. Multivariate data analysis: Global edition. 7 edn. Upper Saddle River: Pearson Higher Education; 2010. p. 89–149.
MacCallum RC, Widaman KF, Zhang S, Hong S. Sample size in factor analysis. Psychol Methods. 1999;4(1):84.
Larsson SC, Bergkvist L, Wolk A. Long-term meat intake and risk of breast cancer by oestrogen and progesterone receptor status in a cohort of Swedish women. Eur J Clin Nutr. 2009;45(17):3042–6.
Inoue-Choi M, Sinha R, Gierach GL, Ward MH. Red and processed meat, nitrite, and heme iron intakes and postmenopausal breast cancer risk in the NIH-AARP Diet and Health Study. Int J Cancer. 2016;138(7):1609–18.
Gooderham NJ, Zhu H, Lauber S, Boyce A, Creton S. Molecular and genetic toxicology of 2-amino-1-methyl-6-phenylimidazo [4, 5-b] pyridine (PhIP). Mutat Res. 2002;506:91–9.
Lauber SN, Ali S, Gooderham NJ. The cooked food derived carcinogen 2-amino-1-methyl-6-phenylimidazo [4, 5-b] pyridine is a potent oestrogen: a mechanistic basis for its tissue-specific carcinogenicity. Carcinogenesis. 2004;25(12):2509–17.
Oster M, Murani E, Metges C, Ponsuksili S, Wimmers K. High-and low-protein gestation diets do not provoke common transcriptional responses representing universal target-pathways in muscle and liver of porcine progeny. Acta Physiol. 2014;210(1):202–14.
Garcia MC, Ray DM, Lackford B, Rubino M, Olden K, Roberts JD. Arachidonic acid stimulates cell adhesion through a novel p38 MAPK-RhoA signaling pathway that involves heat shock protein 27. J Biol Chem. 2009:20936–45.
Olivia MY, Brown JH. G Protein–coupled receptor and RhoA-stimulated transcriptional responses: links to inflammation, differentiation, and cell proliferation. Mol Pharmacol. 2015;88(1):171–80.
Samraj AN, Pearce OM, Läubli H, Crittenden AN, Bergfeld AK, Banda K, Gregg CJ, Bingman AE, Secrest P, Diaz SL. A red meat-derived glycan promotes inflammation and cancer progression. Proc Natl Acad Sci. 2015;112(2):542–7.
Liu Y, Colditz GA, Cotterchio M, Boucher BA, Kreiger N. Adolescent dietary fiber, vegetable fat, vegetable protein, and nut intakes and breast cancer risk. Breast Cancer Res Treat. 2014;145(2):461–70.
Harris HR, Bergkvist L, Wolk A. An estrogen-associated dietary pattern and breast cancer risk in the Swedish Mammography Cohort. Int J Cancer. 2015;137(9):2149–54.
Clemmons DR, Seek MM, Underwood LE. Supplemental essential amino acids augment the somatomedin-C/insulin-like growth factor I response to refeeding after fasting. Metabolism. 1985;34(4):391–5.
Geng J, Li X, Zhou Z, Wu C-L, Bai X, Dai M. EZH2 promotes tumor progression via regulating VEGF-A/AKT signaling in non-small cell lung cancer. Cancer Lett. 2015;359(2):275–87.
Wan X, Wang S, Xu J, Zhuang L, Xing K, Zhang M, Zhu X, Wang L, Gao P, Xi Q. Dietary protein-induced hepatic IGF-1 secretion mediated by PPARγ activation. PloS One. 2017;12(3):e0173174.
Fu L, Shi K, Wang J, Chen W, Shi D, Tian Y, Guo W, Yu W, Xiao X, Kang T. TFAP2B overexpression contributes to tumor growth and a poor prognosis of human lung adenocarcinoma through modulation of ERK and VEGF/PEDF signaling. Mol Cancer. 2014;13:89.
Bishop-Bailey D. PPARs and angiogenesis. Biochem Soc Trans. 2011;39:1601–5.
Sánchez-Chino X, Jiménez-Martínez C, Dávila-Ortiz G, Álvarez-González I, Madrigal-Bujaidar E. Nutrient and nonnutrient components of legumes, and its chemopreventive activity: a review. Nutr Cancer. 2015;67(3):401–10.
Aune D, De Stefani E, Ronco A, Boffetta P, Deneo-Pellegrini H, Acosta G, Mendilaharsu M. Legume intake and the risk of cancer: a multisite case–control study in Uruguay. Cancer Causes Control. 2009;20(9):1605–15.
Kirsh VA, Peters U, Mayne ST, Subar AF, Chatterjee N, Johnson CC, Hayes RB. Prospective study of fruit and vegetable intake and risk of prostate cancer. J Natl Cancer Inst. 2007;99(15):1200–9.
Park Y, Subar AF, Kipnis V, Thompson FE, Mouw T, Hollenbeck A, Leitzmann MF, Schatzkin A. Fruit and vegetable intakes and risk of colorectal cancer in the NIH–AARP diet and health study. Am J Epidemiol. 2007;166(2):170–80.
Silvestre J-S, Théry C, Hamard G, Boddaert J, Aguilar B, Delcayre A, Houbron C, Tamarat R, Blanc-Brude O, Heeneman S. Lactadherin promotes VEGF-dependent neovascularization. Nat Med. 2005;11(5):499–506.
Parodi PW. Dairy product consumption and the risk of breast cancer. J Am Coll Nutr. 2005;24(sup6):556S–68S.
Hawrylewicz E, Huang H, Blair W. Dietary soybean isolate and methionine supplementation affect mammary tumor progression in rats. J Nutr. 1991;121(10):1693–8.
Dong J, Zhang L, He K, Qin L-Q. Dairy consumption and risk of breast cancer: a meta-analysis of prospective cohort studies. Breast Cancer Res Treat. 2011;127(1):23–31.
Fung TT, Schulze MB, Hu FB, Hankinson SE, Holmes MD. A dietary pattern derived to correlate with estrogens and risk of postmenopausal breast cancer. Breast Cancer Res Treat. 2012;132(3):1157–62.
We are grateful for all patients who participated and valued colleagues of Nour-Nejat Hospital and Nemoone Laboratory.
We are thankful to the Drug and Applied Research Center (grant no.: 94-101) and TBZMED for providing financial support (grant no.: T/A/67) to conduct this study planned as an M.Sc. thesis numbered T/A/67. We are also thankful to the Students’ Research Committee and TBZMED for the support of this research.
Ethics approval and consent to participate
The major ethical issues were expressed to every single participant prior to the enrolment. Each patient had signed an informed consent form before the interview. The protocol of this study and questionnaires were reviewed and approved by the Ethical Committee of the Tabriz University of Medical Sciences (TBZMED; Ethic no: IR.TBZMED.REC.1394.401).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Shokri, A., Pirouzpanah, S., Foroutan-Ghaznavi, M. et al. Dietary protein sources and tumoral overexpression of RhoA, VEGF-A and VEGFR2 genes among breast cancer patients. Genes Nutr 14, 22 (2019). https://doi.org/10.1186/s12263-019-0645-7
- Breast cancer