Skip to main content

An 8-week freeze-dried blueberry supplement impacts immune-related pathways: a randomized, double-blind placebo-controlled trial



Blueberries contain high levels of polyphenolic compounds with high in vitro antioxidant capacities. Their consumption has been associated with improved vascular and metabolic health.


The objective was to examine the effects of blueberry supplement consumption on metabolic syndrome (MetS) parameters and potential underlying mechanisms of action.


A randomized double-blind placebo-controlled intervention trial was conducted in adults at risk of developing MetS. Participants consumed 50 g daily of either a freeze-dried highbush blueberry powder (BBP) or a placebo powder for 8 weeks (n = 49). MetS phenotypes were assessed at weeks 0, 4 and 8. Fasting blood gene expression profiles and plasma metabolomic profiles were examined at baseline and week 8 to assess metabolic changes occurring in response to the BBP. A per-protocol analysis was used.


A significant treatment effect was observed for plasma triglyceride levels that was no longer significant after further adjustments for age, sex, BMI and baseline values. In addition, the treatment*time interactions were non-significant therefore suggesting that compared with the placebo, BBP had no statistically significant effect on body weight, blood pressure, fasting plasma lipid, insulin and glucose levels, insulin resistance (or sensitivity) or glycated hemoglobin concentrations. There were significant changes in the expression of 49 genes and in the abundance of 35 metabolites following BBP consumption. Differentially regulated genes were clustered in immune-related pathways.


An 8-week BBP intervention did not significantly improve traditional markers of cardiometabolic health in adults at risk of developing MetS. However, changes in gene expression and metabolite abundance suggest that clinically significant cardiometabolic changes could take longer than 8 weeks to present and/or could result from whole blueberry consumption or a higher dosage. BBP may also have an effect on factors such as immunity even within a shorter 8-week timeframe.

Clinical trial registration, NCT03266055, 2017


With an estimated prevalence of almost 35% in the USA [1] and 19% in Canada [2], metabolic syndrome (MetS) is a condition of great concern in developed countries. By clustering abdominal obesity, insulin resistance (IR), dyslipidemia and hypertension, people with MetS are at higher risk of developing type 2 diabetes and cardiovascular disease [3]. Oxidative stress is a commonality for these risk factors and is suspected of playing a pro-inflammatory role thus promoting their development [4]. Oxidative stress occurs when the body’s antioxidant defenses are not sufficient to counterbalance production of free radicals and other reactive oxygen species (ROS) [4]. Free radicals are highly reactive and have the potential to damage essential biomolecules including lipids, proteins and DNA, which in turn impair their functions [5]. Notably, excess calorie intake leads to increased ROS production, which could eventually generate a pro-oxidant environment susceptible to cause cellular dysfunction promoting obesity-related diseases [4].

Nutrition has the potential to help mitigate issues related to increased ROS production. For example, many fruits are naturally rich in polyphenolic molecules with potent antioxidant and anti-inflammatory activity, which are part of their natural defense system against pathogens, environmental fluctuations and pro-oxidant stressors such as ultraviolet light [6]. Berries are among fruits with the highest polyphenol content, mainly anthocyanins, resulting in their bright color [7,8,9]. Compared with 25 other commonly eaten fruits, wild and cultivated blueberries were previously ranked first and fifth, respectively, for their total phenolic content (respectively 429 and 285 gallic acid equivalent per 100 g of fruit) [9]. Analyses of the effect of lowbush blueberry plant extracts on cultured cells also revealed insulin- and glitazone-like activities, suggesting that blueberries may possess antidiabetic properties [8, 10]. Blueberry supplementation has previously been associated with favorable effects on blood pressure [11, 12], some aspects of brain cognitive function [13, 14] and plasma biomarkers of oxidative stress and inflammation [11]. Blueberry consumption has also been demonstrated to be protective against the risk of myocardial infarction [15] and diabetes [16]. However, these results are inconsistent in the literature, with some clinical trials showing no improvement on lipid profile [11, 17,18,19], insulin resistance [11, 18, 20] or blood pressure [17,18,19,20]. Differences in study designs and populations tested [21] could explain these equivocal results, highlighting the need for further research on blueberry health properties. In addition, inter-individual variability in cardiometabolic responses to nutrition intervention is observed and can be partly attributable to genetic variations [22]. In order to inform future investigations specific to blueberries and polyphenols as they relate to MetS, a better understanding of the physiological mechanism of action for blueberries (polyphenols) on MetS phenotypes is needed.

Accordingly, this 8-week parallel, randomized (1:1) double-blind placebo-controlled clinical trial of adults at risk of MetS had two overarching purposes. First, the study aimed to evaluate the impact of highbush blueberry powder (BBP) on the following MetS parameters: plasma insulin resistance (HOMA-IR) and sensitivity (Matsuda index), plasma lipids, arterial blood pressure, blood glucose, weight, body mass index (BMI) and waist circumference. Second, the study aimed to delineate the mechanism of action of highbush BBP on MetS parameters in overweight/obese adults through analyses of changes in key molecular signalling pathways and metabolic regulatory networks identified using analyses of transcriptomics and metabolomics. It was hypothesized that BBP have a significant effect on cardiometabolic MetS phenotypes and pathways over an 8-week follow-up period.


Study population

Recruitment of study participants took place between September 2017 and November 2018 in the greater Québec City metropolitan area, with follow-up occurring between October 10, 2017 and January 14, 2019. All recruitment and intervention visits took place at the Institute of Nutrition and Functional Foods (INAF), Université Laval. Recruitment ended when the sample size was achieved. A sample size calculation using published data [17] indicated that a total of 28 subjects per group were needed to observe significant changes in fasting plasma triglyceride (TG) and insulin taking into account an anticipated dropout rate of 20%, with a power of 80%. Using group electronic messages sent to university and INAF members as well as newspapers and social media advertisements, a total of 657 individuals contacted the researchers for information on this study. Overall, 110 individuals met email or phone pre-screening criteria and were scheduled for a screening visit at the clinical investigation unit of INAF. Before each visit, participants were asked to follow a 12-h overnight fast and abstain from alcohol for 48 h. They were also asked to avoid vigorous exercise for 24 h prior the visit. Inclusion criteria were as follows: Caucasian men or premenopausal, non-pregnant and non-lactating Caucasian women in good health, having a BMI between 25 and 40 kg/m2 or a waist circumference ≥ 94 cm for men and ≥ 80 cm for women. All participants had to have at least one of the following: fasting plasma TG ≥ 1.35 mmol/L or fasting insulin concentration ≥ 42 pmol/L using our new analytic method and corresponding to a threshold value of 60 pmol/L with the former method that was predictive of a higher risk of cardiovascular disease in the Quebec population [23]. Participants were excluded from the study if they self-reported that they were diagnosed with diabetes, hypercholesterolemia or hypertension or that they were taking medications for these conditions; were allergic or intolerant to blueberries; had a taste aversion to blueberries; were taking medication affecting study parameters; had taken antibiotics, supplements or natural health products on a regular basis over the past 3 months, had undergone surgery in the last 3 months or had planned surgery during the duration of the study; were nicotine users; followed unique dietary patterns (such as a vegan, gluten-free or ketogenic diet); had lost or gained > 5% of their body weight in the last 3 months; or were having more than 2 alcoholic drinks per day on a regular basis. Participants also had to be willing to commit to the study instructions from the run-in period until the end of the intervention. These instructions were to eat a maximum of: 2 portions of berries weekly (including whole fruits, juice, jam and desert but excluding the placebo/blueberry powders); 140 g of food containing cocoa weekly; 1 cup of tea or 4 cups of coffee daily; and 2 alcoholic drinks per week. Red wine and port were prohibited due to their flavonoid content. Participants were also asked to report the use of new medications and to avoid natural health product consumption or changes in lifestyle habits. The inclusion/exclusion criteria were adjusted following the commencement of the trial to enhance recruitment (original criteria: BMI 25.0–40.0 kg/m2 and abdominal obesity; insulin ≥ 60 pmol/L and TG ≥ 1.5 mmol/L); the criteria listed above reflect the revised, most current criteria. There were no other changes to the methods or outcomes after the trial commenced. The inclusion criteria can be found on the clinical trial registry at (NCT03266055). This study was approved by the Ethics Committee of Université Laval. All subjects signed a written informed consent prior their participation to the study.

Study protocol and intervention

