Next Article in Journal
The Effect of Hypo-Hydration on Mood and Cognition Is Influenced by Electrolyte in a Drink and Its Colour: A Randomised Trial
Next Article in Special Issue
Precision Nutrition and Cancer Relapse Prevention: A Systematic Literature Review
Previous Article in Journal
Association of Low Serum l-Carnitine Levels with Peripheral Arterial Stiffness in Patients Who Undergo Kidney Transplantation
Previous Article in Special Issue
Rectal Cancer: 20% Risk Reduction Thanks to Dietary Fibre Intake. Systematic Review and Meta-Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tolerability and Safety of a Nutritional Supplement with Potential as Adjuvant in Colorectal Cancer Therapy: A Randomized Trial in Healthy Volunteers

by
Marta Gómez de Cedrón
1,*,
José Moises Laparra
2,
Viviana Loria-Kohen
3,
Susana Molina
1,
Juan Moreno-Rubio
1,4,
Juan Jose Montoya
5,
Carlos Torres
6,
Enrique Casado
4,
Guillermo Reglero
6,7 and
Ana Ramírez de Molina
1,*
1
Molecular Oncology and Nutritional Genomics of Cancer, IMDEA-Food Institute, CEI UAM+CSIC, Madrid 28049, Spain
2
Molecular Immuno-Nutrition Group, IMDEA-Food Institute, CEI UAM+CSIC, Madrid 28049, Spain
3
Nutrition and Clinical Trials Unit, GENYAL Platform IMDEA-Food Institute, CEI UAM + CSIC, Madrid 28049, Spain
4
Precision Oncology Laboratory (POL), Infanta Sofía University Hospital, Madrid 28702, Spain
5
Department of Radiology, Rehabilitation & Physiotherapy, Faculty of Medicine, Complutense University, 28040 Madrid, Spain
6
Production and Characterization of Novel Foods, CIAL, CEI UAM+CSIC, Madrid 28049, Spain
7
Production and Development of Foods for Health, IMDEA-Food Institute, CEI UAM+CSIC, Madrid 28049, Spain
*
Authors to whom correspondence should be addressed.
Nutrients 2019, 11(9), 2001; https://doi.org/10.3390/nu11092001
Submission received: 23 June 2019 / Revised: 2 August 2019 / Accepted: 21 August 2019 / Published: 24 August 2019

Abstract

:
Bioactive supplements display relevant therapeutic properties when properly applied according to validated molecular effects. Our previous research efforts established the basis to develop a dietary supplement based on a Rosmarinus officinalis supercritical extract. This was enriched in phenolic diterpenes (RE) with proven properties against signaling pathways involved in colon tumorigenesis, and shark liver oil rich in alkylglycerols (AKG) as a bioactive lipid vehicle to improve RE bioavailability and synergize with the potential therapeutic action of the extract. Herein, we have investigated the tolerability and safety of the supplement and the biological and molecular effects from an immuno-nutritional perspective. Sixty healthy volunteers participated in a six week, double-blind, randomized parallel pilot study with two study arms: RE-AKG capsules (CR) and control capsules (CC). Mean age (±SD) of volunteers was 28.32 (±11.39) and 27.5 (±9.04) for the control and the study groups, respectively. Safety of the CR product consumption was confirmed by analyzing liver profile, vital constants, and oxidation markers (LDLox in blood and isoprostanes and thromboxanes in urine). The following were monitored: (1) the phenotyping of plasmatic leukocytes and the ex vivo response of lipopolysaccharide (LPS)-stimulated peripheral blood mononuclear cells (PBMCs); (2) expression of genes associated with immune-modulation, inflammation, oxidative stress, lipid metabolism, and tumorigenesis; and (3) the correlation of selected genetic variants (SNPs) with the differential responses among individuals. The lack of adverse effects on liver profile and oxidation markers, together with adequate tolerability and safe immunological adaptations, provide high-quality information for the potential use of CR as co-adjuvant of therapeutic strategies against colorectal cancer.

Graphical Abstract

1. Introduction

Nowadays, there is an increasing interest in nutritional interventions because of their potential role in preventing colon tumorigenesis. Colorectal cancer is the third most prevalent cancer type in males and the second in females worldwide, displaying an increasing incidence even in traditionally low-risk countries such as Spain. In the last years, intervention strategies boosting the activity of the immune system are preferable, together with vehiculation of bioactive components with proven anticancer efficacy. Promotion of proper intestinal mucosal barrier function, cellular defense, and local or systemic inflammation are the main targets of these interventions [1]. Impairment of the regulatory processes of these aspects represent a potential risk of enhanced growth and metastatic ability of cancer cells due to immune dysfunction [2,3]. In the clinical setting, formulas based on different nutrients (i.e., arginine, glutamine, nucleotides, and/or ω-3 fatty acids) have demonstrated benefits in surgical patients by mean of the reduction of inflammatory responses and/or infections [4,5,6]. These beneficial effects have been recognized by ESPEN (European Society for Clinical Nutrition and Metabolism) [7] and ASPEN (American Society of Parenteral and Enteral Nutrition) [8]. A meta-analysis has also demonstrated their effectiveness in diminishing inflammation and enhancing host immunity, although benefits in the clinical outcome were inconsistent [9].
In this scenario, nutritional bioactive co-adjuvants can also help by influencing innate immune signals, which stem at the intestinal level to drive an adequate maturation of adaptive immunity. Nutritional interventions, in combination with classic pharmacological cancer treatments, represent a promising strategy to increase the number of patients which could respond successfully against immune–metabolic-based diseases derived from chemotherapy. The clinical seek for agents with favorable tolerability and safety profiles which could be incorporated as co-adjuvant players to pharmacological anticancer treatments makes also attractive the bioactive supplements. In cancer patients, activation of type-1 innate and acquired immunities are crucial for tumor immunotherapy to overcome immunosuppression [10]. Immune response(s) depends on a multitude of factors such as the offending agent, the microenvironment of the effector cells, and the host’s capacity to respond. Moreover, it must constantly be adaptive and be able to integrate responses from many different cell types. The first line of defense of the innate immune system is the gut mucosa, which not only acts as the primary barrier, but it can also drive adequate adaptive immune responses by influencing its maturation processes. In this sense, dendritic cells (DCs), monocytes (MC), macrophages (Mθ), and other antigen-presenting cells (APCs) are upstream of T helper (Th) responses. Thus, understanding how naturally-occurring preventive/therapeutic agents selectively influence tolerance and immune cell function can set the stage for a new generation of active cancer immunotherapies. The latter could represent a path forward to develop durable, long-lasting immune responses against cancer [11].
Herein, we have investigated the immune-nutritional tolerability and safety of a supplement based on a Rosmarinus officinalis supercritical fluid extract enriched in phenolic diterpenes (RE) and shark liver oil enriched in alkylglycerols (AKG) as a bioactive lipid vehicle [12] to be used as a potential co-adjuvant in intervention strategies for specific patients suffering from colorectal cancer and/or immune disorders.

2. Study Design and Materials and Methods

2.1. Supplement Composition

As we have previously demonstrated, specific ranges of carnosic acid/carnosol composition of supercritical fluid Rosmarinus officinalis L extracts (RE) display relevant antitumor effects in colon and breast cancers [13]. Moreover, RE synergizes with 5-Fu and sensitizes 5-Fu-resistant colon cancer cells to this drug [14]. However, due to the low bioavailability of the bioactive compounds present in plants [15], and with the aim to obtain therapeutic benefits of these compounds in the clinical setting, we have developed a supplement of RE, approved for human use, together with a vehicle system based on bioactive alkylglycerols (AKGs), also approved for human use (PCT/ES2017/070263). The current strategy was developed to enhance bioavailability of RE and to potentiate the antitumor effects of RE with bioactive alkylglycerols. For this, 60 healthy volunteers were randomized to participate in the study: 30 received the study soft gelatin capsules (CR)—rosemary supercritical extract (11.25 mg of diterpene phenols) and shark liver oil (SLO) enriched in AKGs (150 mg)—and another 30 volunteers received the control capsules (CC). Detailed composition of the study (CR) and control capsules (CC) are indicated in Supplementary Table S1. The study and control capsules had the same coating and appearance, and capsules were consumed once a day for six weeks. The randomization was performed blindly. The company that produced the product was the only one that knew which group corresponded to each product

2.2. Subjects and Study Design

