Dietary protein sources and tumoral overexpression of RhoA, VEGF-A and VEGFR2 genes among breast cancer patients

Background 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. Methods 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. Results 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+. Conclusions 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. Electronic supplementary material The online version of this article (10.1186/s12263-019-0645-7) contains supplementary material, which is available to authorized users.


Introduction
Breast cancer (BC) is the most frequently diagnosed malignancy in females in many countries [1]. In the recent decade, the prevalence rate of BC has been increasing rapidly among Iranian women [2]. Even lately, BC is the main leading cause of female cancer-related mortality [2]. Metastasis is a critical event in malignancy addressed as an indicator of poor prognosis and causes the vast majority of cancer-related deaths [3].
Sufficient epidemiologic evidence has revealed the association between some lifestyle-related risk factors and the BC risk [4]. However, it is far less understood how dietary factors can take part in cancer progression toward the formation of metastasis [5]. 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 [6]. Moreover, in a large prospective cohort study, Cho et al. [7] 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 followup, it has been reported that high red meat intake in early adulthood might increase the risk of BC in later life [8]. 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 [11].
A review has highlighted that signaling pathways leading to cytoskeletal reprogramming are vital for cancerous cellular motility [12]. 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 [13]. 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 [15]. The configuration of active GTP bound allows Rho protein to regulate signal transduction [15]. The phosphatidylinositol-3 kinase (PI3K) activates the RhoA-ROCK signaling pathway [12]. 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 [14]. It is well-established that RhoA overexpressed in breast tumors [16]. However, a few clinical trials indicated that dietary protein intake can affect RhoA expression. Hebels et al. [17] 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 [18]. 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 proand anti-angiogenic factors in healthy tissue, whereas pathologic angiogenesis potentiates tumorigenesis linked to the distorted balance of angiogenic variables grow endothelial cells [21]. 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 [20]. 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 [20]. VEGF-A could also induce lymphangiogenesis [23]. 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 [22]. Overregulated VEGF-A is significantly associated with poor prognosis in BC [22]. 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 [25]. 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 [12]. 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 [5], 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.

Study population
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/m 2 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).

Dietary assessment
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 [31]. 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 [32].
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 reversetranscribed 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′-CCAC TTCGTGATGATTCTGC-3′), and VEGFR2 (F: 5′-CATAGTTGTCGTTGTAGGGTA-3′; R: 5′-CATT TAGTTCAGTTCTTGCT-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.

Statistical analysis
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. [33]. On the other hand, for conducting principal component factor analysis (PCA), Gorsuch [34] and Hair et al. [35] 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/m 2 ); 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 [35]. Factor loadings, the correlation between each variable and a certain component, are presented in the component matrix [35]. 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 [35]. 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 [35]. 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.

Results
General demographic and descriptive characteristics of 177 study participants are shown in 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 clinicalpathological relevant outcomes (Tables 5 and 6).