This study is a randomized double-blind placebo-controlled clinical intervention. Sex-stratified block randomization was applied electronically on the INAF electronic management platform to ensure a 1:1 treatment ratio for each sex to either blueberry (BBP) or placebo powder. Powders were both provided by the US Highbush Blueberry Council. The BBP consisted of a blend of milled freeze-dried highbush blueberries from two cultivars (Vaccinium virgatum (ashei) and Vaccinium corymbosum), in a 1:1 ratio. The placebo powder was isocaloric and of similar aspect and taste as the BBP, but contained a mixture of dextrose, maltodextrin, fructose, citric acid, malic acid, natural and artificial flavor, xanthan gum, silicon dioxide, and FD & C Red 40 and Blue 2 lakes. Both powders were packaged in identical aluminum packets identified with the letters A or B (for blinding) by the US Highbush Blueberry Council before sending packets to the study centre. Sufficient powder packets were provided at each visit to participants. They had to take a total of 50 g of powder daily divided in two packets of 25 g each, taken at least 8 h apart. Fifty grams of freeze-dried BBP is equivalent to about 350 g (2 and 1/3 cups) of fresh blueberries. Participants were asked to dilute the content of each packet in 300 ml water in a plastic cup provided by the study coordinator or to mix it into foods that were already part of their usual diet such as breakfast cereals, milk, yogurt or smoothies. They were asked not to heat the powder or add it to a hot liquid or food. To assess possible adverse events, participants completed a self-administered side effects questionnaire at weeks 4 and 8. Changes in plasma insulin, glucose, lipids and lipoproteins were primary outcomes. Changes in gene expression, metabolites, and blood pressure were secondary outcomes. All other analyses were exploratory.

When eligibility was confirmed, participants were enrolled in the intervention and invited to the research institute for the baseline visit (week − 2) prior to a 2-week run-in period (Fig. 1). In that period, participants had to commit to the abovementioned study directives without taking the powder. At the start of the intervention (week 0), participants were randomized to consume either the BBP or the placebo powder for 8 weeks and scheduled for visits at the research center at 4-week intervals (weeks 4 and 8). The randomization list was created by VGa. Clinical coordinators (VGa, MK, VGu) were responsible for enrolling participants and assigning them to interventions. Compliance to the study protocol and supplementation was documented by the completion of a journal and by returning all powder packets (both empty and remaining) at each visit.

Fig. 1

Graphical representation of the study protocol. Study design graphical representation from recruitment to the end of the supplementation period. The blue line represents the intervention period. Abbreviation: OGTT, oral glucose tolerance test

Dietary assessment and questionnaires

Dietary habits were assessed on weeks 0, 4 and 8 by a web-based, self-administered, past-month food frequency questionnaire (FFQ). This FFQ was previously validated for French-speaking Canadian adults [24]. Briefly, participants were asked to report their consumption of 136 items grouped in 8 categories including dairy products, fruits, vegetables, meat and alternatives, cereals and grain products, beverages, “other foods” and supplements. Powders’ nutritional values were added to intakes reported on the FFQs completed at week 4 and week 8 with consideration of individual compliance percentage. Women and men with intakes lower than 600 or 800 kcal per day or greater than 3500 or 4200 kcal/day, respectively, were considered as under- or over-reporters and as such, their dietary intakes were not included in the analyses [25]. On the final visit, all participants completed a short questionnaire about the powders’ taste acceptability and consumption.

Anthropometric measures

Weight was measured using a BWB-800 electronic scale (Tanita, Arlington Heights, IL) to the nearest 0.1 kg with patients wearing light indoor clothes and no shoes. Height and waist circumference were measured to the nearest millimeter according to procedures recommended at the Airlie conference [26]: waist was defined as the midpoint between the last floating rib and the top of the iliac crest. The mean of three consecutive measures was used for analyses. Systolic (SBP) and diastolic (DBP) blood pressure were measured while sitting on a chair after a 10-min rest. Again, the mean of three measures repeated at 3-min intervals was used for analyses.

MetS phenotype parameters

Blood samples were drawn at each visit. Fasting plasma lipid profile (total cholesterol (Total-C), HDL cholesterol (HDL-C), LDL cholesterol (LDL-C) and triglyceride (TG)), glucose, insulin and glycated hemoglobin (HbA1c) concentrations were measured at the Centre Hospitalier Universitaire de Québec (CHUQ)-Université Laval. Total-C, HDL-C, TG and glucose were measured by enzyme-based assays. LDL-C was calculated with the Friedewald formula [27], and insulin was assessed by chemiluminescence. HbA1c was measured by ion exchange high-performance liquid chromatography.

Participants underwent oral glucose tolerance tests (OGTT) at the beginning (week 0) and at the end (week 8) of the supplementation period. Values for the fasting blood samples were calculated as the mean of times − 15 and 0 min before participants were asked to drink a 75-g glucose solution. Blood was further drawn at 15, 30, 60, 90 and 120 min after glucose solution intake. The homeostatic model assessment of insulin resistance (HOMA-IR) was computed as fasting glucose and fasting insulin product divided by the normalizing factor 22.5 [28]. Matsuda index [29] was calculated as 10,000 divided by square root of (Fasting glucose × Fasting insulin) × (Mean glucose × Mean insulin), with mean glucose and insulin concentrations being calculated as the mean of fasting, 30-, 60-, 90- and 120-min samples. Units of measurement were converted in order to comply with the units used for these equations. Due to high hemolysis in some blood samples, Matsuda index has been calculated for only a limited number of subjects in each group.

Statistical analysis (MetS phenotype parameters)

SAS version 9.4 (SAS Institute, Cary, NC) was used for statistical analyses, which were by originally assigned groups. Each variable was assessed for normality using skewness and kurtosis. Non-normally distributed variables that have been transformed prior statistical analyses are indicated in footnotes below each table. Statistics are presented as means with standard deviations (SD) or ranges. The MIXED procedure in SAS (with compound symmetry as covariance structure) was used to test for differences in dietary intakes and indices of cardiometabolic health including HOMA-IR and Matsuda index. Treatment, time and treatment × time interaction were considered as fixed effects with repeated measurements. Post hoc comparisons among groups were performed using least square means (LS-means) when a significant treatment × time interaction (p < 0.05) was observed. Analyses were performed with and without adjustments for covariates including age, sex and BMI and/or nutrient intakes, as per the footnote of each table. A per-protocol analysis was used. P values < 0.05 were considered significant.


Blood cell gene expression profiles were examined to assess nutrition-related metabolic changes occurring in response to the BBP. Blood samples were collected in PAXgene blood RNA tubes (Qiagen, Valencia, CA, USA) at week 0 and week 8, and were stored at − 80 °C until the analyses. Total RNA was extracted, and samples were sent to the McGill University and Génome Québec Innovation Centre for sequencing. The quality of RNA samples was evaluated with the 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). RNA samples were converted to cDNA with the Illumina NEB stranded mRNA library preparation kit (Illumina, San Diego, CA, USA; rRNA-depleted stranded (HMR)) for sequence library preparation based on the manufacturer’s protocol. Final libraries were sequenced on an Illumina NovaSeq6000 S4 sequencer using paired-end, 100 bp reads. For bioinformatic analyses, raw reads were trimmed for length (n = 50), quality (phred33 score ≥ 30) and adaptor sequence using Trim Galore (v0.6.5), a wrapper tool around Cutadapt (v1.15) and FastQC (v0.11.9). Trimmed reads were pseudo-aligned to the GRCh38 human reference transcriptome using kallisto (v0.46.2), and transcript abundance was estimated with default parameters and 100 bootstraps, and reported in estimated counts [30]. Data normalization and automatic filtering of estimated counts, as well as differential expression analysis were performed with edgeR v3.28.1 [31]. Given the paired nature of our samples, differential transcript expression across the BBP group from week 0 to week 8 was determined using a generalization of a paired t test implemented in the quasi-likelihood framework of edgeR. Differentially expressed transcripts between post- vs. pre-intervention with BBP were considered at a false-discovery rate (FDR)-corrected P value < 0.05. The functional significance of genes showing at least a 25% difference (1.25-fold change) and an unadjusted P value < 0.05 between pre- and post-supplementation states was explored by pathway enrichment analysis using the clusterProfiler v3.16.0 R package [32]. ClusterProfiler implements statistical methods to analyze and visualize functional profiles of genes/gene clusters and produces adjusted P values using the Benjamini-Hochberg procedure (BH-p) for significantly enriched pathways. The following pathway databases were used for functional enrichment analysis: Gene Ontology Biological Processes (GO-BP), Reactome and the Kyoto Encyclopedia of Genes and Genomes (KEGG). The cnetplot function implemented in the clusterProfiler package was used to visualize significantly enriched pathways.