The set-up of the clinical trial consisted of a six week, double-blind, randomized and parallel pilot study with two study arms—RE and alkylglycerol containing capsules (CR) and control capsules (CC)—to evaluate the immunomodulatory effect (by quantifying changes in peripheral blood leukocyte subpopulations and the cytokine profile produced by peripheral blood mononuclear cells (PBMCs) after the ex vivo lipopolysaccharide (LPS) stimulation) (primary outcome), together with the modification of oxidation markers and modulation of the expression of genes related to immune responses, inflammation, oxidative stress, lipid metabolism, and cancer in healthy volunteers. Moreover, association of selected genetic variants (SNPs) with the responses to the intervention were also analyzed (secondary outcomes).
The pilot study protocol was approved by the local Ethics Committee of the IMDEA Food Foundation (IMD PI017), and it was carried out in accordance with The Code of Ethics of The World Medical Association (Declaration of Helsinki). Written informed consent was obtained by all subjects prior to starting the trial. Clinical Trial Registry number: This trial was registered at clinical trials.gov as NCT03492086 http://clinicaltrials.gov/. PATENT: (PCT/ES2017/070263).
Volunteers were recruited from the Campus of Cantoblanco and Parque Científico de Madrid, from the GENYAL Platform, and from the dissemination of the study through the IMDEA Food Foundation web page and other means of communication.
Of the 123 volunteers interested in participating in the study, 63 were excluded before randomization for not meeting the inclusion–exclusion criteria of the study, not attending the screening visits (V0), or because they decided not to participate due to distance or personal reasons. Inclusion criteria: Age between 18 and 55 years; adequate understanding of the study; and willingness to complete the entire treatment. Exclusion criteria: Body mass index (BMI) >30 kg/m2; diagnosis of diabetes mellitus (T2D), hypertension, dyslipidemia, or other cardiometabolic disorders; impaired cognitive function; diagnosed hepatic, renal, or cardiovascular disease; subjects with primary immunodeficiency disorders, consumption of drugs with influence on the immune system, or splenectomy; presence of other pathologies like asthma, food allergies, Crohn’s, myasthenia gravis, or lupus; consumption of vitamins, minerals, supplements of antioxidant extracts, or protein supplements in the two weeks prior to the start of the study; subjects treated with drugs affecting the lipid or glycemic profile (i.e., statins, fibrates, diuretics, corticoids, anti-inflammatories, hypoglycemic agents, or insulin) during the previous 30 days; consumption of anticoagulants or antiplatelet agents, cyclosporine, acetylsalicylic acid, antihistamines, or sedatives; hypersensitivity to rosemary, to its components or other members of the family of lipped plants, or to soybean as excipient of the capsules; allergy or hypersensitivity to fish; habitual smoking or high consumption of alcohol; pregnant or lactating women; and high-intensity physical exercise.
Due to the lack of previous studies carried out with the products under study, a pilot study was carried out with a sample size of 60 volunteers (control group–CC: 18 females and 12 males; and study group–CR: 18 females and 12 males). From this study, we intended to calculate the sample size of subsequent studies, taking into account the main efficacy variable—with a final 60 subjects. The randomization procedure was provided by the Biostatistics Unit of the IMDEA Food Foundation. It was carried out using a randomization table generated by Statistic Software R 2.15 version (www.r-project.org, University of Auckland, Auckland, New Zealand.)
The randomization and masking data were strictly confidential. Once the study was completed, the data file was verified and, after performing the results analysis, the codes were opened to proceed with the interpretation of the same. All volunteers registered a complete nutritional survey: 3-day habitual food registry and questionnaires about physical activity, capsule tolerance, and capsule consumption.
Before the first study visit, a group CC volunteer decided to retire due to having problems to attend to, and another volunteer had to be withdrawn due to increased levels of liver enzymes in the baseline analytical group CC, so that the study was finally initiated by 58 volunteers (35 females and 23 males): 28 in group CC (17 females and 11 males) and 30 in group CR (18 females and 12 males). Throughout the study, one volunteer was withdrawn as a precaution since she had persistent headaches, which turned out to be in group CC, so 57 participants finished (27 in group CC and 30 in group CR). No exclusion was made prior to the analysis, since all volunteers registered a consumption ≥75% of the agreement.
Blood samples were drawn and anthropometric data, bioelectric impedance analysis (BIA), and vital constants were monitored at each visit: V1 (baseline parameters), V2 (control visit), and V3 (final parameters). Morning urine was also collected (V1, V3).
Primary outcome: The immunomodulatory effect of the intervention was assessed by analyzing: (i) The levels of cytokines (IL-1β, IL-2, -4, -5, -10, IFNγ, TNFα) produced by PBMCs after ex vivo LPS stimulation; (ii) the effect on the hematological profile and differential counts of leukocytes; and (iii) the phenotyping of leukocyte subpopulations. Secondary outcomes: (iv) Modification of parameters related to inflammation and oxidative stress; and (v) modulation of the expression of genes associated with immune and inflammatory responses, lipid profile, oxidative stress markers, and cancer (from a panel of 48 specifically-selected genes). Single nucleotide polymorphism (SNP) analysis was included to identify putative associations with the differential responses to the intervention. Figure 1A shows a diagram of the study and Figure 1B displays the study flowchart and main objectives.

2.3. Biochemical Analyses

Triacylglycerol (TG), total cholesterol (TC), low-density (LDL-C) and high-density (HDL-C) cholesterol concentrations were measured by routine laboratory (Laboratorio CQS, Madrid, Spain, which follows the UNE-ISO 15189:2007 directives) methods. Urea, creatinine, hepatic enzymes—glutamyl oxaloacetic transaminase (GOT), glutamic-pyruvatetransaminase (GPT), and gamma glutamyltransferase (GGT)—bilirrubin, and alkaline phosphatase (AP) were also determined. The concentration of LDLox was measured by sandwich enzyme-linked immunosorbent assay (ELISA) by using the monoclonal antibody mAb-4E6 (Mercodia AB, Uppsala, Sweden). Urinary thromboxane B2 (TBX2) and total isoprostanes were quantified by competitive ELISA (Enzo Biochem, Inc., New York, USA and Oxford Biomedical Research, Michigan, USA, respectively).

2.4. Isolation of Peripheral Blood Mononuclear Cells (PBMCs)

Subjects were instructed to fast overnight before each blood collection. Blood samples were collected in heparinized tubes (BD Vacutainer, Franklin Lakes, NJ, USA) at each visit between 08:00 and 10:00 to minimize circadian variations, which processed within 2 h of collection and used for PBMCs isolation. Isolation was carried out under sterile conditions to avoid monocytes activation. Briefly, whole blood was diluted (1:1) with phosphate buffer solution (PBS) and centrifuged by density gradient with Histopaque-1077 (Sigma–Aldrich, Madrid, Spain), according to the manufacturer’s instructions. After collection, PBMCs were washed twice with PBS.

2.5. Multiplex Bead Immunoassay