VEGFR2
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 , and crude fiber (g/day). b Fat intake (g/day), crude fiber (g/day), the age of first childbirth (years), tumor size (cm), age at diagnosis (years), and the number of pregnancies. c Energy intake (kcal/day), dietary fiber (g/day), and waist to hip ratio. d Energy intake (kcal/day), the mean duration of breastfeeding (months), tumor size (cm), and BMI (kg/m 2 ). e Energy intake (kcal/day) and tumor size (cm). f Log transformed and adjusted for fat intake (g/day), the age of menarche (years), and tumor size (cm). g Age of first childbirth (years), the age of menarche (years), and waist circumference (cm). h Log transformed and adjusted for energy intake (kcal/ day), BMI (kg/m 2 ), tumor size (cm), and age of menarche (years). i BMI (kg/m 2 ), tumor size (cm), and carbohydrate intake (g/day). j Log transformed and adjusted for fat intake (g/day) and mean duration of breastfeeding (months). k Mean duration of breastfeeding (months), physical activity rate and OCP usage. l Carbohydrate intake (g/day), the age of menarche (years), and mean duration of breastfeeding (months). m Fat intake (g/day), the age of first childbirth (years), and the number of lactation (n). n Fat intake (g/day), tumor size (cm), and carbohydrate intake (g/day). o Fat intake (g/day) and tumor size (cm). p Dietary fiber intake (g/day) and waist to hip ratio. q Fat intake (g/day). r Fat intake (g/day), red meat intake (g/day), the age of first childbirth (years), waist to hip ratio, and age of menarche (years). s Fat intake (g/day), mean duration of breastfeeding (months), and BMI (kg/m 2 ). t Log transformed and adjusted for energy intake (kcal/day) and mean duration of breastfeeding (months). u Fat intake (g/day), red meat intake (g/day), number of pregnancies (n), and physical activity rate. v Log transformed and adjusted for energy intake (kcal/day), insoluble fiber intake (g/day), and age at diagnosis (years). w Energy intake (kcal/day), insoluble fiber (g/day), and waist to hip ratio. x Fat intake (g/day), BMI (kg/m 2 ), and tumor size (cm). y Fat intake (g/day) and tumor size (cm). z Fat intake (g/day), the age of first childbirth (years), and waist to hip ratio. aa Dietary fiber intake (g/day), fat intake (g/day), and tumor size (cm). bb Mean duration of breastfeeding (months) and tumor size (cm). cc Fat intake (g/day) and age of menarche (years). dd Fat intake (g/day), BMI (kg/m 2 ), and tumor size (cm). ee Red meat intake (g/day), fat intake (g/day), and tumor size (cm). ff Fat intake (g/day), soluble fiber (g/day), and waist to hip ratio. gg Energy intake (kcal), mean duration of breastfeeding (months), and tumor size (cm). hh Red meat intake (g/ day), mean duration of breastfeeding (months), physical activity rate, and tumor size (cm). ii Log transformed and adjusted for fat intake (g/day), waist to hip ratio, and tumor size (cm). jj Log transformed and adjusted for waist circumference (cm). kk Fat intake (g/day), tumor size (cm), and waist circumference (cm). ll Carbohydrate intake (g/day), the age of menarche (years), and waist to hip ratio. mm Fat intake (g/day) and BMI (kg/m 2 )    processed meats (β adj. = 0.191, p = 0.039) correlated significantly with the overexpression of VEGFR2 gene ( Table 2) (Tables 3 and 4).

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 (OR RhoA = 0.24, 95% CI 0.07-0.83; OR VEGF-A = 0.26, 95% CI 0.07-0.97; OR VEGFR2 = 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).

Discussion
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 [7], in a large prospective cohort study, reported that red meat intake strongly elevates BC risk. Lately, Wu et al. [6] 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 [38].   Furthermore, in Sweden, Larsson et al. [37] 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 [17]. 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 [40]. High intake of animal proteins could increase the acid load in the circulation [18], thereby promoting the formation of stress fiber mediated by RhoA in enhancing focal adhesion mechanism [18]. Oster and co-workers [41] 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. [42] 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 proinflammatory factor produced from arachidonic acid and potently participates in platelet aggregation and vascular contraction [42]. Thromboxane A2 is a potent enhancer of RhoA transcription [43]. Red meat is rich in N-glycolylneuraminic acid (Neu5Gc). Samraj et al. [44] 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 trialmanipulated high intake of red meat in the diet of patients with irritable bowel disease, the expression of RhoA increased in colon tissue [17]. 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. [45] 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 [6]. Our results are not consistent with those of a prospective cohort study conducted by Farvid et al. [8] 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 [21]. Molecular evidence indicated that co-expression of VEGF-A and VEGFR2 can associate with poor prognosis and worse clinical outcomes in BC patients [27]. Present results revealed that protein intake from red meat can increase transcript levels of both VEGF-A and VEGFR2 Table 5 Correlations between intake levels of protein from different sources and fold change expression of RhoA, VEGF-A, and VEGFR2 within each category of ALNM status Dietary variables     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 [46]. 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 [10]. 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 [9] which could enhance VEGF-A expression [48]. 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 [9]. 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 [49]. Downregulation of ap-2 can lead to the inhibition of VEGF expression in human H1299 cell line [50]. High transcript levels of PPARγ could stimulate  177). a-c Adjusted for fat (g/day), cholesterol (mg/day), and mean duration of breastfeeding (months). d-f Adjusted for saturated fat (g/day), tumor size (cm), and waist to hip ratio. g-i Adjusted for crude fiber (g/day), folate (μg/day), mean duration of breastfeeding (months), and number of abortion angiogenesis in various carcinoma through increasing VEGF expression [51]. 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 [52], some epidemiologic studies were unsuccessful to reveal associations between legume intake and cancers of the breast [53], prostate [54], and colon [55]. 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 PCAindependent 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 [56]. Lactadherin which is secreted into milk promotes VEGFdependent Akt phosphorylation consequently induces neovascularization [56]. Moreover, hormonecontaining cow's milk may predispose milk consumers to increased IGF-1 and estrogen in timedependent condition linked as a potent risk factor for BC [57]. 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) [58]. 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 [59]. Higher intake of legumes and dairy products may associate with high estradiol levels and increased risk of BC in postmenopausal women [60]. 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.