Plasma samples collected at weeks 0 and 8 were sent to the Analytical Facility for Bioactive Molecules at the Hospital for Sick Children in Toronto, Canada. The quantitative analysis of 630 metabolites from 26 biochemical classes was performed in paired blood samples from 24 participants before and after the BBP supplementation with the MxP® Quant 500 kit for targeted metabolic profiling (Biocrates Life Sciences AG, Innsbruck, Austria). The analysis of metabolites combined flow injection analysis (FIA) with liquid chromatography-based triple quadrupole mass spectrometry (LC-MS/MS) and was performed in an Agilent 1200 series HPLC chromatograph coupled to a SCIEX QTrap 5500 mass spectrometer.

Statistical analysis (metabolomics)

Metabolite data from 24 participants (pre- versus post-intervention with BBP) were processed using the MetaboAnalystR package (v3.0) [33]. First, 113 metabolites with a constant or single value across samples were found and deleted. Non-informative signals were further filtered out based on the interquartile range estimate, and samples were normalized by quantile normalization. Metabolite data were log-transformed and scaled by Pareto scaling (mean-centered and divided by the square root of the standard deviation of each variable). After quantile normalization, one additional metabolite with a constant value was found and deleted. From the original 630 metabolites analyzed, a total of 386 were included in the final statistical analysis. Paired t tests were used to analyze within-subject changes in blood metabolite levels between pre- and post-supplementation states. A volcano plot associated to paired t tests was further used to visualize the most differential metabolite changes between pre- and post-supplementation. A P value < 0.05, along with a count of significant pairs higher than 50% showing at least 25% difference (1.25-fold change), were the criteria used to consider metabolite blood levels to significantly differ between pre- and post-supplementation.

Dimensional reduction

Dimensional reduction was conducted using the partial least squares discriminant analysis (PLS-DA), a supervised algorithm able to reduce the number of metabolites in high-dimensional metabolomics data to produce robust and easy-to-interpret models. This method is able to differentiate the class membership through multivariate regression of a given set of metabolites. In the present study, we used a variation of PLS-DA, the multilevel PLS-DA (mPLS-DA), given its ability to exploit the paired structure of the multivariate data obtained before and after the BBP supplementation in the same group of participants [34, 35]. A sparse mPLS-DA (smPLS-DA) was used to identify the most important metabolites that help discriminate matched study groups [36]. The smPLS-DA algorithm was implemented using the mixOmics R package (v6.12.1) [37]. The variable importance in projection (VIP) coefficients were computed as a weighted sum of squares of the smPLS-DA loadings to depict the relative importance of each metabolite in the classification model. Predictors with large VIP are the most relevant for discriminating class membership [37].


Participants characteristics and adherence