For the analysis of the ex vivo cytokine profile produced by PBMCs after LPS stimulation, isolated PBMC were first incubated for 12 h and then treated with LPS. Supernatants were recovered to determine concentrations of IL-1β, IL-2, IL-4, IL-5, IL-6, IL-8, IL-10, IFNγ, and TNFα using a magnetic bead-based immunoassay (Human High Sensitivity T Cell Magnetic Bead Panel A MAGPIX-Luminex) kit from Millipore Iberica SA (Madrid, 28050, Spain), following the manufacturer’s instructions. A minimum of 50 beads per parameter were analyzed by the MAGPIX-Luminex system (Millipore Iberica SA (Madrid, 28050, Spain). Raw data (median fluorescence intensity, MFI) were analyzed with the xPONENT software 4.1 (Luminex system, Millipore Iberica SA (Madrid, 28050, Spain).

2.6. FACS (Fluorescence Activated Cell Sorting) Analysis of Cluster of Differentiation (CDs) Markers

Analysis of cluster differentiation markers was achieved by flow cytometry with a Gallios analytical flow cytometer (Beckman Coulter, Irving, TX, USA). Aliquots (100 µL) of peripheral blood were mixed with the following fluorochrome-conjugated antibodies: Anti-human CD3-APC, CD4-APCeF750, CD8-PC7, CD14-PC5.5, CD16-PB, CD45-KrO, CD54-FITC, and CD56-PE (Beckman Coulter, Irving, TX, USA), and lysed with the Beckman Coulter Versalyse solution following the manufacturer’s recommendations. To obtain the cell counts for the different populations 100 µL of Flow-Count Fluorospheres (Beckman Coulter, Irving, TX, USA) were added.

2.7. Cell Culture and In Vitro Treatments

Human colon cancer cells SW620 were obtained from American Type Culture Collection (ATCC, Manassas, VA, USA). SW620 cells were cultured in DMEM supplemented with 10% of fetal bovine serum (FBS) and 2 mmol/L glutamine. Cells were maintained under standard conditions of temperature (37 °C), humidity (95%), and carbon dioxide (5%).
For the in vitro validation of the molecular targets found in the clinical trial, we treated SW-620 CRC cells with three different doses corresponding to 0.5 × IC50, 1 × IC50 and 2 × IC50 for 4, being the IC50 value in this cell line of 36.46 ± 7.23 mg/mL as previously described [14].

2.8. RNA Extraction

Total RNA was extracted and purified from homogenized PBMCs with miRNeasy minikit, (Qiagen, Valencia, CA, USA) following the manufacturer’s protocol. Recovered RNA’s concentration and integrity were verified using a Nanodrop ND-1000 spectrophotometer (Nanodrop Technology, Cambridge, UK).

2.9. Gene Expression Analysis

After DNAse I treatment (Invitrogen, Madrid, Spain), 1 g of total RNA was reverse transcribed for 2 h at 37 °C with miScript® II Reverse Transcription kit (Qiagen, Germantown, MD, USA), according to the manufacturer’s guidelines. A Taq-Man Low Density Array (Applied Biosystems, Madrid, Spain) was specifically designed for this experiment, including 47 selected genes (Supplementary Table S2) linked to immune system, inflammation, oxidative stress, lipid metabolism, and cancer-associated genes. Gene-expression assays were performed in a HT–7900 Fast Real-Time PCR. GAPDH was used as endogenous control. RT-StatMiner software (Integromics® Inc., Madison, WI, USA) was used to detect and determine the quality control and differential expression analyses. The Expression Suite Software (Life Technologies, Madrid, Spain) program was used to obtain the Ct data. The ΔCt (Ct gene-Ct GAPDH) was calculated, and then the relative expression (RQ) between visits was calculated (V3−V1) following the 2−ΔΔCt method [16].

2.10. DNA Extraction and Genotyping

Genomic DNA from each participant was isolated from 300 μL of total blood using the QIAamp DNA Blood Mini Kit (Qiagen Sciences, Inc, Germantown, MD, USA) and recovered in 100 μL of nuclease-free water. Concentration and quality were then measured in a nanodrop ND-2000 spectrophotometer (ThermoScientific, Waltham, MA, USA). Genotyping was performed using the QuantStudio 12 K Flex Real-Time PCR System (Life Technologies Inc., Carlsbad, CA, USA) with a TaqMan OpenArray plates. Sixty-four single nucleotide polymorphisms (SNPs) were selected on the basis of their known involvement in different parts of the pathogenic processes of inflammation, immune system, obesity, lipid metabolism, redox homeostasis, and cancer (Supplementary Table S3) and its related phenotypes. The results were analyzed using the TaqMan Genotyper software (Life Technologies Inc., Carlsbad, CA, USA).

2.11. Statistical Analyses

Data were analyzed with the statistical program R Statistical Software version 2.15 (www.r-project.org, University of Auckland, Auckland, New Zealand). Description of qualitative data was made by absolute frequencies and percentages and quantitative data by mean and standard deviation (SD), standard error of the mean (SEM), or 95% confidence intervals (CI), depending on distribution of the data. For comparison of quantitative data between the two groups, Student t-test (parametric test) or Mann-Whitney U or Wilcoxon test (non-parametric test) was used. Chi-square test or Fisher’s test were used to compare proportions of qualitative data between the two groups. Two-way repeated-measures ANOVA was used to evaluate differences of the effect of time (visits), treatment (group), and interaction time x treatment. The p-value was adjusted by sex, age, and BMI. When significant interaction was observed, post-hoc analysis and the Bonferroni correction were applied.
For the statistical analysis of genetic data, X2 test was used to establish the deviation of the genotypic frequencies in controls compared with those expected under the Hardy–Weinberg equilibrium. Three-way ANOVA was used to evaluate the interaction time x treatment x genotype.
All the statistical tests were considered bilateral and a p-value < 0.05 as significant.

3. Results

3.1. Baseline Variables and Safety and Tolerability of Product Consumption

To discard differences after randomization between the control (CC) and intervention (CR) groups and to establish the variable basal levels in the study (V1), the following were monitored: (i) hemogram—leucocytes, hemoglobin, erythrocytes, platelets; (ii) lipid profile—total cholesterol (TC), high density lipoprotein (HDL), low density lipoprotein (LDL), triglycerides (TG); (iii) oxidation-inflammation markers—C reactive protein (CRP), LDLox, isoprostanes, tromboxanes, IFNγ, IL-10, IL-1β, IL-2, IL-4, IL-5, TNFβ-; (iv) anthropometry—BMI, weight, hip and waist circumferences, total fat mass (TFM), total muscle mass (TMM); (v) vital constants—systolic and diastolic blood pressures (SBP, DBP); and (vi) heart rate (HR). No significant differences were found between groups in V1 in any of the variables analyzed, indicating proper randomization (Table 1).
To monitor safety of the study, the levels of liver enzymes (GOT, GPT, GGT, AP), total bilirubin, and creatinine were evaluated at V1, V2, and V3 visits. No significant differences were found in any of the variables, nor in the evolution of hematological profile and vital constants, according to the interaction ‘time x treatment’ (Table 2).
The evolution of inflammation and oxidation parameters (total isoprostanes-Tisos, thromboxanes, TBX2; C reactive protein, CRP; and low density lipoprotein oxidized, LDLox) remained unaltered, except for the levels of LDLox (Table 3A). The post-hoc analysis indicates that LDLox levels in the control group were reduced between visits (p = 0.0289). Besides, a non-significant (p = 0.193) trend was observed to increase LDLox levels in the study group.
The evolution of the lipid profile (total cholesterol, CT; low density lipoprotein, LDL; high density lipoprotein, HDL; and triglycerides, TG) revealed differences in the levels of LDL (Table 3B). The post-hoc analysis showed a significant reduction in the levels of total LDL (p = 0.001) only in the control group (CC). The study group (CR) showed a reduction tendency without reaching levels of significance (p = 0.191).
Evolution of anthropometric parameters and body composition (weight; body mass index, BMI; percentage of fat mass, FM; percentage of muscle mass, MM) was also analyzed. For the percentage of FM, the differences were found at the significance limit (p = 0.0523) (Table 3C). The post-hoc analysis indicates that while the control group (CC) significantly reduced the %FM between visits (p = 0.0016), significant changes could not be quantified for the %FM (p = 0.6983) in the study group (CR). All in all, these results indicate that the observed reduction in LDL and LDLox levels in the control group (CC) seem to be related to the significant reduction in the percentage of fat mass.
During the intervention period, good tolerance to the product was confirmed by registering a questionnaire about the symptoms related to heartburn, diarrhea, nausea, constipation, swollen belly, bad breath, repeated taste, diuretic effect, and presence of severe gastrointestinal symptoms.
There were no observed differences in the study group (CR) during the time of any of the analyzed variables—liver profile: GOT, GPT, GGT, AP, Bil; creatinine; vital constants; oxidation markers: LDLox in blood and isoprostanes and thromboxanes in urine—supporting safety and tolerability of the product.

3.2. Immunomodulatory Effect on Leukocyte Subpopulations

Absolute counts in CD markers—CD45, CD3, CD4, CD8, CD14, CD16, CD56—revealed differences according to the interaction “treatment (CC or CR) x time” (Table 4). The post-hoc analysis indicates a significant increase in CD45+ counts between visits (p = 0.0008) in the study group (CR), meanwhile the control group (CC) had a tendency towards reduction, although not significant (p = 0.7859). CD14+, CD16+, CD56+ significantly increased between visits (V1–V3) in the study group (CR) for CD14+ (p = 0.0001); CD16+ (p = 0.0209) and CD56+ (p = 0.0229). Besides, no significant differences were found in the control group (CC) (p = 0.9636, p = 0.3363, and p = 0.6221, respectively). CD3+, CD4+ and CD8+ cell counts significantly increased in the study group (CR) between visits (p = 0.0001, p = 0.0001, and p = 0.0002, respectively). However, differences were found in the control group (CC) (p = 0.2292, p = 0.1479, and p = 0.4516, respectively).
Immunophenotyping of leukocyte subpopulations, by combined analysis of CD markers, allowed us to identify monocytes (CD56+ CD14+ CD16+ lowCD8+), monocyte/DC-like (CD56 CD4+ CD14+), NK cells (CD56+ CD4), T helper (Th) lymphocytes (CD4+ CD3+) and cytotoxic T lymphocytes (Tcit) (CD8+ CD3+) (Figure 2). Significant changes were found between the treatments: Monocytes (CD56+ CD14+ CD16+ lowCD8+) (p = 0.02) and cytotoxic T lymphocytes (CD8+ CD3+) (p = 0.02) increased in the CR group, while T helper (CD4+ CD3+) decreased (p = 0.02) compared to the CC group. Overall, these results evidence a polarization of the immune adaptive branch towards a phagocytic and effector cytotoxic response. Also, bearing in mind the key role of CD4+ T-helper cells in the generation of CD8+ T-cell memory subsets, the cellular response observed supports a non-specific effector response.

3.3. Cytokine Profile After Ex Vivo LPS Stimulation of PBMC

Further studies have tried to evaluate variations in the cellular response capacity to lipopolysaccharide (LPS) from Gram negative bacteria, which potentially occurs from loss of intestinal barrier caused by colon cancer malignancy [17]. For this reason, we compared the cellular response to the ex vivo LPS stimulation of PBMCs from CC and CR groups. Paired comparisons of the relative variation of each cytokine as a function of the treatment applied were performed (Figure 3). The analysis revealed an increased cellular capacity in response to the LPS stimulation regardless IFNγ (4-fold, p = 0.027) and IL-6 (1.9-fold, p = 0.011), while decreased IL-8 (6.6-fold, p = 0.001) levels of the CR group compared to control group. These observations are in accordance with the development of adequate tolerance responses to LPS, which has been previously reported to occur in coordination with the production of IFNγ [18].
To explore the potential relationship of the cytokine profile with the treatment applied, principal components analysis (PCA) plots were also used for the concurrent display of changes in cell cultures subjected to the treatment with CR or CC. The PCA allows to reduce the number of variables and to identify patterns in the data. By expressing the data in such a way, this analysis highlights the similarities and differences of the data according to the treatment (Supplementary Figure S1). PCA analysis, where two components extracted account for 98.3% of the variance, indicates perturbations in the cytokines secreted by PBMCs cultures obtained from patients receiving CR. The two-dimensional chart of PC1 (Eigenvalue, 1.21) and PC2 (Eigenvalue, 0.79) shows the variability of the observations and how the two dimensions change, suggesting the influence of CR on the phenotypic adaptation responses of PBMCs. Collectively, the ex vivo responses to LPS stimulation of the PBMC cultures established from healthy donors that received CR or CC product, suggest that the likely physiologic impact of CR product on immunocompetent cells can be extended successfully to patients and may substantially potentiate immune response(s) in vivo.
The trend in the effector CD4+ CD3+ cell population associated to the cytokine profile in the individuals administered with CR does not allow us to assume increases in the population of regulatory T cells—Treg—with immunosuppressive character which, generally, in an immature state, are associated with tumors. Moreover, the production of IFNγ favoring a Th1 cellular response usually produces cytotoxic macrophages and effector cytotoxic T cells that would suppose a potential favorable response against tumor development.

3.4. Modulation of the Expression of Genes Related to Immunomodulation, Inflammation, Oxidative Stress, and Cancer

Inflammation and immunity coexist at different steps throughout the progression of tumors. Thus, chronic inflammation has been linked to various stages of the tumorigenic process, including proliferation, transformation, invasion, angiogenesis, and/or metastasis. In the tumor microenvironment, the distribution of the immune system subpopulations, together with cytokines and chemokines, dictates the direction of the response towards a protective inflammation or towards tumor promotion. Here, reactive oxygen species (ROS) can also contribute to the development and progression of numerous diseases, including cancer, by modifying DNA, proteins, and/or lipids, and thus increasing the risk of mutagenesis. Moreover, recent studies emphasize obesity as a risk factor in the development of cancer. For these reasons, genes related to inflammation and inflammatory response, immunomodulation, oxidative stress and antioxidant response, fatty acid metabolism, adipogenesis, obesity, and molecular targets of rosemary in cancer are represented in this study. Supplementary Table S2 summarizes the pathways and genes analyzed in this study.
Gene expression analyses show differences in the evolution between groups as a result of the treatment, except JAK1, NFE2L2, CHKA, and BMP2 genes, which significantly diminished in the CR group (Figure 4A). However, the healthy nature of donors could bias these effects. For this reason, we performed in vitro studies in colon cancer cells (SW-620) treated with RE to confirm and validate the specific effect of RE on decreasing the expression levels of JAK1, NFE2L2, CHKA, and BMP2 genes. As shown in Figure 4B, RE diminished the expression levels of these genes, supporting its usefulness as co-adjuvant in colon cancer therapy (further discussed).

3.5. Associations Between Different Genetic Variants (Single Nucleotide Polymorphisms, SNPs) and Responses to the Dietary Intervention with the Capsules

Many studies with formulas to enhance immune system may mask their effects due to the fact of grouping different types of patients. For this purpose, we selected 64 single nucleotide polymorphisms (SNPs) related to inflammation, immune system, obesity, lipid metabolism, redox homeostasis, and cancer to address putative associations between responses to the study variables. All polymorphisms were in Hardy–Weinberg equilibrium except for rs662799, rs5443, and rs6131 SNPs (Supplementary Table S3). The associations between SNP genotypes and changes in the study variables regarding the interaction treatment (CC/CR) x visit were analyzed. Given the small sample size of individuals in the study, the evaluation of the effect of the SNPs genotypes on the evolution of the variables was examined using the dominant model. Figure 5 shows the parameters with a significant interaction (p < 0.05), or in the limit of the significance between treatment-gene.
Study group carriers homozygous for the major allele (CC) for IFNγ rs2069727 SNP displayed significant reduced levels of IL-5 (p = 0.000361) (Figure 5A). IFNγ rs2069727 SNP could be involved in immune responses mediated by macrophage activation. Study group volunteers genotyped AG+AA (carriers of the minor allele) for ALOX5 rs7913948 SNP displayed significant reduced levels of IL-5 (p = 0.004827) (Figure 5B). ALOX5 rs7913948 is involved in the synthesis of pro-inflammatory and anti-inflammatory molecules. A significant interaction was found between FABP2 rs1799883 SNP and absolute CD45+ counts. The post-hoc analysis indicates that study group carriers of the minor allele (CT, TT) for FABP2 rs1799883 SNP displayed significant increased numbers in absolute cell counts for CD45+ marker (p = 0.01885) (Figure 5C). This SNP has been described to be involved in the intracellular transport of long-chain fatty acids, and thus, it may affect cell growth and proliferation of PBMC positive for CD45+ marker.
Finally, an interaction close to the limit of significance (p = 0.0839) was found between ESR1 rs2234693 and absolute CD14+ counts. The post-hoc analysis revealed that CR group carriers of the minor allele (CT, CC) for ESR1 rs2234693 SNP displayed significant increased numbers in absolute cell counts for CD14+ marker (p < 0.00001) (Figure 5D).

4. Discussion

This study demonstrates the safety and tolerability of a dietary supplement (CR) based on a supercritical rosemary extract (RE) and alkylglycerols used as bioactive vehicle. These safety features of the dietary supplement open new avenues to its potential use as co-adjutant in different diseases. Rosemary extracts (RE) have been shown to exhibit antitumoral effects, both in vitro and in vivo, in a wide range of cancer types [13,14]. Individual components of RE, such as carnosic acid, carnosol, urosolic acid, and rosmarinic acid, have been described as bioactive agents, although the complete extract displays higher bioactivities due to the synergistic effects of the different compounds. Rosemary extracts synergize with several antitumor agents [14] used in clinics, frequently not only enhancing their effects, but also inhibiting their metabolism and reducing the appearance of resistance [14]. Notably, RE is already approved as an agent for human health by the FDA and the EFSA.
Solubility in the physiological media has been identified as a major drawback for bioactive compounds, which failed to display beneficial effects in clinical studies. Herein, shark liver oil (SLO) enriched in bioactive AKGs (i.e., stimulation of hematopoiesis and immunological defenses and antitumor and anti-metastatic activities [19,20,21,22]) was selected as a bioactive lipid carrier of rosemary extract. AKGs were shown effective for improving bioavailability and effectivity of the extract [23]. The anticancer effects of AKGs have been associated with the recruitment and activation of macrophages as a primary antitumor defense (innate immune system response) and to the increased production of cytokines such as IL-12 and IFN-γ, enhancing the cytotoxic activity of NK cells and cytotoxic T lymphocytes, together with an increase proliferation of pre-activated T cells and NK cells [24]. These effects could be reflected in the significant increase of the monocyte, as well as cytotoxic T cell population (Figure 2). This cellular adaptation supports the ex vivo increased IFN γ production by established PBMC primary cultures from patients receiving the CR (Figure 3).
These effects can synergize with the TLR4 inhibition by carnosic acid [25] from rosemary extract when used as complementary therapeutics. Also, the ‘redox’ control of TLR3 activation by antioxidants from rosemary can modulate an adequate immune system response towards a monocyte immunosurveillance phenotype [26]. The molecular signaling mediated by TLR3 involves STAT3 inducing a Th1 type immune response and a MθM1 polarization of the CD56CD4+CD14+ population. Although there are controversial results regardless the role of AKG or shark oil containing AKG in the promotion [27] or not [28] of a MθM1 phenotype, our results seem to indicate that the combination of RE with AKG could contribute to a more balanced M1/M2 of the Mθ population [29]. These effects could involve the induction of IFN-γ production from naive and activated T and NK cells, enhancing the cytotoxic activity of NK cells and cytotoxic T lymphocytes (Figure 3) favoring host defense and protection. Similarly, the association of IFN-γ with the CR group of treatment could predispose immune response(s) towards a more tough control of cell growth [30,31]. The observed reduction in the expression of JAK1 supports the potential of CR to attenuate local inflammation in low grade inflammatory diseases (i.e., metabolic syndrome, obesity, and type 2 diabetes) that could predispose to cancer [32,33].
In addition, NFE2L2 expression levels significantly diminished in the study group (CR) (p = 0.031) compared to control group. Nrf2 is a transcription factor activated by oxidative stress, leading to the expression of antioxidant and detoxifying enzymes and drug transporters (MRPs), as well as anti-apoptotic proteins and proteasomes. Nevertheless, several reports support a role of Nrf2 in oncogenesis. Oncogenes K-Ras, B-Raf, and Myc mediate the transcription and amplification of Nrf2 in cancer cells, contributing to the oncogenic process by diminishing ROS levels [34]. Thus, the observed downregulation in the study group (CR) suggests an overall improvement towards a reduction of oxidative stress and/or inflammation. Moreover, in patients with cancer, downregulation of Nrf2 is expected to contribute to the reduction of drug resistance and an overall therapy improvement.
BMP2 expression levels significantly decreased in the study group (CR) (p = 0.008) compared to control (CC). BMP2 has recently been described to contribute to the promotion and development of colon cancer stem cells, together with an increase in drug resistance. Thus, the nutritional intervention based of RE-AKG (CR) suggests a positive influence on controlling the pool of colon cancer stem cells (CSC) [35].
CHKA expression levels also decreased in the study group (CR), although not significantly (p = 0.083). Chokα is a well described oncogene which stimulates proliferation, transformation, invasion, and metastasis, and it also potentiates tumorigenesis of other oncogenes such as RhoA [36]. It is overexpressed in different human tumors, such as breast, lung, prostate, and colon cancer, and it has been proposed a biomarker with prognostic value of the evolution in early stages of non-small-cell lung (NSCL) cancer [37]. Moreover, Chokα activity has been shown to be increased in colon adenocarcinomas relative to normal tissue, and it has been proposed as an important target in cancer therapy [38,39]. Thus, the observed decrease in the study group (CR) compared to control support its benefits in patients with cancer.
Importantly, the observed results in gene expression regulation were also confirmed in SW-620 colon cancer cells treated only with RE, reflecting the specific role of the phytochemicals present in the RE extract.
The association between different genetic variants (single nucleotide polymorphisms, SNPs) and the response to the dietary intervention was analyzed in 64 single nucleotide polymorphisms (SNPs), which have been related to inflammation, immune system, obesity, lipid metabolism, redox homeostasis, and cancer. Given the small sample size to perform the evaluation by genotypes, the results were analyzed by the dominant model.
Study group homozygous carriers for the major allele (CC) of IFNγ rs2069727 SNP, displayed significant reduced levels of IL-5 (p = 0.000361). Study group (CR) volunteers genotyped AG+AA (carriers of the minor allele) for ALOX5 rs7913948 SNP displayed significant reduced levels of IL-5 (p = 0.004827). A significant interaction was also found between FABP2 rs1799883 SNP and absolute CD45+ counts. Study group carriers of the minor allele (CT, TT) for FABP2 rs1799883 SNP displayed significant increased numbers in absolute cell counts for CD45+ marker. To the best of our knowledge, there is a lack of studies in relation to the biological significance of these findings. These results reflect the genetic contribution to the differential responses to this nutritional intervention and warrants further attention and future studies.
One of the main limitations of the study is the age and healthy status of the volunteers, as impaired immune function has been implicated in the declining health and higher incidence of cancer in the elderly [40]. Nevertheless, as shown in this study, the proposed underlying mechanism(s) seems to imply innate immune ‘Toll-like’ receptors (TLRs) driving the polarization of myeloid cells. It has been suggested that TLRs signaling contributes to the pathogenesis of age-related neurodegenerative diseases [41]. Thus, age-related changes compromise the first cellular line of defence, such as neutrophils and macrophages, as well as cytotoxic T lymphocytes activity. In this scenario, the regulatory influence of the developed product increasing a surveillance monocyte population it is not expected to cause ‘deleterious’ health effects. The impaired or improved efficacy of the product will be genetically and physiologically determined by age and physiological status of cancer patients, affecting bone marrow myeloid cell production and immune-senescence. This effect can result similar as that on chemotherapeutic agents, which efficacy or toxicity will be determined by metabolic capacities affected by age. Another limitation of this pilot study is the small sample size, and more studies are needed to confirm the results obtained in this study. These aspects warrant further studies to establish the efficacy of the product in the clinical practice.
In summary, our results support the safety and tolerability of a dietary supplement based on rosemary extract and alkylglycerols as vehicle that is reflected moderate changes in plasmatic leukocyte populations associated to the study intervention.
The lack of immunosuppressive effects, together with a polarized leukocyte population towards phagocytic and effector cytotoxic T cells in response to CR consumption, supports its inclusion in anticancer strategies. A summary of the potential main biological benefits of the RE-AKG-based formula found in this study that encourages its usefulness as co-adjuvant in colon cancer therapy is presented (Figure 6).

5. Conclusions

Herein are investigated the biological and molecular effects of a nutritional supplement based on a Rosmarinus officinalis supercritical extract rich in phenolic diterpenes (RE) and shark liver oil rich alkylglycerols (AKG) as potential adjuvant in the therapy of patients suffering from cancer or immune disorders. 60 healthy volunteers participated in a six week, double-blind, randomized parallel pilot study with two study arms: RE-AKG capsules (CR) and control capsules (CC). Safety, tolerability and biological effect of product consumption was confirmed. Our results suggest positive immuno-nutritional effects by mean of (i) activation of innate immune responses, (ii) polarization of this response towards effector phagocytic and cytotoxic natural killer cells (NK) and (iii) an anti-inflammatory cytokine profile which potentially might sustain the immune response towards effector anti-tumor cells. Gene expression modulation supports its potential usefulness in cancer patients. Defined SNPs (Single Nucleotide Polymorphisms) were also associated to differential responses regardless immunological variables.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6643/11/9/2001/s1, Table S1. Control Capsules (CC) and Study Capsules (CR) composition; Table S2. Table of the selected genes and pathways analyzed in the study; Table S3. Table of the selected SNPs with frequencies and Hardy–Weinberg equilibrium (HWE) analysis; Figure S1. Principal Components Analysis (PCA) plots of the evolution (V1 and V3) of the ex vivo cytokines produced by LPS stimulated PBMCs for CR and CC. Positively identified cytokines with major influence in PCA were IFNγ, and IL-6 and IL8, with a major weight of PC1 (Eigenvalue, 1.21) versus the PC2 (Eigenvalue, 0.79) represented by other cytokines.

Author Contributions

The authors’ responsibilities were as follows: M.G.d.C., V.L.-K., G.R., and A.R.d.M. conceived and designed the study; V.L.-K. collected the biological material; M.G.d.C. and S.M. performed the experiments; M.G.d.C., J.M.L., V.L.-K., and A.R.d.M. interpreted the data. M.G.d.C. wrote the first draft of the manuscript. All authors critically revised the manuscript for important intellectual content and approved the final manuscript.

Funding

This work was supported by the Spanish Ministry of Science (Plan Nacional I+D+i AGL2016-76736-C3), Regional Government of Community of Madrid (P2018/BAA-4343-ALIBIRD2020-CM), Ramón Areces Foundation and EU Structural Funds.

Acknowledgments

We thank all nutritionists and technical staff from GENYAL Platform from IMDEA Food Institute for their support with the development of the clinical trials.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AKGAlkylglycerols
APAlkaline phosphatase
APCAntigen-presenting cell
BilBilirubin
BMIBody mass index
BMP2Bone morphogenetic protein 2
CHKACholine kinase A
CRPC reactive protein
CCControl capsules
CRIntervention capsules
DCDendritic cell
GGTGamma glutamyltransferase
GOTGlutamyl oxaloacetic transaminase
GPTGlutamic-pyruvatetransaminase
HDLHigh-density lipoprotein
IFNInterferon
ILInterleukin
JAK1Janus kinase 1
LDLLow-density lipoprotein
LPSLipopolysaccharide
Macrophages
MCMonocytes
NFE2L2Nuclear factor (erythroid-derived 2)-like 2
PBMCPeripheral blood mononuclear cells
PCAPrincipal component analysis
RERosmarinus officinalis L. extract

References

  1. Suchner, U.; Kuhn, K.S.; Furst, P. The scientific basis of immunonutrition. Proc. Nutr. Soc. 2000, 59, 553–563. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Senatore, F.J.; Murali, S.; Dasanu, A.C. Accelerated Premalignant Polyposis and Second Colon Cancers: Incriminating Immunosuppression, Radiotherapy, and Systemic Chemotherapy Agents. J. Gastrointest. Cancer 2016, 47, 152–156. [Google Scholar] [CrossRef] [PubMed]
  3. Watt, S.K.; Hasselbalch, H.C.; Skov, V.; Kjær, L.; Thomassen, M.; Kruse, T.A.; Burton, M.; Gögenur, I. Whole Blood Gene Expression Profiling in patients undergoing colon cancer surgery identifies differential expression of genes involved in immune surveillance, inflammation and carcinogenesis. Surg. Oncol. 2018, 27, 208–215. [Google Scholar] [CrossRef] [PubMed]
  4. Andrews, F.J.; Griffiths, R.D. Glutamine: Essential for immune nutrition in the critically ill. Br. J. Nutr. 2002, 87 (Suppl. S1), S3–S8. [Google Scholar] [CrossRef]
  5. Grimble, G.K.; Westwood, O.M. Nucleotides as immunomodulators in clinical nutrition. Curr. Opin. Clin. Nutr. Metab. Care 2001, 4, 57–64. [Google Scholar] [CrossRef] [PubMed]
  6. Weiss, G.; Meyer, F.; Matthies, B.; Pross, M.; Koenig, W.; Lippert, H. Immunomodulation by perioperative administration of n-3 fatty acids. Br. J. Nutr. 2002, 87 (Suppl. S1), S89–S94. [Google Scholar] [CrossRef]
  7. Arends, J.; Bachmann, P.; Baracos, V.; Barthelemy, N.; Bertz, H.; Bozzetti, F.; Fearon, K.; Hütterer, E.; Isenring, E.; Kaasa, S.; et al. ESPEN guidelines on nutrition in cancer patients. Clin. Nutr. 2017, 36, 11–48. [Google Scholar] [CrossRef] [PubMed]
  8. August, D.A.; Huhmann, M.B.A.S.P.E.N. clinical guidelines: Nutrition support therapy during adult anticancer treatment and in hematopoietic cell transplantation. JPEN J. Parenter. Enter. Nutr. 2009, 33, 472–500. [Google Scholar] [CrossRef]
  9. Song, G.M.; Tian, X.; Liang, H.; Yi, L.J.; Zhou, J.G.; Zeng, Z.; Shuai, T.; Ou, Y.X.; Zhang, L.; Wang, Y. Role of Enteral Immunonutrition in Patients Undergoing Surgery for Gastric Cancer: A Systematic Review and Meta-Analysis of Randomized Controlled Trials. Medicine (Baltim.) 2015, 94, e1311. [Google Scholar] [CrossRef]
  10. Jain, A.; Pasare, C. Innate Control. of Adaptive Immunity: Beyond the Three-Signal Paradigm. J. Immunol. 2017, 198, 3791–3800. [Google Scholar] [CrossRef]
  11. Gaudino, S.J.; Kumar, P. Presenting Cells and T Cells Impacts Intestinal Homeostasis, Bacterial Infections, and Tumorigenesis. Front. Immunol. 2019, 6, 360. [Google Scholar] [CrossRef] [PubMed]
  12. Corzo-Martínez, M.; Vázquez, L.; Arranz-Martínez, P.; Menéndez, N.; Reglero, G.; Torres, C.F. Production of a bioactive lipid-based delivery system from radish liver oil by enzymatic glycerolysis. Food Bioprod. Process. 2017, 100, 311–322. [Google Scholar] [CrossRef]
  13. Gonzalez-Vallinas, M.; Reglero, G.; de Molina, A.R. Rosemary (Rosmarinus officinalis L.) Extract as a Potential Complementary Agent in Anticancer Therapy. Nutr. Cancer 2015, 67, 1221–1229. [Google Scholar] [PubMed]
  14. Gonzalez-Vallinas, M.; Molina, S.; Vicente, G.; de la Cueva, A.; Vargas, T.; Santoyo, S.; García-Risco, M.R.; Fornari, T.; Reglero, G.; Ramírez de Molina, A. Antitumor effect of 5-fluorouracil is enhanced by rosemary extract in both drug sensitive and resistant colon cancer cells. Pharm. Res. 2013, 72, 61–68. [Google Scholar] [CrossRef] [PubMed]
  15. Gonzalez-Vallinas, M.; González-Castejón, M.; Rodríguez-Casado, A.; Ramírez de Molina, A. Dietary phytochemicals in cancer prevention and therapy: A complementary approach with promising perspectives. Nutr. Rev. 2013, 71, 585–599. [Google Scholar] [CrossRef]
  16. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  17. Ahmad, R.; Kumar, B.; Chen, Z.; Chen, X.; Müller, D.; Lele, S.M.; Washington, M.K.; Batra, S.K.; Dhawan, P.; Singh, A.B. Loss of claudin-3 expression induces IL6/gp130/Stat3 signaling to promote colon cancer malignancy by hyperactivating Wnt/beta-catenin signaling. Oncogene 2017, 36, 6592–6604. [Google Scholar] [CrossRef]
  18. Chen, J.; Ivashkiv, L.B. IFN-gamma abrogates endotoxin tolerance by facilitating Toll-like receptor-induced chromatin remodeling. Proc. Natl. Acad. Sci. USA 2010, 107, 19438–19443. [Google Scholar] [CrossRef]
  19. Pugliese, P.T.; Jordan, K.; Cederberg, H.; Brohult, J. Some biological actions of alkylglycerols from shark liver oil. J. Altern. Complement. Med. 1998, 4, 87–99. [Google Scholar] [CrossRef]
  20. Lewkowicz, P.; Banasik, M.; Głowacka, E.; Lewkowicz, N.; Tchórzewski, H. Effect of high doses of shark liver oil supplementation on T cell polarization and peripheral blood polymorphonuclear cell function. Pol. Merkur. Lek. 2005, 18, 686–692. [Google Scholar]
  21. Deniau, A.L.; Mosset, P.; Pédrono, F.; Mitre, R.; Le Bot, D.; Legrand, A.B. Multiple beneficial health effects of natural alkylglycerols from shark liver oil. Mar. Drugs 2010, 8, 2175–2184. [Google Scholar] [CrossRef] [PubMed]
  22. Deniau, A.L.; Mosset, P.; Le Bot, D.; Legrand, A.B. Which alkylglycerols from shark liver oil have anti-tumour activities? Biochimie 2011, 93, 1–3. [Google Scholar] [CrossRef] [PubMed]
  23. Molina, S.; Moran-Valero, M.I.; Martin, D.; Vázquez, L.; Vargas, T.; Torres, C.F.; Ramirez de Molina, A.; Reglero, G. Antiproliferative effect of alkylglycerols as vehicles of butyric acid on colon cancer cells. Chem. Phys. Lipids 2013, 175–176, 50–56. [Google Scholar] [CrossRef] [PubMed]
  24. Iannitti, T.; Palmieri, B. An update on the therapeutic role of alkylglycerols. Mar. Drugs 2010, 8, 2267–2300. [Google Scholar] [CrossRef] [PubMed]
  25. Park, M.Y.; Mun, S.T. Carnosic acid inhibits TLR4-MyD88 signaling pathway in LPS-stimulated 3T3-L1 adipocytes. Nutr. Res. Pract. 2014, 8, 516–520. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Yoshida, S.; Shime, H.; Funami, K.; Takaki, H.; Matsumoto, M.; Kasahara, M.; Seya, T. The Anti-Oxidant Ergothioneine Augments the Immunomodulatory Function of TLR Agonists by Direct Action on Macrophages. PLoS ONE 2017, 12, e0169360. [Google Scholar] [CrossRef] [PubMed]
  27. Yamamoto, N.; St Claire, D.A., Jr.; Homma, S.; Ngwenya, B.Z. Activation of mouse macrophages by alkylglycerols, inflammation products of cancerous tissues. Cancer Res. 1988, 48, 6044–6049. [Google Scholar] [PubMed]
  28. Vitorino, D.C.; Buzzachera, C.F.; Curi, R.; Fernandes, L.C. Effect of chronic supplementation with shark liver oil on immune responses of exercise-trained rats. Eur. J. Appl. Physiol. 2010, 108, 1225–1232. [Google Scholar] [CrossRef]
  29. Mantovani, A.; Sica, A.; Sozzani, S.; Allavena, P.; Vecchi, A.; Locati, M. The chemokine system in diverse forms of macrophage activation and polarization. Trends Immunol. 2004, 25, 677–686. [Google Scholar] [CrossRef]
  30. Ucer, U.; Scheurich, P.; Bartsch, H.; Berkovic, D.; Ertel, C.; Pfizenmaier, K. Specific membrane receptors for human interferon-gamma (IFN-gamma). Behring Inst. Mitt. 1987, 81, 88–97. [Google Scholar]
  31. Cassetta, L.; Pollard, J.W. Cancer immunosurveillance: Role of patrolling monocytes. Cell Res. 2016, 26, 3–4. [Google Scholar] [CrossRef]
  32. Stuart, E.; Scheurich, P.; Bartsch, H.; Berkovic, D.; Ertel, C.; Pfizenmaier, K. Therapeutic inhibition of Jak activity inhibits progression of gastrointestinal tumors in mice. Mol. Cancer 2014, 13, 468–474. [Google Scholar] [CrossRef]
  33. An, H.J.; Choi, E.K.; Kim, J.S.; Hong, S.W.; Moon, J.H.; Shin, J.S.; Ha, S.H.; Kim, K.P.; Hong, Y.S.; Lee, J.L.; et al. INCB018424 induces apoptotic cell death through the suppression of pJAK1 in human colon cancer cells. Neoplasma 2014, 61, 56–62. [Google Scholar] [CrossRef] [Green Version]
  34. DeNicola, G.M.; Karreth, F.A.; Humpton, T.J.; Gopinathan, A.; Wei, C.; Frese, K.; Mangal, D.; Yu, K.H.; Yeo, C.J.; Calhoun, E.S.; et al. Oncogene-induced Nrf2 transcription promotes ROS detoxification and tumorigenesis. Nature 2011, 475, 106–109. [Google Scholar] [CrossRef]
  35. Kim, B.R.; Oh, S.C.; Lee, D.H.; Kim, J.L.; Lee, S.Y.; Kang, M.H.; Lee, S.I.; Kang, S.; Joung, S.Y.; Min, B.W. BMP-2 induces motility and invasiveness by promoting colon cancer stemness through STAT3 activation. Tumour Biol. 2015, 36, 9475–9486. [Google Scholar] [CrossRef]
  36. Ramirez de Molina, A.; Penalva, V.; Lucas, L.; Lacal, J.C. Regulation of choline kinase activity by Ras proteins involves Ral-GDS and PI3K. Oncogene 2002, 21, 937–946. [Google Scholar] [CrossRef]
  37. Ramirez de Molina, A.; Sarmentero-Estrada, J.; Belda-Iniesta, C.; Tarón, M.; Ramírez de Molina, V.; Cejas, P.; Skrzypski, M.; Gallego-Ortega, D.; de Castro, J.; Casado, E.; et al. Expression of choline kinase alpha to predict outcome in patients with early-stage non-small-cell lung cancer: A retrospective study. Lancet Oncol. 2007, 8, 889–897. [Google Scholar] [CrossRef]
  38. Janardhan, S.; Srivani, P.; Sastry, G.N. Choline kinase: An important target for cancer. Curr. Med. Chem. 2006, 13, 1169–1186. [Google Scholar] [CrossRef]
  39. Ramirez de Molina, A.; Gallego-Ortega, D.; Sarmentero-Estrada, J.; Lagares, D.; Gómez Del Pulgar, T.; Bandrés, E.; García-Foncillas, J.; Lacal, J.C. Choline kinase as a link connecting phospholipid metabolism and cell cycle regulation: Implications in cancer therapy. Int. J. Biochem. Cell Biol. 2008, 40, 1753–1763. [Google Scholar] [CrossRef]
  40. Jackaman, C.; Tomay, F.; Duong, L.; Abdol Razak, N.B.; Pixley, F.J.; Metharom, P.; Nelson, D.J. Aging and cancer: The role of macrophages and neutrophils. Ageing Res. Rev. 2017, 36, 105–116. [Google Scholar] [CrossRef]
  41. Calvo-Rodríguez, M.; de la Fuente, C.; García-Durillo, M.; García-Rodríguez, C.; Villalobos, C.; Núñez, L. Aging and amyloid β oligomers enhance TLR4 expression, LPS-induced Ca2+ responses, and neuron cell death in cultured rat hippocampal neurons. J. Neuroinflamm. 2017, 14, 24. [Google Scholar] [CrossRef]
Figure 1. (A) Diagram of the study. (B) Study flowchart, objectives, and analyzed parameters. Abbreviations: CC, control capsules; CR, intervention capsules; BMI, body mass index; PBMC, peripheral blood mononuclear cell; SNPs, single nucleotide polymorphisms.
Figure 1. (A) Diagram of the study. (B) Study flowchart, objectives, and analyzed parameters. Abbreviations: CC, control capsules; CR, intervention capsules; BMI, body mass index; PBMC, peripheral blood mononuclear cell; SNPs, single nucleotide polymorphisms.
Nutrients 11 02001 g001
Figure 2. Combined CD analysis to identify leukocyte subpopulations. The phenotyping of leukocyte subpopulations, by combined analysis of CD markers, allows to differentiate among monocytes (CD56+ CD14+ CD16+ lowCD8+), macrophages (CD56 CD4+ CD14+), NK cells (CD56+ CD4), T helper (Th) lymphocytes (CD4+ CD3+), and cytotoxic T lymphocytes (Tcit) (CD8+ CD3+).
Figure 2. Combined CD analysis to identify leukocyte subpopulations. The phenotyping of leukocyte subpopulations, by combined analysis of CD markers, allows to differentiate among monocytes (CD56+ CD14+ CD16+ lowCD8+), macrophages (CD56 CD4+ CD14+), NK cells (CD56+ CD4), T helper (Th) lymphocytes (CD4+ CD3+), and cytotoxic T lymphocytes (Tcit) (CD8+ CD3+).
Nutrients 11 02001 g002
Figure 3. Paired comparisons of the relative variation of each cytokine as a function of the treatment applied. The analysis revealed an increased cellular capacity in response to the LPS stimulation regardless of IFNγ (4-fold, p = 0.027) and IL-6 (1.9-fold, p = 0.011), while decreased IL-8 (6.6-fold, p = 0.001) levels of the CR group compared to control group.
Figure 3. Paired comparisons of the relative variation of each cytokine as a function of the treatment applied. The analysis revealed an increased cellular capacity in response to the LPS stimulation regardless of IFNγ (4-fold, p = 0.027) and IL-6 (1.9-fold, p = 0.011), while decreased IL-8 (6.6-fold, p = 0.001) levels of the CR group compared to control group.
Nutrients 11 02001 g003
Figure 4. Modulation of the expression of genes related to immunomodulation, inflammation, oxidative stress, and cancer. (A) Genes whose expression were statistically different between groups (CC or CR), regardless of evolution between visits, are shown. JAK1, NFE2L2, and BMP2 were significantly decreased in the study group (CR) compared to control (CC) (p = 0.02, p = 0.031, and p = 0.008, respectively). CHKA was diminished in the study group (CR) between visits, although not significantly (p = 0.083). (B) In vitro validation of the effect of rosemary extract in downregulating the expression of JAK1, NE2L2, BMP2, and CHKA in SW-620 colon cancer cells after treatment with three different doses corresponding to 0.5 × IC50, 1 × IC50, and 2 × IC50 for 4, being the IC50 value in this cell line of 36.46 ± 7.23 mg/mL, as previously described [14]. *, **, *** = p-value < 0.05, 0.01 and 0.005 respectively.
Figure 4. Modulation of the expression of genes related to immunomodulation, inflammation, oxidative stress, and cancer. (A) Genes whose expression were statistically different between groups (CC or CR), regardless of evolution between visits, are shown. JAK1, NFE2L2, and BMP2 were significantly decreased in the study group (CR) compared to control (CC) (p = 0.02, p = 0.031, and p = 0.008, respectively). CHKA was diminished in the study group (CR) between visits, although not significantly (p = 0.083). (B) In vitro validation of the effect of rosemary extract in downregulating the expression of JAK1, NE2L2, BMP2, and CHKA in SW-620 colon cancer cells after treatment with three different doses corresponding to 0.5 × IC50, 1 × IC50, and 2 × IC50 for 4, being the IC50 value in this cell line of 36.46 ± 7.23 mg/mL, as previously described [14]. *, **, *** = p-value < 0.05, 0.01 and 0.005 respectively.
Nutrients 11 02001 g004
Figure 5. Associations between genetic variants (SNPs) and responses to the nutritional intervention. (A) IL-5 levels by genotype and treatment for INFγ rs2069727; (B) IL-5 levels by genotype and treatment for ALOX5 rs7913948, (C) Absolute CD45+ counts by genotype and treatment for FABP2 rs1799883, (D) Absolute CD14+ counts by genotype and treatment for ESR1 rs2234693.
Figure 5. Associations between genetic variants (SNPs) and responses to the nutritional intervention. (A) IL-5 levels by genotype and treatment for INFγ rs2069727; (B) IL-5 levels by genotype and treatment for ALOX5 rs7913948, (C) Absolute CD45+ counts by genotype and treatment for FABP2 rs1799883, (D) Absolute CD14+ counts by genotype and treatment for ESR1 rs2234693.
Nutrients 11 02001 g005
Figure 6. Proposed model of potential main biological benefits of RE-AKG based formula found in this study that encourages its usefulness as co-adjuvant in colon cancer therapy.
Figure 6. Proposed model of potential main biological benefits of RE-AKG based formula found in this study that encourages its usefulness as co-adjuvant in colon cancer therapy.
Nutrients 11 02001 g006
Table 1. Baseline parameters (V1) (validation of randomization).
Table 1. Baseline parameters (V1) (validation of randomization).
Group CC (n = 28)Group CR (n = 30)p
Age28.32 ± 11.3927.5 ± 9.040.963
Hemogram
Leucocytes (10.3/µL)6.42 ± 1.626.33 ± 1.830.856
Hemoglobin (g/dL)13.57 ± 1.213.61 ± 1.430.896
Eritrocytes (10.6/µL)4.52 ± 0.354.57 ± 0.420.669
Platelets (10.3/µL)231.64 ± 42.35236.1 ± 47.380.846
Lipidic Profile
TC (mg/dL)174.21 ± 25.90184.26 ± 29.780.319
HDL (mg/dL)58.64 ± 10.2861.64 ± 9.570.255
LDL (mg/dL)101 ± 26.02108.25 ± 27.000.302
TG (mg/dL)72.82 ± 33.2071.8 ± 31.070.904
Oxidation-inflammation markers
CRPus (mg/dL)0.12 ± 0.150.09 ± 0.010.308
LDL-ox (U/L)36.91 ± 11.0037.09 ± 11.080.95
Isoprostanes (ng/mL)1.06 ± 0.551.23 ± 0.910.853
Tromboxanes (ng/mL)983.05 ± 1149.571193.59 ± 1662.670.783
IFNg (pg/mL)6.80 ± 8.005.27 ± 5.920.619
IL 10 (pg/mL)621.51 ± 1065.24926.88 ± 1692.320.987
IL 1β (pg/mL)1361.55 ± 833.31005.92 ± 825.920.05
IL 2 (pg/mL)0.44 ± 0.240.45 ± 0.220.73
IL 4 (pg/mL)2.95 ± 1.553.19 ± 1.770.363
IL 5 (pg/mL)0.44 ± 0.220.84 ± 1.220.844
TNFα (pg/mL)1940.55 ± 546.881671.57 ± 666.930.107
Anthropometric data
Weight (kg)63.11 ± 8.1065.57 ± 8.310.258
BMI (kg/m2)22.4 ± 2.3422.54 ± 2.190.984
Waist (cm)78.18 ± 7.3380.96 ± 9.960.23
Hip (cm)98.69 ± 5.0197.38 ± 8.000.456
T Fat Mass (%) 27.36 ± 8.1125.6 ± 9.520.455
T Muscle Mass (%) 32.2 ± 6.0533.92 ± 7.080.497
Vital constants
SBP116 ± 12.12118.03 ± 15.350.577
DBP68.79 ± 6.9769.17 ± 8.910.856
Data represents mean ± SD. A p-value > 0.05 was considered non-significant. CC (Control capsules), CR (Intervention capsules), TC (Total cholesterol), HDL (High Density Lipoprotein), LDL (Low Density Lipoprotein), TG (Triglycerydes), BMI (Body Mass Index), SBP (Systolic Blood Pressure), DBP (Diastolic Blood Pressure).
Table 2. Safety and tolerability of product consumption was assessed by monitoring evolution between visits of (A) liver enzymes (glutamyl oxaloacetic transaminase, GOT; glutamic-pyruvatetransaminase, GPT; gamma-glutamyltransferase, GGT; bilirubin; and alkaline phosphatase, AP) and creatinine; (B) vital constants (systolic blood pressure, SBP; diastolic blood pressure, DBP; and heart rate, HR); and (C) evolution of hematologic profile.
Table 2. Safety and tolerability of product consumption was assessed by monitoring evolution between visits of (A) liver enzymes (glutamyl oxaloacetic transaminase, GOT; glutamic-pyruvatetransaminase, GPT; gamma-glutamyltransferase, GGT; bilirubin; and alkaline phosphatase, AP) and creatinine; (B) vital constants (systolic blood pressure, SBP; diastolic blood pressure, DBP; and heart rate, HR); and (C) evolution of hematologic profile.
(A) Liver Enzymes, Total Bilirubin, Creatinine
CC (Control)CR (Study)Effect *
V1V2V3V1V2V3VisitGroupVisit × Group
GOT (UI/L)18.26 (0.73)18.26 (1.05)18.07 (0.75)19.7 (0.86)18.63 (0.84)20.57 (0.9)0.24160.15120.1512
GPT (UI/L)17.74 (1.4)16.15 (1.21)15.7 (1.31)18.83 (1.47)16.2 (1.26)18.63 (1.42)0.00650.3680.0794
GGT (UI/L)17.99 (2.18)16.01 (1.5)15.24 (1.04)18.63 (1.3)17.46 (1.14)18.36 (1.21)0.01870.28670.1182
AP (UI/L)57.07 (4.21)56.19 (3.39)57.07 (3.66)54.2 (3.37)54.83 (2.67)56.97 (2.89)0.26370.7410.4248
TBil (mg/dL)0.77 (0.08)0.7 (0.08)0.71 (0.08)0.79 (0.06)0.7 (0.06)0.65 (0.05)0.0380.94840.4705
Creatinine (mg/dL)0.79 (0.02)-0.83 (0.02)0.8 (0.02)-0.85 (0.02)<0.0010.56610.5346
(B) Vital Constants
CC (Control)CR (Study)Effect *
V1 V3V1 V3VisitGroupVisit × Group
SBP (mmHg)115.85 (2.28) 110.69 (2.31)118.03 (2.8) 115.37 (2.63)0.00910.15880.3909
DBP (mmHg)69 (1.21) 69.88 (1.7)69.17 (1.65) 69.28 (1.64)0.60830.99660.6844
HR (beat/min)71.38 (2.46) 70.08 (2.34)67.21 (2.59) 66.1 (1.94)0.35590.12540.9132
(C) Hematological Profile
CC (Control)CR (Study)Effect *
V1 V3V1 V3VisitGroupVisit × Group
Erythrocytes (10.6/µL)4.53 (0.07) 4.42 (0.07)4.57(0.08) 4.46 (0.08)<0.0010.57530.8106
Erythrocytes (%)40.84 (0.57) 39.73 (0.58)40.95 (0.7) 39.94 (0.67)<0.0010.59270.8402
Hemoglobin (g/dL)13.63 (0.23) 13.47 (0.22)13.61 (0.26) 13.52 (0.25)0.07920.73240.6625
Leukocytes (10.3/µL)6.4 (0.32) 6.01 (0.29)6.33 (0.33) 6.5 (0.32)0.61780.62660.1667
Lymphocytes (10.3/µL)2.15 (0.09) 2.14 (0.1)2.2 (0.09) 2.3 (0.14)0.4080.44490.3247
Neutrophils (10.3/µL)3.53 (0.28) 3.2 (0.22)3.41 (0.28) 3.47 (0.23)0.46450.83030.2574
Monocytes (10.3/µL)0.53 (0.03) 0.5 (0.03)0.51 (0.02) 0.5 (0.02)0.25640.74440.4594
Basophils (10.3/µL)0.06 (0) 0.05 (0)0.07 (0) 0.06 (0)<0.0010.09620.9151
Platelets (10.3/µL)231.33 (8.3) 217.59 (7.78)236.1 (8.65) 228.67 (8.26)0.00480.47150.3886
* Effect: Visit indicates evolution in the same group from V1 to V3; group indicates the differences between groups (CC: control capsules), CR: intervention capsules) regardless of time; and visit × group indicates the differences in evolution (V1 to V3) between groups as a result of the treatment (CC, CR). The p-value of the interaction was adjusted by sex, age, and BMI. Data expressed as mean (SEM). A p-value > 0.05 was considered non-significant.
Table 3. Evolution between visits of (A) inflammation and oxidation parameters (LDLox, total isoprotanes, tromboxanes TBX2, and C-reactive protein); (B) lipid profile (total cholesterol, TC; high density lipoprotein, HDL; triglycerides, TG; and low density lipoprotein, LDL); and (C) anthropometric data (weight; body mass index, BMI; fat mass, FM; muscle mass, MS; and waist circumference, Wc).
Table 3. Evolution between visits of (A) inflammation and oxidation parameters (LDLox, total isoprotanes, tromboxanes TBX2, and C-reactive protein); (B) lipid profile (total cholesterol, TC; high density lipoprotein, HDL; triglycerides, TG; and low density lipoprotein, LDL); and (C) anthropometric data (weight; body mass index, BMI; fat mass, FM; muscle mass, MS; and waist circumference, Wc).
(A) Markers of Inflammation and Oxidation
CC (Control)CR (Study)Effect *
V1V3V1V3VisitGroupVisit × Group
LDLox (U/L)36.91 (2.12)34.15 (1.81)37.09 (2.02)38.66 (2.14)0.59020.36530.0196
T isop (ng/mg creat)0.87 (0.08)0.84 (0.06)0.8 (0.09)0.95 (0.07)0.31290.89480.1045
TXB2 (ng/mg creat)713.1 (128.16)709.39 (120.37)713.74 (178.88)879.44 (187.51)0.55060.73630.4338
CRP (mg/dL)0.11 (0.03)0.09 (0.02)0.09 (0.02)0.07 (0.01)0.08260.17250.7098
(B) Lipid Profile
CC (Control)CR (Study)Effect *
V1V3V1V3VisitGroupVisit × Group
TC (mg/dL)174.96 (5.02)171.05 (5.02)184.26 (5.44)185.77 (5.93)0.60620.07230.1745
HDL (mg/dL)58.41 (2)65.67 (2.65)61.64 (1.75)65.59 (2.7)<0.0010.62350.1491
TG (mg/dL)73.19 (6.5)70.74 (5.9)71.8 (5.67)75.2 (7.73)0.85270.81380.4571
LDL (mg/dL)101.91 (5.02)91.23 (5.29)108.25 (4.93)105.13 (5.6)<0.0010.10780.0434
(C) Anthropometric Data
CC (Control)CR (Study)Effect *
V1V3V1V3VisitGroupVisit × Group
Weight (kg)63.11 (1.59)62.57 (1.58)65.57 (1.52)65.35 (1.58)0.02220.02760.3273
BMI (kg/m2)22.3 (0.45)22.11 (0.42)22.54 (0.4)22.46 (0.39)0.0031<0.0010.3137
FM (%)26.95 (1.57)26.02 (1.56)25.6 (1.74)25.5 (1.71)0.01250.2280.0523
MM (%)32.42 (1.19)32.92 (1.2)33.92 (1.29)33.92 (1.3)0.07270.07610.0787
Wc(cm)78.11 (1.43)77.75 (1.37)80.96 (1.82)79.14 (1.19)0.22670.09540.4164
* Effect: Visit indicates evolution in the same group from V1 to V3; group indicates the differences between groups (CC, CR) regardless of time; and visit × group indicates the differences in evolution (V1 to V3) between groups as a result of the treatment (CC, CR). The p-value of the interaction was adjusted by sex, age, and BMI. Data expressed as mean (SEM). A p-value > 0.05 was considered non-significant.
Table 4. Evolution in absolute counts of cluster differentiation (CD) markers.
Table 4. Evolution in absolute counts of cluster differentiation (CD) markers.
Cluster Differentiation Markers Counts
CC (Control)CR (Study)Effect *
V1V3V1V3VisitGroupVisit × Group
CD45+ (cs/mL)6458.99 (457.29)6231.24 (474.55)6260.18 (403.99)7764.86 (378.69)0.04910.25780.0207
CD3+ (cs/mL)∆1615.72 (124.1)1744.94 (140.34)1581.5 (93.34)2012.5 (161.31)0.00180.55770.0678
CD4+ (cs/mL)∆960.67 (66.18)1045.5 (77.89)902.11 (59.85)1149.33 (89.98)<0.0010.90080.0714
CD8+ (cs/mL) ∆526.94 (68.13)555.56 (63.82)524.28 (36.67)667.11 (62.62)0.00470.50250.0563
CD14+ (cs/mL)473.06 (39.68)471.22 (40.84)409.94 (28.53)565.72 (42.81)0.01440.78580.0126
CD16+ (cs/mL)3485.56 (362.15)3183.78 (310.67)3253.83 (274.68)3935.72 (224.94)0.3650.33930.0318
CD56+ (cs/mL)429.61 (51.53)407.78 (46.1)415.33 (40.05)516.11 (50.95)0.24170.49040.0716
* Effect: Visit indicates the evolution in the same group from the beginning to the end of the intervention; group indicates the differences between groups regardless of time; and visit × group indicates the differences in evolution between groups as a result of the treatment received (CC or CR). The p-value of the interaction was adjusted by sex, age, and BMI. Data expressed as mean (SEM). A p-value > 0.05 was considered non-significant. Analysis performed on a subsample (n = 18) for group CC and (n = 18) for group CR. ∆ Counts on CD45+.

Share and Cite

MDPI and ACS Style

Gómez de Cedrón, M.; Laparra, J.M.; Loria-Kohen, V.; Molina, S.; Moreno-Rubio, J.; Montoya, J.J.; Torres, C.; Casado, E.; Reglero, G.; Ramírez de Molina, A. Tolerability and Safety of a Nutritional Supplement with Potential as Adjuvant in Colorectal Cancer Therapy: A Randomized Trial in Healthy Volunteers. Nutrients 2019, 11, 2001. https://doi.org/10.3390/nu11092001

AMA Style

Gómez de Cedrón M, Laparra JM, Loria-Kohen V, Molina S, Moreno-Rubio J, Montoya JJ, Torres C, Casado E, Reglero G, Ramírez de Molina A. Tolerability and Safety of a Nutritional Supplement with Potential as Adjuvant in Colorectal Cancer Therapy: A Randomized Trial in Healthy Volunteers. Nutrients. 2019; 11(9):2001. https://doi.org/10.3390/nu11092001

Chicago/Turabian Style

Gómez de Cedrón, Marta, José Moises Laparra, Viviana Loria-Kohen, Susana Molina, Juan Moreno-Rubio, Juan Jose Montoya, Carlos Torres, Enrique Casado, Guillermo Reglero, and Ana Ramírez de Molina. 2019. "Tolerability and Safety of a Nutritional Supplement with Potential as Adjuvant in Colorectal Cancer Therapy: A Randomized Trial in Healthy Volunteers" Nutrients 11, no. 9: 2001. https://doi.org/10.3390/nu11092001

APA Style

Gómez de Cedrón, M., Laparra, J. M., Loria-Kohen, V., Molina, S., Moreno-Rubio, J., Montoya, J. J., Torres, C., Casado, E., Reglero, G., & Ramírez de Molina, A. (2019). Tolerability and Safety of a Nutritional Supplement with Potential as Adjuvant in Colorectal Cancer Therapy: A Randomized Trial in Healthy Volunteers. Nutrients, 11(9), 2001. https://doi.org/10.3390/nu11092001

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

Article Metrics

Back to TopTop