1. Introduction
Northern wild rice (
Zizania palustris, WLD) is an annual plant naturally growing in shallow water in Great Lake and Canadian Prairie regions. WLD was a traditional food of Indigenous people in North America and known as Indian rice or manoomin [
1]. WLD is not a type of rice, but it has been considered as a kind of healthy whole grain [
2]. WLD grain contains significantly more contents of proteins, fiber, vitamins and antioxidants compared to refined white rice (WHR) [
3].
Previous studies demonstrated that WLD reduced hypercholesterolemia, hypertriglyceridemia and the extents of atherosclerotic lesions in low-density lipoprotein (LDL) receptor-knockout mice compared to WHR [
4,
5,
6]. WHR is the most commonly consumed rice product in worldwide regions. The refining process removed most nutrients from the out layer of brown rice. Intake of WHR is associated with increased risks for type 2 diabetes (T2D) [
7,
8]. T2D is characterized by hyperglycemia, insulin resistance, obesity and chronic low-grade inflammation [
9]. Our recent study demonstrated that the replacement of WHR with WLD reduced hyperglycemia, hyperlipidemia, insulin resistance and inflammatory markers in high fat (HF) diet-induced insulin resistant mice [
10]. The underlying mechanism for the metabolic and anti-inflammatory effects of WLD remains to be determined.
Gastrointestinal tract or gut is the major organ in the body to absorb nutrients from foods and plays critical roles in metabolism and inflammation [
11]. In addition, gut is the largest reservoir for microbes in the body. Most of gut bacteria are not harmful to the body, and many of them are even beneficial to health since they help to generate healthy nutrients from foods, modulate metabolism, and inhibit inflammation and the overgrowth of harmful microbes in the gut. The diversities and composition of bacteria in gut was affected by HF diet [
12]. Gut bacteria are essential for HF diet-induced obesity [
13]. Multiple groups of bacteria promote the digest of insoluble fiber from foods and generate short chain fatty acids (SCFAs) in gut. SCFAs are involved in the modulation of metabolism and inflammation in the body [
14].
The present study examined the metabolites in WLD versus WHR, and the effects of the supplementation of two dosages of WLD to HF diet on metabolism, vascular inflammation, pro-inflammatory cytokines, gut microbiota and fecal SCFAs in mice compared to WHR-supplemented HF diet. Correlations between fecal microbiota, SCFAs and circulatory metabolic or inflammatory cytokines in the mice were further investigated.
2. Materials and Methods
2.1. Dietary Components
WLD was obtained from Floating Leaf Wild Rice Inc. (Springfield, MB). Asian white rice was obtained from local grocery store. Carbohydrate-free HF diet powder (cat#: D12492px11) was received from Research Diets (New Brunswick, NJ) which contains 47% of fat (91% from lard) and 35% of protein (weight/weight or w/w). WLD and WHR were grinded and passed through a 0.5 mm mesh. After the supplementation of 26 g% weight of rice powder, experimental HF diet was pelleted, which contained 35 g% fat accounting for 60% of the total calories, 26 g% protein for 20% of calories and 26 g% carbohydrate for the remaining 20% of calories in the diet.
2.2. Animals
C57 BL/6J mice (male, 6 weeks of age) were obtained from the Jackson Laboratory (Bar Harbor, ME). Mice were held in stainless steel cages in air-conditioned room and received regular rodent chow for one week for stabilization. The mice were randomly divided into 3 groups (n=5/group) after the stabilization. WHR diet group received HF diet supplemented with 26 g% of WHR (w/w). WTWD diet group was fed with HF diet supplemented with 13 g% of WHR and 13 g% of WLD. WLD diet group received HF diet supplemented with 26 g% of WLD. The dietary intervention lasted for 12 weeks.
2.3. Animal Monitoring and Sample Collection
Body weights and food intake were recorded at the onset and before the end of dietary intervention. Blood was collected from the saphenous vein to measure blood biochemical variables. At the end of the intervention, mice were euthanized via isoflurane (5%, v/v) inhalation. Abdominal aortae were surgically harvested and submerged in ice-cold Hank’s balance salt solution (HBSS) in culture dishes for
ex vivo monocyte adhesion assay as previously described [
15]. The protocol of the animal experiment was approved by the Animal Management and Protocol Committee at the University of Manitoba.
2.4. Measurements of Circulatory Glucose, Triglycerides and Cholesterol
The levels of plasma glucose, total cholesterol and triglycerides of mice after an overnight fasting were analysed using Sekisui Diagnostics SL reagents (Charlottetown, PE, Canada).
2.5. Measurements of Plasma Insulin and Pro-Inflammatory Cytokines
The levels of insulin, tumor necrosis factor-α (TNFα), plasminogen activator inhibitor-1 (PAI-1) and monocyte chemotactic protein-1 (MCP-1) antigens in plasma were measured using enzyme-linked immunosorbent assay (ELISA) kits from EMD Millipore (Billerica, MA for insulin), BD Bioscience (San Diego, CA for TNFα), Thermo Scientific (Ottawa, ON for MCP-1), and Oxford Biomedical Research (Oxford, MI for PAI-1). Homeostatic model assessment-insulin resistance (HOMA-IR) was calculated from plasma glucose and insulin in simultaneously collected blood samples using a mouse-specific formula [
16].
2.6. Monocyte Adhesion Assay
Mouse abdominal aortae were isolated at the end of dietary regimen and submerged in ice-cold HBSS in 100 mm culture dishes. Fluorescently labeled WEHI-274.1 mouse monocytes (1x10
5) were added to each dish and incubated in Roswell Park Memorial Institute 1640 medium at room temperature with freshly isolated aorta for 30 min on a rotating mixer. Non-adhered monocytes were carefully removed via washes using ice-cold HBSS twice. Monocytes adhered on the surface of aorta were counted under microscopy (10X magnification). The averages of adhered monocytes on the intima of an aorta in five fields were applied for data analysis [
15].
2.7. Fecal Bacteria DNA Extraction and 16S rRNA Gene Sequencing
Feces were collected from cages of mice individually hosted for overnight. DNA was extracted from feces using Power Fecal DNA Isolation Kit (QIAGEN, Germantown, MD) and quantified using a NanoDrop 2000 spectrophotometer (Thermo Scientific). Fecal DNA was amplified using primers targeting the V4-V5 region of bacterial DNA sequence [515F (5′-GTGYCAGCMGCCGCGGTAA) and 926R (5′-CCGYCAATTYMTTTRAGTTT)]. The DNA amplicons were normalized using Charm Biotech 96-well normalization kit for 16S rRNA gene sequencing on an Illumina MiSeq sequencer in Integrated Microbiome Resource at the Dalhousie University as previously described [
17].
2.8. Bioinformatics Analyses of Gut Microbiota
The fastq file containing raw gut microbiome data was demultiplexed and then was trimmed to remove primers. Trimmed reads were imported into Quantitative Insights into Microbial Ecology 2 (QIIME2) platform. Diversity metrics (Core-metrics-phylogenetic) within QIIME2 were used to evaluate α- and β-diversity of gut microbiota [
17].
2.9. Analysis of Fecal SCFAs
Fecal fatty acids were extracted using propyl chloroformate and derivatized with a reaction system containing propanol, pyridine and water. The analysis of SCFAs was conducted on an Agilent 7890A gas chromatography coupled with an Agilent 5975A inert XL EI/CI mass spectrometry in Microbiome Insights (Vancouver, BC, Canada) through customer service as previously described [
18].
2.10. Extraction and Metabolomics Sample Analysis of WR
WLD and WHR were milled separately and passed through a 60 mesh filter. Rice powders (100 mg in triplicates) were extracted with 2 mL of methanol: water (3:2, v/v) and vortexed for 2 min, and then sonicated for 40 min at 4˚C. After a centrifugation (14,000 rpm) for 10 minutes at 4˚C using HERMLE 16/018 rotor (Hermle Labortechnik GmbH, Siemensstr, Wehingen, Germany), the supernatant was dried under nitrogen and stored at -80°C. Prior to injections, dried extract were reconstituted with 200 µL of water: acetonitrile (4:1, v/v) containing 300 ng norvaline. Analyses were conducted in a high performance liquid chromatography (HPLC) system (1290 Infinity Agilent Ltd., Santa Clara, CA) coupled with a 6538 UHD Accurate LC-Quadrupole Time-Of-Flight (QTOF)-mass spectrometry (MS) (Agilent Technologies, Santa Clara, CA) with dual electro-spray ionization (ESI) source. A 2.1 mm x 100 mm Agilent ZORBAX SB-Aq column (Agilent Technologies) was maintained at 60°C for chromatographic separation of samples using water (A) and acetonitrile (B) containing 0.1% formic acid. Run time for the reactions was set at 10 minutes with gradient of 0±6 min 2% B; 6±8.50 min 60% B; 8.50±8.60 min 2% B and 8.60±10 min 2% B were used to re-equilibrate the column. For each individual sample analysis, 2 μL of extract was injected with a flow rate maintained at 0.7 mL/min.
2.11. Metabolomics Data Acquisition
Metabolomic data acquisitions were completed in both positive (+) and negative (-) ESI modes. MS spectra were collected within the range of 50±1700 m/z and known references masses of 121.0508 and 922.0097 (ESI+) and 112.9860 and 1033.9880 (ESI-) utilized during all runs. The liquid chromatography-QTOF-MS metabolomics data workflow followed standard protocols and utilized multiple algorithms incorporated in Agilent Mass Hunter Qualitative (MHQ, 7.01 and Mass Profiler Professional software through custom service in the Canadian Centre for Agri-Food Research in Health and Medicine [
10]. Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to assign potential physiological functions.
2.12. Statistics
Quantitative data was presented as means±standard deviation (SD). Probabilities from data >2 groups were analyzed using the one-way analysis of variance assay (ANOVA), followed by Kruskal-Wallis test paired with post-hoc Tukey test or Pairwise Wilcoxon test for probability between two groups. Liner discriminant analysis Effect Size (LEfSe) of gut microbiota was analyzed using Galaxy module. Correlations between two sets of variables were assessed using linear regression analysis. OriginPro 2021 software was used for plotting and statistical analysis. Probability <0.05 were considered as significant.
4. Discussion
The results of the present study demonstrated that the supplementation of WLD to HF diet dose-dependently attenuated circulatory metabolic and pro-inflammatory markers in mice compared to those receiving HF+WHR diet. The abundances of probiotic L. gasseri species bacteria in the feces of mice fed with HF+WLD diet were significantly greater than that in mice fed with HF+WHR or HF+WTWD diet. The abundance of fecal propionic acid in HF+WLD diet-fed mice was significantly higher than HF+WHR diet-fed mice. The fecal abundances of both L. gasseri and propionic acid negatively correlated with cholesterol, pro-inflammatory cytokines, FPG or HOMA-IR in the mice. The results of metabolomic analysis demonstrated that WLD contained substantially higher abundances of metabolites with capacity to modulate glucose metabolism compared to WHR, including AMP, homoarginine, Glc-1, 6-BP, Glu-MI and L-glutamate.
Previous studies by Moghadasian et al. demonstrated that WLD treatment induced attenuation of atherosclerosis was associated with increased abundances of fecal bacteria, including
Lactobacillus, and fatty acids, including butyric acid or longer chain fatty acids, in LDL receptor-deficient mice [
19]. Hou et al reported that Asian wild rice (Zizania latifolia) intake reduced liver steatosis, insulin resistance and gut dybiosis in HF diet-fed mice [
20]. The present study for the first time examined the effects of WLD on SCFAs including acetic acid and propionate acid and relationship between SCFA and gut microbiome in HF diet-fed mice. The findings of the present study demonstrated that WLD diet significantly increased the abundances of fecal propionic acid, and the relative abundance of fecal propionic acid negatively correlated with metabolic and inflammatory markers in HF-diet induced insulin resistant mice.
The results of LEfSe analysis indicated the feces of WLD diet-fed mice was enriched with
Lactobacillaceae family and
Lactobarcillus genus bacteria. The abundances of
L. gasseri species bacteria in the feces of mice fed with WLD diet were significantly higher than that in mice receiving WHR or WTWD diet. Although the negative correlation between the abundances of fecal
L. gasseri and propionate acid did not reach a significant level, but both
L. gasseri and propionate acid negatively correlated to FPG, cholesterol and MCP-1 in the mice. Several species of
Lactobacillus, including
L. gasseri, have been recognized as potent probiotics for the management of diabetes and obesity in clinical trials [
21,
22]. The findings of the present study suggest that the increases of the abundance of
L. gasseri in the gut may contribute to metabolic benefits and anti-inflammatory effects of WLD. Relationship between
L. gasseri and the production of propionate acid is warrant to be investigated.
In addition, the abundances of fecal Bacteroides genus bacteria and Ruminococcaceae family bacteria in WLD diet-fed mice were significantly higher than that in WHR diet-fed mice. The abundances of fecal isobutyric acid positively correlated with Bacteroidetes phylum bacteria in the mice, but there was no detectable correlation between fecal isobutyric acid and metabolic or inflammatory variables in mice.
Previous studies demonstrated that fiber-rich diet reduced FPG and insulin resistance in T2D patients, which was associated with enhancement of a variety of gut bacteria capable to generate SCFAs and modulate glucose metabolism [
23,
24]. Our group previously reported that fiber-rich brown rice and germinated brown rice attenuated FPG, lipids and inflammatory markers, in HF diet-fed mice [
13]. WLD contains 3-times higher level of fiber than that in WHR [
25,
26]. The precise role of fiber in the beneficial effect of WLD diet on metabolism, inflammation, gut microbiota and SCFA production are warrant to be further clarified in subsequent experiments.
WLD contains near 20-fold higher levels of AMP compared to WHR. AMP is a known agonist of AMP-activated protein kinase (AMPK). AMPK plays crucial modulatory roles in glucose, lipid metabolism, insulin resistance, and inflammation [
27]. Previous study demonstrated that the colonization
of Lactobacillus rhamnosus GG upregulated the phosphorylated AMPK in colon of mice [
28]. Propionic acid reduced gluconeogenesis in HepG2 hepatocytes via the activation of AMPK pathway [
29]. Propionic acid inhibited insulin induced
de novo lipogenesis and increased glucose uptake in primary rat adipocytes [
30]. Treatment of human adipose tissue explants with propionic acid resulted in downregulations of inflammatory cytokines and upregulations of lipoprotein lipase and glycose transporter [
31]. Besides, WLD also contains abundant amounts of L-glutamate [
32], AMPA [
33], KAPA [
34,
35], homoarginine [
36,
37], and Glc-1, 6-BP [
38,
39] which are involved in the regulation of glucose metabolism or insulin resistance, pyridoxamine or vitamin 6 [
40,
41], RGG [
42] and hyaluronic acid [
43], which were known to benefit in the management of diabetes, inflammation or diabetic complications. Those metabolites enriched in WLD may partially contribute to the anti-diabetic and anti-inflammatory benefits of WLD in HF diet induced insulin resistant mice.
The anti-insulin resistance and anti-inflammatory effects of WLD in HF diet-fed mice were consistent with previous reports [
10]. The present study added a half-dosage of WLD (WTWD) diet in the dietary regimen to determine the effectiveness of a lower and more feasible dosage of WLD on metabolism, inflammation, gut microbiota and SCFA production in mice. The results demonstrated that the WTWD diet, containing half-dosage of WLD compared to the WLD diet, significantly reduced triglycerides, total cholesterol, insulin resistance, pro-inflammatory cytokines and monocyte adhesion compared to WHR diet. The metabolic and anti-inflammatory effects of the WTWD diet were relatively weaker than the WLD diet as expected. The results suggest that weaker, but effective, metabolic and anti-inflammatory benefits may be achieved using a lower dosage of WLD (representing 10% of daily calorie intake) supplemented in HF diet in mice, which may be useful for future regimen design in human trials.
The limitation of the present study includes that the effects of WLD on the metabolism, inflammation and gut microbiota were only assessed in male, but not in female, mice. Future investigations may be required to compare the effects of WLD in female animals. The present study did not investigate the biological activities of compounds enriched in WLD. The effect and regimen for using WLD in diabetic patients need to be verified in randomized controlled clinical trials.
In conclusion, WLD supplementation in HF diet reduced FPG, lipids, insulin, HOMA-IR, monocyte adhesion and inflammatory cytokines related to monocyte adhesion compared to that in mice receiving WHR diet, which was associated with increased abundances of probiotic bacteria, L. gasseri, and propionic acid in feces of mice receiving WLD diet. The abundance of fecal L. gasseri and propionic acid negatively correlated with the metabolic and pro-inflammatory cytokines in the peripheral circulation of the mice. Multiple compounds enriched in WLD may independently modulate glucose metabolism and inflammation. The findings suggest that oral administration of WLD may improve metabolism and mitigate chronic inflammation in mice through multiple pathways, including, but not limited to, the modulation of gut microbiota, SCFA production and AMPK activation in high fat diet fed mice.
Abbreviations Used
AMP: adenosine monophosphate.
AMPA: aminomethylphosphonic acid.
AMPK: AMP-activated protein kinase.
ANOVA: analysis of variance assay.
ELISA: enzyme-linked immunosorbent assay.
FPG: fasting plasma glucose.
G1c-1: 6-BP: α-D-Glucose 1, 6-biphosphate.
GlcNAc-1P: N-Acetyl-alpha-D-glucosamine 1-phosphate.
Glu-MI: 6-(alpha-D-Glucosaminyl)-1D-myo-inositol.
HBSS: Hank’s balance salt solution.
HF: high fat (diet).
HOMA-IR: Homeostatic model assessment-insulin resistance.
KAPA: 8-amino-7-oxononanoate.
MCP-1: monocyte chemotactic protein-1.
PAI-1: plasminogen activator inhibitor-1.
RGG: α-L-rhamnopyranosyl-(1-2)-β-D-galactopyranosyl-(1)-β-D-glucopyranoside.
SCFA: short chain fatty acid.
TNFα: tumor necrosis factor-α.
WHR: white rice.
WLD: wild rice.
WTWD: 13g % of WHR and 13g % of WLD supplemented in HF diet.