Of the 110 individuals screened for eligibility, 59 participants were included in this randomized double-blind, placebo-controlled intervention trial, with analyses conducted on a final sample of 49 individuals who consumed the BBP (n = 25) or the placebo powder (n = 24) (Fig. 2). Characteristics of participants at week 0 who completed the study protocol are presented in Table 1. Participants were adults between the ages of 22 and 53. Adherence to the study protocol was robust, with an overall powder intake compliance rate of 92.7 ± 7.5% and no significant differences between groups (BBP: 92.4 ± 7.7% and placebo: (92.9 ± 7.3%). Results from the side effects questionnaire indicated that participants experienced some minor side effects with both the BBP and placebo powder (Supplementary Figure 1), but these side effects were generally tolerable with only one participant dropping out of the study as a result of persistent abdominal discomfort.

Fig. 2

CONSORT 2010 Flow Diagram

Table 1 General baseline characteristics of participants

Dietary intake

As presented in Table 2, participants had dietary intakes of approximately 2000 kcal/day. There were no significant between-treatment differences. Time effects were observed for carbohydrates, total sugars and fibre. These results relate to the nutritional content of the BBP and placebo powder, which are comprised of similar quantities of carbohydrates and total sugar with a higher fibre content in the BBP compared to the placebo powder (Supplementary Table 1); thus, these results relate to the compliance to the intervention protocols. There was a significant treatment × time interaction for fibre whereby the BBP group had significantly greater changes (increases) in fibre intake compared with the placebo group from week 0 to weeks 4 and 8.

Table 2 Effects of treatment and/or time on nutritional intake

MetS Phenotypes

Cardiometabolic risk factors at different time points during the intervention are presented in Table 3. A treatment effect was observed only for plasma TG levels. This difference was no longer significant after adjustments for age, sex, BMI and baseline values (data not shown). These results were also no longer significant after adjusting for calories, fibre or total sugar intake, individually, therefore suggesting that the addition of calories, fibre or total sugar from the interventions contributed to the results (data not shown). Treatment × time interactions were all non-significant.

Table 3 Indices of cardiometabolic health over time by intervention group

Changes from week 0 to week 8 in plasma concentrations of cardiometabolic risk factors were also investigated. As presented in Supplementary Table 2, changes were not significantly different between groups either in the unadjusted or the age, sex and BMI-adjusted models.

Data from the OGTT-derived markers of insulin resistance are presented in Table 4. For HOMA-IR, a threshold value of 2.5 is used to define individuals with IR [38]. A time effect was observed for HOMA-IR, with the mean being lower than the defined threshold value for BBP and placebo at week 0 but increased above the threshold after 8 weeks in both groups. For Matsuda index, a value of less than 4.3 predicts IR [38]. Only the placebo group had a mean value below this threshold at 8 weeks. After adjusting for calories, fibre or total sugar intake, individually, the time effect remained significant for Matsuda index (P < 0.04), was no longer significant for HOMA-IR (P > 0.05). No treatment × time interactions were revealed. All results remained unchanged after adjustment for age, sex and BMI (data not shown).

Table 4 Effects of treatment and/or time on insulin resistance (HOMA-IR) and insulin sensitivity (Matsuda index)

Transcriptomics: effect of gene expression

Results from the transcriptomics analyses are presented in Supplementary Table 3 and Figs. 3 and 4. The assessment of differential transcript expression revealed a total of 49 genes differentially expressed at FDR-adjusted P value < 0.05, which were differentially expressed from week 0 to week 8 in the BBP intervention group (Supplementary Table 3). Statistically significant fold changes in differentially expressed genes ranged from − 1.4 (NT5C3A, NM_001002010) to 1.4 (DDX11L5, NR_051986), with 12 genes showing at least a 1.25-fold change (Fig. 3). Individual changes in gene expression for the top differentially up- and downregulated genes are shown in Fig. 4. Pathway enrichment analysis revealed that the 259 most differentially expressed genes were clustered into immune-related pathways, as shown by the top-five significantly enriched GO-BP categories (Fig. 5a), including “defense response to virus” (n = 16 genes, BH-p = 3.0 × 10-5), “response to virus” (n = 18, BH-p = 3.7 × 10-5), “cellular extravasation” (n = 8, BH-p = 3.6 × 10-4), “response to lipopolysaccharide” (n = 15, BH-p = 3.2 × 10-3) and “response to molecule of bacterial origin” (n = 15, BH-p = 4.1 × 10-3). This is consistent with the five significantly enriched pathways found in the Reactome database (Fig. 5b), including “interferon alpha/beta signaling” (n = 8 genes, BH-p = 1.4 × 10-3), “interferon signaling” (n = 10, BH-p = 3.3 × 10-2), “antimicrobial peptides” (n = 7, BH-p = 3.3 × 10-2), “alpha defensins” (n = 3, BH-p = 3.3 × 10-2) and “caspase-mediated cleavage of cytoskeletal proteins” (n = 3, BH-p = 4.8 × 10-2). Finally, an immune-related KEGG pathway was also found to be significantly enriched, the “NOD-like receptor signaling pathway” (n = 10, BH-p = 9.8 × 10-3; network plot not shown).

Fig. 3

Global gene expression change between pre- and post-supplementation states in the blueberry group. MA plot shows the log2 average abundance of transcripts in counts per million mapped reads (log CPM) on the x-axis and the log2-fold change (log FC) on the y-axis. Non-significant genes are represented by grey dots. Over- and under-expressed genes (FC > 1.25) with unadjusted significant differences (paired t test P value < 0.05) are coloured in green and red, respectively. Significant differentially expressed genes from paired t tests (FDR-adjusted P value < 0.05) and showing at least a 1.25 FC are labelled with gene names. The dashed lines represent 1.25 FC

Fig. 4

Top differentially expressed genes between pre- and post-supplementation states in the blueberry group. Box and whisker plots show median, first, and third quartiles, and maximum and minimum values for the 24 sample pairs before (Pre) and after (Post) the blueberry supplementation. The three transcripts which exhibited the most significant (FDR-adjusted P value < 0.05) over- and under-expression derived from paired t tests (Post vs Pre) are shown on the top and bottom rows, respectively. Green and red lines stand for increasing or decreasing gene expression levels between pre- and post-supplementation states within individual paired samples

Fig. 5

Network plots of enriched terms following the blueberry supplementation. The network plot depicts the linkages among differentially expressed gene clusters and functional enriched terms in the Gene Ontology Biological Processes (GO-BP) (a) and Reactome (b) pathway databases. The size of the grey dots is proportional to the number of genes in the enriched pathway (from 3 to 16) and the red-to-green color gradient of gene dots represents the direction of the gene expression fold change following the blueberry supplementation from down- to up-regulation, respectively

Metabolomics: effect of BBP on metabolites

The paired t test analysis demonstrated that following BBP supplementation, a total of 50 metabolites had significantly different blood levels (paired t test unadjusted P value < 0.05), as compared to pre-supplementation levels. This number was reduced to 35 metabolites after applying the statistical significance criteria (> 1.25-fold change and > 50% of significant metabolite counts) (Fig. 6 and Supplementary Table 4). None of the metabolites analyzed showed significant differences following multiple testing correction. Metabolites showing a significant reduction of their blood levels following the BBP supplementation are shown on the top-left corner of the volcano plot, and those showing a significant increase are shown on the top-right corner of the paired volcano plot (Fig. 6). The top-five under- and overabundant metabolites are shown in paired boxplots in Fig. 7. Among the under-abundant metabolites, those showing a more significant decrease are ornithine (Orn), hypoxanthine, diacylglycerol DAG(16:1/18:2), ceramide Cer(d16:1/24:0) and indoxyl sulfate (Ind-SO4). On the other hand, triglycerides TG(16:0/30:2), TG(20:2/34:2), TG(16:1/32:0), TG(14:0/34:3), TG(16:1/32:0) are among the metabolites that have undergone a significant increase following the BBP supplementation.

Fig. 6

Impact of BBP supplementation on blood metabolite levels. Volcano plot of paired comparisons between metabolite blood levels in pre- and post-supplementation groups. On the x-axis, a count of significant sample pairs is shown. On the y-axis, the minus logarithm of paired t test P values is shown. Blood levels of a given metabolite were considered significantly different between pre- and post-supplementation states when the paired t test P value was < 0.05, the change in metabolite blood levels was higher than 25% (> 1.25-fold change), and the count of significant pairs was higher than the 50% of the total count of pairs. Each dot represents a metabolite. Metabolites showing statistically significant changes following the blueberry supplementation are depicted as blue dots on the right (increase) and left (decrease) top corners. Top-ten significantly different metabolites are labelled. Orn ornithine, DG diacylglycerol, Cer ceramide, Ind-SO4 indoxyl sulfate, TG triglyceride, HipAcid hippuric acid

Fig. 7

Top metabolites showing significant changes following BBP supplementation. Box and whisker plots show median, first, and third quartiles and maximum and minimum values for the 24 sample pairs before (Pre) and after (Post) the blueberry supplementation. The five metabolites which exhibited the most significant increases and decreases following the supplementation are shown on the top and bottom rows, respectively. Green and red lines stand for increasing or decreasing metabolite blood levels between pre- and post-supplementation states within individual paired samples. Orn ornithine, DG diacylglycerol, Cer ceramide, Ind-SO4 indoxyl sulfate, TG triglyceride, HipAcid hippuric acid

Results from smPLS-DA strengthened findings from t test analysis, by identifying shared subsets of metabolites. The score plot derived from smPLS-DA shows the complete separation of pre- and post-supplementation groups, without overlap (Fig. 8). Component 1 was primarily responsible for group discrimination, accounting for a 13% of the variance, with component 2 accounting for 2.7%. The ten metabolites associated with the first component and underlying the discrimination between pre- and post-supplementation groups are shown in the loadings panel of the first component (Fig. 8). Most of these metabolites, whose relative contribution to group discrimination, or loading weight, is highlighted in Fig. 8, have been previously identified as differentially abundant between pre- and post-supplementation in the t test analysis (Fig. 7). The top five metabolites exhibiting significant reductions include ornithine, hypoxanthine, DAG(16:1/18:2), Cer(d16:1/24:0) and Ind-SO4. The top five metabolites exhibiting significant increases in response to BBP supplementation include TG(16:0/30:2), TG(20:2/34:2), TG(16:1/32:0), TG(14:0/34:3) and HipAcid.

Fig. 8

Identification of metabolites using Sparse multilevel partial least squares-discriminant analysis (smPLS-DA). A bi-dimensional score plot is shown on the left panel. The score plot reveals the distinct blood metabolomic profile between pre- (red dots) and post-supplementation (green dots) paired participants. The two principal components of the smPLS-DA model along with their corresponding variance in group discrimination are shown on y- and x-axes, respectively. The loading plot representing the top 10 metabolites selected on the first component of the smPLS-DA model is shown on the right. Horizontal bars represent the loading weights of each metabolite. Most important metabolites in group discrimination are ordered according to their loading weights, from bottom to top. Bar colour represents either an increase following supplementation (green bars) or decrease following supplementation (red bars). TG triglyceride, Cer ceramide, DG diacylglycerol, SM sphingolipid, Orn ornithine


This randomized placebo-controlled clinical trial investigated the effects of an 8-week BBP intervention on metabolic markers of cardiometabolic health, transcriptomics and metabolomics in adults at risk of developing MetS. Compared to placebo, BBP intake had no major, significant effect on cardiometabolic health or glycemic response after 8 weeks. However, transcriptomics and metabolomics data indicate significant changes occurring in response to BBP, thus demonstrating that BBP may have potential health-related effects after an intervention longer than 8 weeks and/or with whole blueberry consumption and/or a higher dosage.

Results from studies assessing the impact of highbush blueberries on health have been variable. With a 28-day supplementation of 11 g BBP daily, flow-mediated dilation improved and SBP was significantly reduced in men [39]. Similarly, a reduction in SBP and DBP was observed in a study of comparable size and identical duration and dose as in the present study, but composed mainly of women with the MetS [11]. However, with half of the BBP dose used by Basu et al. [11] and in the present study, no improvement in SBP or DBP following a 6-week BBP supplementation were observed in a sample of men [19]. Through a meta-analysis of studies with similar designs to the present RCT [40], the overall body of evidence suggests that highbush BBP does not have a significant effect on SBP (− 0.28 mmHg [95% CI: − 1.11, 0.56]) and DBP (− 0.50 mmHg [95% CI: − 1.24, 0.24]) [40]. It is however possible that whole blueberries have a significant effect on BP, as a meta-analysis of RCTs reported a favourable effect of 2- to 24-week duration of supplementation with whole berries (including studies on cranberries, bilberries, blueberries, whortleberries, elderberries or raspberries), on SBP (− 2.72 mmHg [95% CI: − 5.32, − 0.12]) [41].

In contrast to findings made in human interventions, numerous animal studies reported benefits of blueberries alone or mixed with other berries on weight management [42,43,44], markers of glucose/insulin resistance or sensitivity [10, 42,43,44,45,46], lipid profile [42, 47, 48], vascular health [49, 50], kidney function [51] or inflammatory markers [10, 42,43,44, 50,51,52,53]. Results from studies assessing other MetS phenotype parameters have been consistent with our findings overall. Using a similar design and intervention protocol, Basu and colleagues reported no difference between the placebo and the blueberry group for serum glucose concentration, lipid profile, HbA1c and HOMA-IR [11]. Similarly, in a study providing blueberry smoothies for a 6-week period, no impact on lipid, fasting glucose and insulin levels or insulin sensitivity were observed [18]. In these two clinical trials, obese participants with MetS were recruited and similar dose of BBP were given as in the present study. Moreover, a meta-analysis reported no effect of berry supplementation on total-C, HDL-C and TG levels, but the berry supplementation did lead to a significant reduction in LDL-C (− 0.21 mmol/L [95% CI: − 0.34, − 0.07]), fasting glucose (− 0.10 mmol/L [95% CI: − 0.17, − 0.03]), HbA1c [− 0.20% [95% CI: − 0.39, − 0.01]) and BMI (− 0.36 kg/m2 [95% CI: − 0.54, − 0.18]) [41]. Results reported herein are comparable to those reported in previous studies and in a meta-analysis including healthy and metabolically deteriorated individuals [40]. It is also important to keep in mind that although at least one cardiometabolic parameter was improved in most human intervention studies with BBP (e.g. blood pressure or inflammation biomarkers), the effect size remained quite small. Moreover, variability in baseline cardiometabolic parameters may contribute to differing results. For example, in the present study baseline BP was within the normal range overall and no significant BP changes were observed, but in a sample of participants with hypertension, BBP may be more likely to affect BP.

Additionally, fasting insulin and glucose levels are limited measures of insulin sensitivity. The gold standard for its assessment is the euglycemic clamp but because it is a strenuous procedure, 120- or 180-min OGTTs are more often used to get information on body glucose disposal efficiency after a test meal or a glucose load [54] from which various indices have been derived to get a better idea of glucose and insulin dynamics [38, 54]. One of them, HOMA-IR, is very common and the simplest to compute [28]. Another, the Matsuda index, is a good approximation of whole-body insulin sensitivity [29, 54]. Many studies have reported on these outcomes in response to blueberry interventions. Similar to the results presented herein, a 12-week tart cherry juice supplementation had no effect on HOMA-IR or insulin concentration [55]. Furthermore, a meta-analysis of pomegranate supplementation reported no change on mean HOMA-IR with a pooled estimated effect of − 0.04 [95% CI: − 0.53, 0.46] [56]. Moreover, a red wine polyphenol supplementation study for 8 weeks did not improve HOMA-IR, fasting glucose and insulin concentrations nor Matsuda index following a mixed-meal test in obese participants [57].

A possible explanation for the absence of effect of the BBP on cardiometabolic risk factors in the present study may be attributable to the extra carbohydrates and sugars provided by the powders, being added to the usual diet of study participants. While the analysis of energy intake did not reveal significant changes between baseline and treatment weeks, participants of both groups consumed significantly more carbohydrates at weeks 4 and 8 in comparison to week 0, of which about two third were sugars. Higher carbohydrate intake has been associated with lower HDL-C and higher TG concentrations [58]. The net increase in the intake of carbohydrates from the powders might have masked the potential lowering effects of BBP antioxidants on cardiometabolic parameters if the participants did not compensate by reducing their consumption of carbohydrates from other foods and beverages. In addition, the heterogeneity of participants whereby some had elevated TG, some had elevated fasting insulin and others had both high TG and high insulin levels may have obscured the results and reduced the probability of detecting significant associations. We cannot rule out the possibility that cardiometabolic parameters are differentially affected depending on an individual’s profile at the entry in the study; this should be further investigated.

The transcriptomics and metabolomics analyses allowed for the identification and understanding of activated metabolic pathways following the BBP dietary intervention. BBP had a significant impact on 35 individual metabolites and the expression of 49 genes, thus providing other important insights for future research exploring health-related outcomes resulting from BBP supplementation. It is possible that BBP had an effect on markers of inflammation, immunity and oxidative stress, but these outcomes were not investigated in the present study. Further exploration of the roles of each of these metabolites and genes can help guide future research endeavours. As a targeted metabolomics profiling, we acknowledge that this study is restricted to a limited number of metabolites and biochemical classes; although the targeted profile was comprehensive, exploring microbiota-derived metabolites other than secondary bile acids, indoles or branched-chain amino acids would have been of interest to the present study. Nevertheless, some interesting findings are worth highlighting. For example, ornithine, produced via the urea cycle, is a non-essential amino acid [59]. The substrate action of ornithine can lead to excessive polyamine synthesis which plays a role in modulating the development of certain types of cancer [59, 60]. This protective role has been supported by previous, in vitro, preclinical and clinical research on the effects of blueberries and other berries [61]. Indeed, future research is needed to explore this further and investigate the metabolic relevance of a reduction in plasma ornithine concentrations following blueberry supplementation. There was also a significant reduction in hypoxanthine following BBP supplementation. Circulating concentrations of hypoxanthine, a purine molecule which is a by-product of adenosine triphosphate catabolism, are elevated following an ischemic event [62]. On the other hand, the increase in plasma hippuric acid, we observed has been previously reported following anthocyanin supplementation [20, 63,64,65]. The significant increases in several TG metabolites following BBP supplementation was, to our knowledge, a novel finding. However, this finding was perhaps not surprising given the well-established relationship between carbohydrate intake and plasma TG [66]. Therefore, this finding of increased TG metabolites further relates to our dietary results demonstrating a significant increase in carbohydrates and total sugars throughout the intervention (due to the nutritional breakdown of the intervention products: Supplementary Table 1). While previous research has demonstrated that blueberries play a role in combatting oxidative stress [11], we did not observe significant changes in gene expression or metabolites related to oxidative stress pathways or 8-iso-PGF2α plasma levels [67].

Pathway analyses can determine if differentially expressed genes are part of predefined physiological networks more than what would be expected by chance alone. This allows for the generation of mechanistic hypotheses and identification of putative mechanisms [68]. It was interesting to find that differentially regulated genes were clustered into immune-related pathways, thus suggesting that the BBP supplementation could have anti-inflammatory effects. This finding is consistent with a formerly conducted transcriptomic analysis following a 4-week, 1L/day blueberry-apple juice dietary intervention, whereby the researchers also found gene expression changes in immune-response pathways alongside signalling pathways for apoptosis, cell adhesion and lipid metabolism [69]. This further relates to previous research indicating that blueberries can have immunomodulatory effects and reduce oxidative stress in adults with MetS [70], and that a blueberry green tea polyphenol soy complex could have a potential protective role against viral infections in athletes [71]. Moreover, a randomized controlled trial of 38 g/day BBP for 6 weeks resulted in a significant increase in natural killer cells, a type of peripheral lymphocyte playing a key role in the immune response [72]. It is important to note that our analysis was performed in whole blood samples, with gene expression analysis being consequently carried out on circulating immune cells, which are the most relevant cells to the immune system [73] and thus differential gene expression related to immunity may be expected. However, previous research has demonstrated that white blood cells, mainly peripheral blood mononuclear cells, adequately reflect the expression of the majority of genes in a metabolically relevant tissue such as the skeletal muscle tissues [74]. The identified changes in key molecular signalling pathways and metabolic regulatory networks can ultimately be used to identify relevant biological and metabolic pathways that can become targets for therapy and may assist researchers in developing new hypotheses related to the health benefits of blueberries. However, although it is also possible that a nutritional intervention may alter the composition of different cell types in blood, which then causes the changes in gene expression, the BBP intervention did not seem to induce any significant change in immune cell count (data not shown).

Taken together, while BBP did not appear to impact specific cardiometabolic risk factors related to MetS within the given time frame of the study, the significant changes in gene expression and metabolites indicate the potential for BBP to impact various health outcomes, particularly those related to immunity. MetS-related health outcomes may occur beyond an 8-week intake period or from an intervention with whole blueberries rather than BBP or a higher dosage. For example, a long-term BBP supplementation study in individuals with MetS investigated the benefits of taking 13 g or 22 g BBP on a daily basis (representing respectively a quarter and a half of the dose given in the present study), for a period of 6 months found significant improvements of endothelial function and arterial stiffness in subjects consuming the highest dose of BBP [20].

The present study had some limitations. First, it included premenopausal women, but menstrual cycle hormonal fluctuation can influence cardiometabolic biomarkers especially those associated to the lipid profile [75] and endothelial function [76]. Since the only constraint regarding powder intake was not to heat it, it is possible that consuming it with food such as milk could have had an effect on its antioxidant properties due to the matrix effect [77]. Furthermore, milled freeze-dried highbush BBP was used as a surrogate of whole fresh blueberries for practical reasons. Although this method has been reported to preserve blueberry antioxidants [78, 79], we cannot discard the possibility that fresh blueberries may exert more potent cardiometabolic benefits than those freeze-dried and milled into powder. Nonetheless, the strengths of the present study include its randomized, placebo-controlled, double-blind design that limits the influence of confounding factors, which may bias estimates of treatment effects. Compliance to treatment was high and study directives were well followed by participants. However, given the overall significant increase in carbohydrates and total sugars over time, it is not possible to differentiate between an independent BBP effect or an overall macronutrient effect on metabolic parameters within the present study. Powders were of similar aspect and taste and, unless comparing BBP and placebo powder side to side, participants were probably unable to discover to which treatment they were randomized. Also, it is important to highlight that blueberry supplementation studies in healthy but overweight individuals at risk of developing MetS are less frequent thus highlighting the novelty of this work. Most studies have been conducted on subjects older than those in the present study and/or having a MetS diagnosis. Moreover, in addition to reporting changes for traditional biomarkers of cardiometabolic health, this trial also includes analyses on glucose and insulin resistance/sensitivity indices, which are better indicators than fasting values alone. The literature was especially lacking human interventions reporting these measurements. Lastly, to our knowledge, this was the first study to assess changes in gene expression resulting from BBP supplementation. Finally, this study is primarily generalizable to men and women at risk of MetS.


In conclusion, no significant effects on plasma cardiometabolic risk factors were observed herein following an 8-week intervention with BBP. However, given the significant changes in gene expression and several metabolites following BBP supplementation, it is plausible that a longer follow-up and/or higher dosage could have resulted in significant improvements in plasma MetS parameters. Future BBP interventional clinical trials should seek to explore health outcomes related to the changes in gene expression and metabolites reported in the current trial, while providing an intervention > 8 weeks, intervening with whole blueberries, and/or providing a higher dosage of BBP.

Availability of data and materials

Please contact authors for data requests.



Blueberry powder


Body mass index


Diastolic blood pressure


Food frequency questionnaire


Insulin resistance


Metabolic syndrome


Reactive oxygen species


Systolic blood pressure




  1. 1.

    Aguilar M, Bhuket T, Torres S, Liu B, Wong RJ. Prevalence of the metabolic syndrome in the United States, 2003-2012. JAMA. 2015;313(19):1973–4.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Riediger ND, Clara I. Prevalence of metabolic syndrome in the Canadian adult population. Can Med Assoc J. 2011;183(15):E1127–34.

    Article  Google Scholar 

  3. 3.

    O’Neill S, O’Driscoll L. Metabolic syndrome: a closer look at the growing epidemic and its associated pathologies: metabolic syndrome. Obes Rev. 2015;16(1):1–12.

    Article  PubMed  Google Scholar 

  4. 4.

    Rani V, Deep G, Singh RK, Palle K, Yadav UCS. Oxidative stress and metabolic disorders: pathogenesis and therapeutic strategies. Life Sci. 2016;148:183–93.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Kehrer JP, Klotz L-O. Free radicals and related reactive species as mediators of tissue injury and disease: implications for health. Crit Rev Toxicol. 2015;45(9):765–98.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Kasote DM, Katyare SS, Hegde MV, Bae H. Significance of antioxidant potential of plants and its relevance to therapeutic applications. Int J Biol Sci. 2015;11(8):982–91.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Wu X, Beecher GR, Holden JM, Haytowitz DB, Gebhardt SE, Prior RL. Concentrations of anthocyanins in common foods in the United States and estimation of normal consumption. J Agric Food Chem. 2006;54(11):4069–75.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Martineau LC, Couture A, Spoor D, Benhaddou-Andaloussi A, Harris C, Meddah B, et al. Anti-diabetic properties of the Canadian lowbush blueberry Vaccinium angustifolium Ait. Phytomedicine. 2006;13(9-10):612–23.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Wolfe KL, Kang X, He X, Dong M, Zhang Q, Liu RH. Cellular antioxidant activity of common fruits. J Agric Food Chem. 2008;56(18):8418–26.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    DeFuria J, Bennett G, Strissel KJ, Perfield JW II, Milbury PE, Greenberg AS, et al. Dietary blueberry attenuates whole-body insulin resistance in high fat-fed mice by reducing adipocyte death and its inflammatory sequelae. J Nutr. 2009;139(8):1510–6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Basu A, Du M, Leyva MJ, Sanchez K, Betts NM, Wu M, et al. Blueberries decrease cardiovascular risk factors in obese men and women with metabolic syndrome. J Nutr. 2010;140(9):1582–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Johnson SA, Figueroa A, Navaei N, Wong A, Kalfon R, Ormsbee LT, et al. Daily blueberry consumption improves blood pressure and arterial stiffness in postmenopausal women with pre- and stage 1-hypertension: a randomized, double-blind, placebo-controlled clinical trial. J Acad Nutr Diet. 2015;115(3):369–77.

    Article  PubMed  Google Scholar 

  13. 13.

    Bowtell JL, Aboo-Bakkar Z, Conway ME, Adlam A-LR, Fulford J. Enhanced task-related brain activation and resting perfusion in healthy older adults after chronic blueberry supplementation. Appl Physiol Nutr Metab. 2017;42(7):773–9.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Miller MG, Hamilton DA, Joseph JA, Shukitt-Hale B. Dietary blueberry improves cognition among older adults in a randomized, double-blind, placebo-controlled trial. Eur J Nutr. 2018;57(3):1169–80.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

    Cassidy A, Mukamal KJ, Liu L, Franz M, Eliassen AH, Rimm EB. High Anthocyanin intake is associated with a reduced risk of myocardial infarction in young and middle-aged women. Circulation. 2013;127(2):188–96.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Muraki I, Imamura F, Manson JE, Hu FB, Willett WC, van Dam RM, et al. Fruit consumption and risk of type 2 diabetes: results from three prospective longitudinal cohort studies. BMJ. 2013;347(aug28 1):f5001.

    Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Stull AJ, Cash KC, Johnson WD, Champagne CM, Cefalu WT. Bioactives in blueberries improve insulin sensitivity in obese, insulin-resistant men and women. J Nutr. 2010;140(10):1764–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Stull AJ, Cash KC, Champagne CM, Gupta AK, Boston R, Beyl RA, et al. Blueberries improve endothelial function, but not blood pressure, in adults with metabolic syndrome: a randomized, double-blind, placebo-controlled clinical trial. Nutrients. 2015;7(6):4107–23.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Riso P, Klimis-Zacas D, Del Bo’ C, Martini D, Campolo J, Vendrame S, et al. Effect of a wild blueberry (Vaccinium angustifolium) drink intervention on markers of oxidative stress, inflammation and endothelial function in humans with cardiovascular risk factors. Eur J Nutr. 2013;52(3):949–61.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Curtis PJ, van der Velpen V, Berends L, Jennings A, Feelisch M, Umpleby AM, et al. Blueberries improve biomarkers of cardiometabolic function in participants with metabolic syndrome—results from a 6-month, double-blind, randomized controlled trial. Am J Clin Nutr. 2019;109(6):1535–45.

    Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Cutler BR, Petersen C, Anandh Babu PV. Mechanistic insights into the vascular effects of blueberries: evidence from recent studies. Mol Nutr Food Res. 2016;61:1600271.

    Article  Google Scholar 

  22. 22.

    Peña-Romero AC, Navas-Carrillo D, Marín F, Orenes-Piñero E. The future of nutrition: nutrigenomics and nutrigenetics in obesity and cardiovascular diseases. Crit Rev Food Sci Nutr. 2018;58(17):3030–41. Taylor & Francis.

    Article  PubMed  Google Scholar 

  23. 23.

    Scarsella C, Alméras N, Mauriège P, Blanchet C, Sauvé L, Dewailly E, et al. Prevalence of metabolic alterations predictive of cardiovascular disease risk in the Québec population. Can J Cardiol. 2003;19(1):51–7.

    PubMed  Google Scholar 

  24. 24.

    Labonté M-È, Cyr A, Baril-Gravel L, Royer M-M, Lamarche B. Validity and reproducibility of a web-based, self-administered food frequency questionnaire. Eur J Clin Nutr. 2012;66(2):166–73.

    Article  PubMed  Google Scholar 

  25. 25.

    Willett W. Nutritional epidemiology. 3rd ed. New York: Oxford University Press; 2013.

    Google Scholar 

  26. 26.

    Loham T, Roche A, Martorel R. Standardization of anthropometric measurements. The Airlie (VA) Consensus Conference. 1988. p. 39–80.

  27. 27.

    Friedewald WT, Levy RI, Fredrickson DS. Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem. 1972;18(6):499–502.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985;28(7):412–9.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Matsuda M, DeFronzo RA. Insulin sensitivity indices obtained from oral glucose tolerance testing: comparison with the euglycemic insulin clamp. Diabetes Care. 1999;22(9):1462–70.

    CAS  Article  Google Scholar 

  30. 30.

    Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol. 2016;34(5):525–7.

    CAS  Article  Google Scholar 

  31. 31.

    McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012;40(10):4288–97.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Yu G, Wang LG, Han Y, He QY. ClusterProfiler: an R package for comparing biological themes among gene clusters. OMICS A J Integr Biol. Mary Ann Liebert, Inc. 2012;16:284–7.

    CAS  Article  Google Scholar 

  33. 33.

    Chong J, Xia J. MetaboAnalystR: an R package for flexible and reproducible analysis of metabolomics data. Bioinformatics. 2018;34:4313–4.

    CAS  Article  Google Scholar 

  34. 34.

    Cao MD, Giskeødegård GF, Bathen TF, Sitter B, Bofin A, Lønning PE, et al. Prognostic value of metabolic response in breast cancer patients receiving neoadjuvant chemotherapy; 2012.

    Book  Google Scholar 

  35. 35.

    Westerhuis JA, van Velzen EJJ, Hoefsloot HCJ, Smilde AK. Multivariate paired data analysis: multilevel PLSDA versus OPLSDA. Metabolomics. 2010;6(1):119–28.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Lê Cao KA, Boitard S, Besse P. Sparse PLS discriminant analysis: biologically relevant feature selection and graphical displays for multiclass problems. BMC Bioinformatics. 2011;12:253.

    Article  Google Scholar 

  37. 37.

    Rohart F, Gautier B, Singh A, Lê Cao KA. mixOmics: an R package for ‘omics feature selection and multiple data integration. PLoS Comput Biol. 2017;13:e1005752 Public Library of Science, Schneidman D, editor.

    Article  Google Scholar 

  38. 38.

    Gutch M, Kumar S, Razi S, Gupta K, Gupta A. Assessment of insulin sensitivity/resistance. Indian J Endocrinol Metab. 2015;19(1):160–4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Rodriguez-Mateos A, Istas G, Boschek L, Feliciano RP, Mills CE, Boby C, et al. Circulating anthocyanin metabolites mediate vascular benefits of blueberries: insights from randomized controlled trials, metabolomics, and nutrigenomics. J Gerontol A. 2019;74(7):967–76.

    CAS  Article  Google Scholar 

  40. 40.

    Zhu Y, Sun J, Lu W, Wang X, Wang X, Han Z, et al. Effects of blueberry supplementation on blood pressure: a systematic review and meta-analysis of randomized clinical trials. J Hum Hypertens. 2017;31(3):165–71.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Huang H, Chen G, Liao D, Zhu Y, Xue X. Effects of berries consumption on cardiovascular risk factors: a meta-analysis with trial sequential analysis of randomized controlled trials. Sci Rep. 2016;6.

  42. 42.

    Lee I-C, Kim DY, Choi BY. Antioxidative activity of blueberry leaf extract prevents high-fat diet-induced obesity in C57BL/6 mice. J Cancer Prev. 2014;19(3):209–15.

    Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Aranaz P, Romo-Hualde A, Zabala M, Navarro-Herrera D, Ruiz de Galarreta M, Gil AG, et al. Freeze-dried strawberry and blueberry attenuates diet-induced obesity and insulin resistance in rats by inhibiting adipogenesis and lipogenesis. Food Funct. 2017;8(11):3999–4013.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Wu T, Gao Y, Guo X, Zhang M, Gong L. Blackberry and blueberry anthocyanin supplementation counteract high-fat-diet-induced obesity by alleviating oxidative stress and inflammation and accelerating energy expenditure. Oxidative Med Cell Longev. 2018;2018:1–9.

    CAS  Google Scholar 

  45. 45.

    Elks CM, Terrebonne JD, Ingram DK, Stephens JM. Blueberries improve glucose tolerance without altering body composition in obese postmenopausal mice: blueberries and postmenopausal obesity. Obesity. 2015;23(3):573–80.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Liu W, Mao Y, Schoenborn J, Wang Z, Tang G, Tang X. Whole blueberry protects pancreatic beta-cells in diet-induced obese mouse. Nutr Metab. 2019;16(1):34.

    CAS  Article  Google Scholar 

  47. 47.

    Vendrame S, Daugherty A, Kristo AS, Klimis-Zacas D. Wild blueberry (Vaccinium angustifolium )-enriched diet improves dyslipidaemia and modulates the expression of genes related to lipid metabolism in obese Zucker rats. Br J Nutr. 2014;111(2):194–200.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Jiao X, Wang Y, Lin Y, Lang Y, Li E, Zhang X, et al. Blueberry polyphenols extract as a potential prebiotic with anti-obesity effects on C57BL/6 J mice by modulating the gut microbiota. J Nutr Biochem. 2019;64:88–100.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Kalea AZ, Clark K, Schuschke DA, Kristo AS, Klimis-Zacas DJ. Dietary enrichment with wild blueberries (Vaccinium angustifolium) affects the vascular reactivity in the aorta of young spontaneously hypertensive rats. J Nutr Biochem. 2010;21(1):14–22.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Mykkänen OT, Huotari A, Herzig K-H, Dunlop TW, Mykkänen H, Kirjavainen PV. Wild blueberries (Vaccinium myrtillus) alleviate inflammation and hypertension associated with developing obesity in mice fed with a high-fat diet. PLoS ONE. 2014;9:e114790.

    Article  Google Scholar 

  51. 51.

    Nair AR, Elks CM, Vila J, Del Piero F, Paulsen DB, Francis J. A blueberry-enriched diet improves renal function and reduces oxidative stress in metabolic syndrome animals: potential mechanism of TLR4-MAPK signaling pathway. PLoS ONE. 2014;9:e111976.

    Article  Google Scholar 

  52. 52.

    Vendrame S, Daugherty A, Kristo AS, Riso P, Klimis-Zacas D. Wild blueberry (Vaccinium angustifolium) consumption improves inflammatory status in the obese Zucker rat model of the metabolic syndrome. J Nutr Biochem. 2013;24(8):1508–12.

    CAS  Article  PubMed  Google Scholar 

  53. 53.

    Vendrame S, Tsakiroglou P, Kristo AS, Schuschke DA, Klimis-Zacas D. Wild blueberry consumption attenuates local inflammation in the perivascular adipose tissue of obese Zucker rats. Appl Physiol Nutr Metab. 2016;41(10):1045–51.

    CAS  Article  PubMed  Google Scholar 

  54. 54.

    Muniyappa R, Lee S, Chen H, Quon MJ. Current approaches for assessing insulin sensitivity and resistance in vivo: advantages, limitations, and appropriate usage. Am J Physiol-Endocrinol Metab. 2008;294(1):E15–26.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Chai SC, Davis K, Wright RS, Kuczmarski MF, Zhang Z. Impact of tart cherry juice on systolic blood pressure and low-density lipoprotein cholesterol in older adults: a randomized controlled trial. Food Funct. 2018;9(6):3185–94.

    CAS  Article  PubMed  Google Scholar 

  56. 56.

    Huang H, Liao D, Chen G, Chen H, Zhu Y. Lack of efficacy of pomegranate supplementation for glucose management, insulin levels and sensitivity: evidence from a systematic review and meta-analysis. Nutr J. 2017;16(1):67.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Woerdeman J, Del Rio D, Calani L, Eringa EC, Smulders YM, Serné EH. Red wine polyphenols do not improve obesity-associated insulin resistance: a randomized controlled trial. Diabetes Obes Metab. 2018;20(1):206–10.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Ma Y, Li Y, Chiriboga DE, Olendzki BC, Hebert JR, Li W, et al. Association between carbohydrate intake and serum lipids. J Am Coll Nutr. 2006;25(2):155–63.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Sivashanmugam M, Jaidev J, Umashankar V, Sulochana KN. Ornithine and its role in metabolic diseases: an appraisal. Biomed Pharmacother. 2017;86:185–94.

    CAS  Article  Google Scholar 

  60. 60.

    Heller JS, Chen KY, Kyriakidis DA, Fong WF, Canellakis ES. The modulation of the induction of ornithine decarboxylase by spermine, spermidine and diamines. J Cell Physiol. 1978;96(2):225–34.

    CAS  Article  PubMed  Google Scholar 

  61. 61.

    Neto CC. Cranberry and blueberry: evidence for protective effects against cancer and vascular diseases. Mol Nutr Food Res. 2007;51(6):652–64.

    CAS  Article  PubMed  Google Scholar 

  62. 62.

    Farthing DE, Farthing CA, Xi L. Inosine and hypoxanthine as novel biomarkers for cardiac ischemia: from bench to point-of-care. Exp Biol Med. 2015;240(6):821–31.

    CAS  Article  Google Scholar 

  63. 63.

    Rodriguez-Mateos A, Del Pino-García R, George TW, Vidal-Diez A, Heiss C, Spencer JPE. Impact of processing on the bioavailability and vascular effects of blueberry (poly)phenols. Mol Nutr Food Res. 2014;58(10):1952–61.

    CAS  Article  PubMed  Google Scholar 

  64. 64.

    de Ferrars RM, Cassidy A, Curtis P, Kay CD. Phenolic metabolites of anthocyanins following a dietary intervention study in post-menopausal women. Mol Nutr Food Res. 2014;58(3):490–502.

    CAS  Article  PubMed  Google Scholar 

  65. 65.

    Nieman DC, Kay CD, Rathore AS, Grace MH, Strauch RC, Stephan EH, et al. Increased plasma levels of gut-derived phenolics linked to walking and running following two weeks of flavonoid supplementation. Nutrients. 2018;10:1718 Multidisciplinary Digital Publishing Institute.

    Article  Google Scholar 

  66. 66.

    Parks EJ, Hellerstein MK. Carbohydrate-induced hypertriacylglycerolemia: historical perspective and review of biological mechanisms. Am J Clin Nutr. 2000;71(2):412–33.

    CAS  Article  PubMed  Google Scholar 

  67. 67.

    Andrianjafimasy M, Zerimech F, Akiki Z, Huyvaert H, Le Moual N, Siroux V, et al. Oxidative stress biomarkers and asthma characteristics in adults of the EGEA study. Eur Respir J. 2017;50(6):1701193.

    CAS  Article  PubMed  Google Scholar 

  68. 68.

    Rudkowska I, Paradis A-M, Thifault E, Julien P, Tchernof A, Couture P, et al. Transcriptomic and metabolomic signatures of an n-3 polyunsaturated fatty acids supplementation in a normolipidemic/normocholesterolemic Caucasian population. J Nutr Biochem. 2013;24(1):54–61.

    CAS  Article  PubMed  Google Scholar 

  69. 69.

    van Breda SGJ, Wilms LC, Gaj S, Jennen DGJ, Briedé JJ, Helsper JP, et al. Can transcriptomics provide insight into the chemopreventive mechanisms of complex mixtures of phytochemicals in humans? Antioxid Redox Signal. 2014;20(14):2107–13.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  70. 70.

    Nair AR, Mariappan N, Stull AJ, Francis J. Blueberry supplementation attenuates oxidative stress within monocytes and modulates immune cell levels in adults with metabolic syndrome: a randomized, double-blind, placebo-controlled trial. Food Funct. 2017;8(11):4118–28.

    CAS  Article  PubMed  Google Scholar 

  71. 71.

    Ahmed M, Henson DA, Sanderson MC, Nieman DC, Gillitt ND, Lila MA. The protective effects of a polyphenol-enriched protein powder on exercise-induced susceptibility to virus infection. Phytother Res. 2014;28(12):1829–36.

    CAS  Article  PubMed  Google Scholar 

  72. 72.

    McAnulty LS, Collier SR, Landram MJ, Whittaker DS, Isaacs SE, Klemka JM, et al. Six weeks daily ingestion of whole blueberry powder increases natural killer cell counts and reduces arterial stiffness in sedentary males and females. Nutr Res. 2014;34(7):577–84.

    CAS  Article  PubMed  Google Scholar 

  73. 73.

    Charles A, Janeway J, Travers P, Walport M, Shlomchik MJ. The components of the immune system. Immunobiology: the immune system in health and disease 5th edition [Internet]. New York: Garland Science; 2001. [cited 2021 Jan 28]; Available from:

    Google Scholar 

  74. 74.

    Rudkowska I, Raymond C, Ponton A, Jacques H, Lavigne C, Holub BJ, et al. Validation of the use of peripheral blood mononuclear cells as surrogate model for skeletal muscle tissue in nutrigenomic studies. OMICS. 2011;15(1-2):1–7.

    CAS  Article  PubMed  Google Scholar 

  75. 75.

    Schisterman EF, Mumford SL, Sjaarda LA. Failure to consider the menstrual cycle phase may cause misinterpretation of clinical and research findings of cardiometabolic biomarkers in premenopausal women. Epidemiol Rev. 2014;36(1):71–82.

    Article  PubMed  Google Scholar 

  76. 76.

    Kawano H, Motoyama T, Kugiyama K, Hirashima O, Ohgushi M, Yoshimura M, et al. Menstrual cyclic variation of endothelium-dependent vasodilation of the brachial artery: possible role of estrogen and nitric oxide. Proc Assoc Am Physicians. 1996;108(6):473–80.

    CAS  PubMed  Google Scholar 

  77. 77.

    Serafini M, Testa MF, Villaño D, Pecorari M, van Wieren K, Azzini E, et al. Antioxidant activity of blueberry fruit is impaired by association with milk. Free Radic Biol Med. 2009;46(6):769–74.

    CAS  Article  PubMed  Google Scholar 

  78. 78.

    Michalska A, Łysiak G. Bioactive compounds of blueberries: post-harvest factors influencing the nutritional value of products. Int J Mol Sci. 2015;16(8):18642–63.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  79. 79.

    Nemzer B, Vargas L, Xia X, Sintara M, Feng H. Phytochemical and physical properties of blueberries, tart cherries, strawberries, and cranberries as affected by different drying methods. Food Chem. 2018;262:242–50.

    CAS  Article  PubMed  Google Scholar 

Download references


The authors gratefully acknowledge participants’ contribution to our study. We also recognize the work of the clinical investigation unit staff and of Maximilien Franck for their involvement in the study coordination, recruitment of the subjects and data collection.


This work was supported by the US Highbush Blueberry Council (USHBC). The funders were not involved in the study design, data analysis, or interpretation of results. Marie-Claude Vohl is a Canada Research Chair in Genomics Applied to Nutrition and Metabolic Health. André Marette is holding a Pfizer research Chair in the pathogenesis of insulin resistance and cardiovascular diseases. Michèle Rousseau received a studentship from the INITIA Foundation. Justine Horne was supported through postdoctoral fellowships from the Canadian Institutes of Health Research, Centre Nutrition, santé et société (NUTRISS) and the Institut sur la nutrition et les aliments fonctionnels (INAF).

Author information




Conceptualization and design of the experiments: VGa, CC, PC, DR, AM and MCV; clinical supervision: PC, VGa, VGu and MK; data analysis: MR, FG and JTM; formal analysis and writing: MR, JH and JTM; supervision: MCV. Review and editing: all authors. MR, JH, JTM, FG, CC, GP, PC, VGa, VGu, MK, DR, AM and MCV critically revised the manuscript. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Marie-Claude Vohl.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Université Laval (2017-135 R-3) and registered at (NCT03266055). All subjects signed a written informed consent prior their participation to the study.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table 1

. Nutrition Facts: Blueberry Powder and Placebo Powder

Additional file 2: Supplementary Figure 1

. Proportion of participants who experienced side effects following powder consumption

Additional file 3: Supplementary Table 2

. Effect of treatment on changes in cardiometabolic variables of interest

Additional file 4: Supplementary Table 3

. Transcriptomics results for significant (p<0.05) changes in gene expression from week 0 to week 8 in BBP group

Additional file 5: Supplementary Table 4

. Metabolomics results for significant (p<0.05) changes in metabolites from week 0 to week 8 in BBP group

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Rousseau, M., Horne, J., Guénard, F. et al. An 8-week freeze-dried blueberry supplement impacts immune-related pathways: a randomized, double-blind placebo-controlled trial. Genes Nutr 16, 7 (2021).

Download citation


  • Metabolic syndrome
  • Blueberry
  • Overweight/obesity
  • Nutrition
  • Immunity
  • Transcriptomics
  • Metabolomics
  • Gene expression