Next Article in Journal
Graptopetalum paraguayense Extract Ameliorates Proteotoxicity in Aging and Age-Related Diseases in Model Systems
Next Article in Special Issue
Pre-Administration of Berberine Exerts Chemopreventive Effects in AOM/DSS-Induced Colitis-Associated Carcinogenesis Mice via Modulating Inflammation and Intestinal Microbiota
Previous Article in Journal
Pesticide Toxicity Footprints of Australian Dietary Choices
Previous Article in Special Issue
Modulation of Adipocyte Metabolism by Microbial Short-Chain Fatty Acids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Use of Mushrooms and Spirulina Algae as Supplements to Prevent Growth Inhibition in a Pre-Clinical Model for an Unbalanced Diet

1
Institute of Biochemistry and Nutrition, The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot 7610001, Israel
2
Koret School of Veterinary, The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem, Rehovot 7610001, Israel
*
Author to whom correspondence should be addressed.
Nutrients 2021, 13(12), 4316; https://doi.org/10.3390/nu13124316
Submission received: 3 November 2021 / Revised: 26 November 2021 / Accepted: 27 November 2021 / Published: 29 November 2021
(This article belongs to the Special Issue Dietary Bioactives, Gut Microbiota, and Human Health)

Abstract

:
Today’s eating patterns are characterized by the consumption of unbalanced diets (UBDs) resulting in a variety of health consequences on the one hand, and the consumption of dietary supplements in order to achieve overall health and wellness on the other. Balanced nutrition is especially crucial during childhood and adolescence as these time periods are characterized by rapid growth and development of the skeleton. We show the harmful effect of UBD on longitudinal bone growth, trabecular and cortical bone micro-architecture and bone mineral density; which were analyzed by micro-CT scanning. Three point bending tests demonstrate the negative effect of the diet on the mechanical properties of the bone material as well. Addition of Spirulina algae or Pleurotus eryngii or Agaricus bisporus mushrooms, to the UBD, was able to improve growth and impaired properties of the bone. 16SrRNA Sequencing identified dysbiosis in the UBD rats’ microbiota, with high levels of pro-inflammatory associated bacteria and low levels of bacteria associated with fermentation processes and bone related mechanisms. These results provide insight into the connection between diet, the skeletal system and the gut microbiota, and reveal the positive impact of three chosen dietary supplements on bone development and quality presumably through the microbiome composition.

1. Introduction

Malnutrition refers to deficiencies (undernutrition), excesses (over-nutrition), or imbalances in a person’s intake of energy and/or nutrients [1]. Micro-nutrient deficiency is another form of malnutrition which until recently was only associated with undernutrition, however today it is well established as a deficiency that can also happen in overweight and obese people. This phenomenon occurs due to consumption of energy-dense, micronutrient-poor foods and drinks and named: the hidden hunger [2].
The popular diets produced and consumed for the past 50 years are no longer nutritionally optimal; they typically contain high amounts of calories, sugar, fat and additives and are depleted in dietary fibers, essential amino acids and micronutrients [3]. This kind of eating pattern, is nowadays common also in low- and middle-income countries, which are also suffering from hunger conditions [4].
The deleterious effects of unbalanced diet (UBD) consumption are well established, and include obesity, dyslipidemia, heart disease, high blood pressure, and impact on daily performance, behavior and mood [5,6]. Yet, there are not enough studies discussing the effect of such a diet on skeletal development and quality [7,8]. Many factors including genes, ethnicity, gender, hormones, physical activity level and diet have a significant impact on bone development, quality, and peak bone mass (PBM) achievement with sexual maturation [9]. Twenty to forty percent of the variance in PBM is determined by lifestyle choices; thus, although heredity is the major contributor, optimization of lifestyle factors is an important strategy to maintain bone health in later life [9].
The effect of diet on skeletal tissue and bone health may be mediated through the microbial communities. The gut microbiota (GM) is acquired at birth and can be modulated by a number of environmental factors such as antibiotic treatments, physiologic state, diseases and of course diet [10]. Alterations in the microbiome have already been associated with variety of conditions, such as inflammatory bowel diseases (IBD), obesity, cardiovascular disease [10,11], and bone health [12]. So far, four potential mechanisms have been suggested as potential linkers between the GM and the skeleton: (1) nutrient absorption and vitamins biosynthesis, (2) the immune system, (3) microbe associated molecular patterns (MAMPs) across the endothelial barrier and (4) signaling molecules/Hormones alteration [13].
In addition to the increased consumption of UBD, the consumption of dietary supplements is in a constant rise as well; according to the Council for Responsible Nutrition report from 2019, 77% of Americans, more than 170 million people, informing they consume dietary supplements [14]. People are interested in easy and natural ways to improve their health status. Among the dietary supplements are vitamins, minerals, protein, amino acids, herbs and botanical [14,15].
Spirulina algae (SP) are filamentous cyanobacteria belonging to the Oscillactoriaceae family. SP is mainly known by its high content of protein, providing all the essential amino acids with a great bioavailability. Additionally, two essential fatty acids are found in the algae: linoleic (18:2 n-6) and γ-linolenic (18:3 n-6), which are very limited in the common diet. SP is also considered as good source for a variety of vitamins, minerals and functional ingredients such as phycocyanin, β-carotene and polysaccharides [16,17,18].
Nutrient composition is relatively similar in different mushrooms. Pleurotus eryngii (PE) and Agaricus bisporus (AB), commonly known as the king oyster mushroom and the white button mushroom, respectively, are two of the most widely cultivated mushrooms [19]. They are low in fat, containing protein rich in essential amino acids and fibers. Polysaccharides are the main part of the carbohydrate fraction of the mushrooms and have a variety of beneficial health properties [19,20]. Mushrooms are also good sources of several vitamins, particularly from the B complex and minerals as phosphorus, potassium and magnesium. The bioactive substances found in mushrooms are polysaccharides, mainly β-glucans and secondary metabolites as polyphenols [20,21].
Lowering blood glucose, improving plasma lipid concentration and antioxidant activity are just a few of the health benefits associated with algae and medicinal mushrooms [22,23,24,25]. In spite of the rising popularity and the intensive research of these supplements, the impact on bone development and quality has not been investigated yet.
In the present study, we established a model that mimics the modern UBD and tested its effects on the growth and skeletal development in young rats. Moreover, we evaluated the impact of botanical dietary supplements, Spirulina algae and Pleurotus eryngii and Agaricus bisporus mushrooms when added to this diet.

2. Materials and Methods

2.1. Animals

Forty-eight three weeks old Female Sprague Dawley (SD) post weaning rats, were purchased from Harlan Laboratories (Rehovot, Israel) and housed under environmentally controlled conditions. All procedures were approved by the Hebrew University Animal Care Committee (permit number AG-18-15530-2). After 5 days of adaptation to control chow diet, we randomly divided the rats into six groups, with eight rats in each group: (1) Control (CON); (2) Unbalanced diet (UBD); (3) Unbalanced diet + Spirulina algae (UBD + SP); (4) Unbalanced diet + Pleurotus eryngii mushroom (UBD + PE); (5) Unbalanced diet + Agaricus bisporus mushroom (UBD + AB); and (6) Unbalanced diet + SP + PE + AB (UBD + MIX). At the end of the experiment, after 6 weeks, the rats were anesthetized with isoflurane and blood and caecum samples were collected and stored at −80 °C. Next, their femora, tibias, and vertebra were harvested. The femora and lumbar vertebrae were manually cleaned of soft tissue, wrapped in saline soaked gauze and stored at −20 °C until analysis (mechanical/micro-CT testing). Tibiae were fixed immediately after sacrifice for histological analysis.

2.2. Diet Preparation and Composition

During the postweaning period rats from the CON group were fed a diet based on the recommendation of the American Institute of Nutrition (AIN-93G) that was formulated for the growth phase of rodents: 16% fat, 63.5% carbohydrate and 20.5% protein with the addition of multi vitamins and multi minerals based on AIN-93G as well (Table 1) [26]. The UBD groups were fed a diet high in fat and low in protein, vitamins and minerals, to mimic the popular modern unbalanced diet: 25% fat, 65% carbohydrate, and 10% protein with only 50% of the AIN-93G recommendation for vitamins and minerals for rodents in growth phase (Table 2). The entire diet was homogenized, shaped as dumplings and frozen at −20 °C. The UBD group was the basis for the supplemented diets, with either Spirulina algae, Pleurotus eryngii mushroom, Agaricus bisporus mushroom or all three. The supplements were added to the UBD as dry powder proportionally to the daily intake recommendations for humans (Table 3). The UBD + MIX diet contain the combined dosage of all three supplements.

2.3. Histological Staining of Growth-Plate (GP) Sections

Safranin-O staining was used to examine the tibial GPs. Tibia samples were fixed overnight in 4% paraformaldehyde (PFA, Sigma, St Louis, MO, USA) at 4 °C followed by 3 weeks of decalcification in 0.5 M EDTA pH 7.4. Afterwards, the samples were dehydrated by a 1-hour wash with increased concentrations of ethanol (70%, 85%, 95% and 2 × 100%), transferred into Histo-Clear (Bar-Naor, Israel) twice for 2 h each and embedded in paraffin blocks. By using a Leica (Agentec, Yakum, Israel) microtome, we prepared 5 μm transverse tissue sections for histological staining [27].
The sections were deparaffinized by heating the slides, followed by two washes with xylen, rehydration by washes with decreasing ethanol concentrations (2 × 100%, 95%, 85%, and 70%) and a final wash with distilled water. Weigert’s iron hematoxylin solution, a fast- green solution and acetic acid were used for Safranin-O staining. The sections were dried and mounted with DPX mounting for histology.

2.4. Imaging and Measurement of GPs

Stained transverse sections of tibiae were viewed by an eclipse E400 Nikon light microscopy with ×4, ×10, ×20 or ×40 objectives, using light filters. Images were taken by a high -resolution camera (Olympus DP 71) controlled by Cell A software (Olympus, Japan). The thickness of the entire GPs, the proliferative zone (PZ) and the hypertrophic zone (HZ) as well as the number of cells were measured using the Cell A software with a measuring tool feature. Measurements were performed on these sections from 5 different animals from each group. In each slide, 10 random locations throughout the GPs were selected and measured [28,29].

2.5. Bone Microarchitecture

Right femora were scanned using a Skyscan 1174 (Skyscan, Bruker, Belgium) X-ray computed micro-tomography device. Images were obtained at 50 kV X-ray tube voltage and 800 μA current, using a 0.25 mm aluminum filter, 4000 ms exposure time, and 15 μm optical resolution. For each specimen, a series of 900 projection images were obtained (a rotation step of 0.4°, averaging 2 frames, for a total 360° rotation). A stack of 2-D X-ray shadow projections was reconstructed to obtain images using NRecon software (Skyscan). Next, images were subjected to morphometric analysis using CTAn software (CT Analyser 1.13.5.1, Skyscan). Morphometric parameters were calculated as suggested by recent guidelines for bone microstructure assessment. To analyze the diaphyseal cortical region, 200 slices, centered at the mid diaphysis, equivalent to 2.764 mm, were chosen. Global grayscale threshold levels for the cortical region were between 71 and 255. For the trabecular region, a total of 150 slices, equivalent to 2.073 mm of the bone, were selected, and adaptive grayscale threshold levels between 58 and 255 were used. Two phantoms with known density (0.25 and 0.75 g/cm3) were scanned under the same conditions of the femora samples allowing to measure the cortical diaphysis BMD (bone mineral density); quantification were carried out using CTAn software [8,30].
The 3rd–5th lumbar vertebrae were scanned and analyzed as well. The spatial resolution was 18 μm and the total rotation was 180°, except that all other parameters were identical to the femoral scans. The region chosen for the analysis was manually selected and consisted of 120 slices of the 5th vertebra starting from the proximal end-plate. Adaptive grayscale threshold levels between 66 and 255 were selected for the analysis of the trabecular vertebral region.
The length of the femora and 3rd–5th segments of the lumbar vertebrae were measured using the Micro-CT device prior to the scans. By using the Amira software (v.6.4, FEI, Hillsboro, OR, USA) reconstructed scans were volume-rendered to visualize the 3D morphology and BMD variation using visualization, of the selected sample from each group.

2.6. Three Point Bending for Bone Mechanical Analysis

Right femora were tested using an Instron mechanical tester (Model 3345). Each bone was placed within a custom- built saline containing testing chamber and on two supports having rounded profiles (2 mm in diameter), so that the supports were in touch with the posterior aspect of the diaphysis. The distance between the stationary supports was set to 10 mm, to ensure that the relatively tubular portion of the mid-diaphysis rests on these supports. A pronged loading device was applied to the anterior surface of the bones, precisely in the middle between the two supports. First, an initial preload of 0.1 N was applied to hold the bone in place; following that, the prong was advanced at a constant rate of 600 µm/min, loaded up to the fracture point, identified by a sudden >20% decrease in load [8]. Force-displacement data were collected by the Instron software BlueHill (version 2.0, Instron Corporation, Norwood, Massachusetts, USA) at 10 Hz. The resulting force-displacement curves were used to calculate bone stiffness, bone yield point, load of fracture, maximal load and area under the curve was measured to calculate the total energy to fracture (E to F) [31].

2.7. Macro Ad Micro-Nutrient Analysis

Macro nutrient analysis of the three supplements was performed by analytical laboratory AMINOLAB ((Ness Ziona, Israel). Analyses of fatty acids (FAs), amino acids (AAs) and minerals composition were conducted by the Interdepartmental Equipment Unit, The Robert H. Smith Faculty of Agriculture, Food and Environment, The Hebrew University of Jerusalem (Israel). Qualitative analysis of FAs was performed by gas chromatography (GC); liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used for AAs quantitative analysis after hydrolysis of the samples; and last, Inductively Coupled Plasma Optical Emission spectroscopy (ICP-OES) and Inductively Coupled Plasma Mass Spectrometry (ICP-MS) were performed after microwave-assisted nitric acid digestion of the samples in order to assess the mineral content in the supplements.

2.8. Caecum Sample Collection and DNA Extraction

At the day of the sacrifice, caecum from all 48 rats were collected and stored at −80 °C until the day of the extration. 0.25 gr of defrosted caecum content were extracted using QIAamp PowerFecal kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions.

2.9. 16SrRNASequencing

16SrRNA Sequencing were performed by the University of Illinois at Chicago Core for Research Informatics (UICCRI), using single-end reads with a length of 272 bp with coverage of ~38,500 per sample. Sequencing trimming was done by cutadapt, Amplicon Sequence Variants (ASVs) were inferred using DADA2 and taxonomy was classified against the SILVA 132_16S database. A biom files for each experiment was received.

2.10. Bioinformatics Analysis

Bioinformatics analysis of genus level was conducted by the Bioinformatics Unit, Weizmann Institute of Science (Rehovot, Israel). The threshold for significant differential expression (DE) was set at: padj for multiple testing by Benjamin-Hochberg procedure ≤ 0.05, |log2FoldChange| ≥ 1 and max count ≥ 30 [32]. Principal coordinate analysis (PCoA) was calculated by using Bray-Curtis distances matrix. PCoA is an ordination analysis that simplifies high dimensional data into 2D plot. The axes of the plot are ranked in order of importance: differences along the first principal component axis (PC1) are more important than differences along the second principal component axis (PC2). In order to evaluate the diversity in each of the tested caecum samples Alpha diversity was calculated using the Shannon index. Differential abundance analysis was carried out using Deseq2 on the counts produced by Washington’s lab [33].

2.11. Statistical Analysis

All data are expressed as mean ± SD. The significance of differences between groups was determined using JMP 14.0.0 Statistical Discovery Software (SAS Institute 2000) by one-way analysis of variance. Differences between groups were further evaluated by a Tukey-Kramer HSD test and considered significant at p < 0.05.

3. Results

In order to evaluate the influence of unbalanced diet (UBD) and different nutritional supplements (Spirulina algae, Pleurotus eryngii mushroom and Agaricus bisporus mushroom) added to this diet, on skeletal development, a 6-week long experiment was conducted on SD rats after weaning (3 weeks old). This time line is mimicking the human growth period until sexual maturation [8].

3.1. The Positive Effect of Nutritional Supplements on Growth Pattern

The effect of the diets on growth was first determined in terms of body weight and length. Body weight represents whole body growth rate, including energy balance, whereas body length is an indicator for longitudinal bone growth.
The UBD group gained less weight compared to the CON group, while the supplemented groups, UBD + SP, UBD + PE and UBD + AB, gained more weight than the UBD group and managed to partially rescue growth inhibition. The UBD + MIX group showed the highest weight gain from all UBD groups, and by the end of the experiment did not differ from the CON groups (Figure 1A; Supplementary Table S1 for statistical differences). Body length was measured from the tip of the nose to the end of the tail. Results showed the same patterns as in body weight—UBD rats were shorter than the CON rats, and by the end of the experiment differ also from the UBD supplemented groups. No statistical difference was found between UBD + AB, UBD + MIX and CON groups in the last length measurement (Figure 1B,D; Supplementary Table S2 for statistical differences). Femur and 3rd–5th lumbar vertebrae length have been measured in the end of the experiment as well. All four supplemented groups demonstrated higher results in these measurements compared with the UBD group (Figure 1E,F).
Rats from the UBD group consumed fewer calories than all the other. This modification happened subsequent to the slowdown in their growth. Deceleration in growth occurred at the 5th day, while the decline in caloric intake started on day 15th of the experiment (Figure 1C; Supplementary Figure S1). A slowdown in growth was also observed in the supplemented groups through the same period; however, rats did not reduce their caloric intake which allows them to subsequently improve their growth. Thus, rats which received supplementation consumed more calories than rats consuming UBD, throughout the experiment, and demonstrated growth patterns approximating the CON group (Figure 1; Supplementary Figure S1).
Since longitudinal bone growth originates from the GP [34], and abnormal impact of the UBD on bone elongation was observed, we performed histological analysis. Safranin-O staining which stains the cartilage in a noticeable red-orange was performed. We found co-responding results with the rats’ length. Both GP thickness and total number of cells were negatively affected by the unbalanced diet and improved thanks to the different supplementations (Figure 2).

3.2. The Positive Effect of Nutritional Supplements on Bone Quality

Consumption of a balanced diet is one of the main environmental factors necessary for adequate skeletal development and fulfilment of its full genetic potential [35]. Some of the critical properties encompassed by the term bone quality are trabecular architecture, cortical geometry and mechanical integrity [36]. Micro-CT and three-point bending techniques were used to estimate these parameters.
Analysis of the femora trabecular bone demonstrated major alterations to bone microarchitecture in rats from the UBD group. Lower values of BV/TV% and Tb.N and higher values of Tb.Sp were found in rats from the UBD group compared with those from the CON group. Twenty-five to thirty-seven percentage of improvement occurred in these parameters in the supplemented groups. Tb.Th that was also negatively affected by the UBD consumption was fully rescued by the supplements to the CON level (Table 4; Figure 3).
Vertebral bone is for the most part a trabecular bone compartment, and surprisingly, analysis of the 5th vertebra yielded milder effects of the diets (Table 4). These differences could be implied by a different growth pattern between the spine and the limbs, or by different mechanical forces subjected between these two bones.
Cortical bone analysis revealed low values in parameters such as Ct.Ar/Tt.Ar and Cs.Th in the UBD groups which were improved in the supplemented groups. BMD, which reflects the mineralization level of the bones, showed lower values in rats consuming UBD. The UBD + AB and UBD + MIX diets improved the mineral density of the cortical bone the most, reaching to the same level as the CON group (Table 4; Figure 3). The skeleton plays a critical structural role in bearing functional loads, and failure to do so results in fractures; thus, it is valuable to evaluate the mechanical properties of the bones. We found that all tested mechanical parameters were lower in bones of rats fed UBD. The linear slope which represents the stiffness of the bent object and E to F which defined as the work that must be done in order to fracture the bone were lower in the UBD group compared to the CON group, and tended to be improve in the supplemented groups. Moreover, yield point (load at which the load-deflection curve ceases to be linear, indicating micro-damage accumulation in the bone material) and maximal load were lower in the UBD group than in the CON group, and partially improved in the supplemented groups (Table 4).

3.3. Diet Composition Analysis of the Supplements

To further comprehend the mechanism by which the supplements commence their beneficial effect, diet composition was investigated [37]. The results in Figure 4 displayed the percentage of carbohydrates, proteins and fats found in the supplements. The main component of Spirulina is protein and of mushrooms carbohydrate. Fat fraction in the three supplements is the smallest, and almost undetectable in Spirulina.

3.3.1. Treatment of Amino Acids and Mineral Deficiencies

Amino acids content was evaluated in the three supplements by LC-MS/MS [38]. Except tryptophan, which destructs during hydrolysis all essential AA were found within the nutritional supplements (Supplementary Table S3). The AA content in the two mushrooms is quite similar, with some advantage to PE, while greatest amounts were found in SP powder for both essential and non-essential AA (Supplementary Table S3). Although the supplements contained all the essential AA, calculation of AA content in kg diet relative to the required minimum presented by the nutritional research council (NRC) [39] (values presented in %), showed poor values in the supplemented groups (Supplementary Table S4). When calculation of the essential AA consumed by the rats throughout the experiment was done; higher amounts in the supplemented groups as compared with the UBD group were found, most probably due to their higher food intake. According to the NRC requirements, only in phenylalanine + tyrosine (tyrosine may supply up to 50 percent of phenylalanine and thus calculated together) rats in SP, AB and MIX supplementations groups reached the necessary minimum. Comparison to the CON group still demonstrated major deficiencies (Figure 5A).
Evaluation of the mineral content in SP, PE and AB powders was performed by ICP-OES and ICP-MS techniques [40]. Compared to the mushrooms, spirulina was found being rich specially in sodium, chloride, calcium and magnesium. The highest amount of potassium was found in PE mushroom, while AB was abundant in iron. None lethal trace of lead, tin, cadmium and mercury were found in the three supplements as well (Supplementary Table S5). The essential mineral content in kg diet relative to the required minimum presented by the NRC (values presented in %) was calculated and as in AA, none of the dietary supplements have been able to improve the deficiencies caused by the UBD. Levels above the requested minimum according to NRC were in sodium, chloride and zinc, regardless of the added powders (Supplementary Table S6). When comparing the amount of total essential minerals consumed by the rats during all 6 weeks of experiment, similar to the amino acids: the supplemented groups consumed more minerals compare with the UBD group, however still could not reach the CON or the NRC levels. Zinc is the only mineral that the supplemented rats consumed more of it compared to both UBD and NRC (Figure 5B).

3.3.2. Functional Fatty Acids

Analysis of the supplement and the main source of fat in the diets, soy oil, were made by GC [41]. The results suggest a resemblance between the soy oil and the two mushrooms, with 18:2 cis 9,12 being the main FA followed by 18:1 cis 9 and 16:0. Both 16:0 and 18:2 cis 9,12 were also found in Spirulina; however, in a different ratio, 16:0 was the main FA detected and it represents about 50% of Spirulina FA. 18:3 cis 6,9,12 and 16:1 cis 9 were also found in spirulina, while 18:3 cis 6,9,12 being the second most abundant FA (Figure 6). The presence of these FA in the diet, although in minute amounts may generate modifications in the cellular or microbiome level.

3.4. Microbiome Analysis

In order to assess the specific changes in the rats’ intestinal microbiota, 16SrRNA sequencing was preformed from caecum samples of 9 weeks old rats.
Alpha diversity was calculated using the Shannon index in order to evaluate the diversity in each of the tested caecum samples. Significantly decreased was found in the UBD group compared with the CON group. Supplementation of PE, AB or Mix was able to significantly increase the alpha diversity; however only trends of improvement were found with SP supplementation (Figure 7A).
Principal coordinate analysis (PCoA) was used with the aim to evaluate the ecological distance between samples and produce a graphical configuration in a low-dimension. Figure 7B presents PCoA plot based on the Bray-Curtis distances matrix, which was calculated on the genus level. The calculated distances in the CON samples was much smaller; the group is mostly clustered together and separated from the other experimental groups. All four supplementation groups led to a smaller variance in the microbiome composition inside the group and between the groups as well. Samples from these groups are more clustered together compared with the UBD group but not as in the CON and separation between the supplemented groups and the CON group was still observed.
On an effort to point out the bacteria that might take part in mediating the effect of UBD consumption on the skeletal system we focused on the significantly bacteria changes in the UBD compared with the CON group and examined the influence the different supplements had on the abundance of these bacteria. Genera abundance was first assessed using a heat map (Figure 8A).
The genera Mucispirillum and Lachnospiraceae_NK4A136_group differed significantly between the UBD and the CON groups: increased and decreased respectively with UBD consumption. We have not seen any changes in these bacteria in the supplemented groups compare to both the UBD and the CON groups.
We saw an opposite effect in the UBD and the supplemented groups in five different bacteria: Butyricimonas, Parabacteroides, Rhodospirillales /other, Parasutterella and Escherichia /Shigella which were elevated in the UBD compared with the CON group, and were reduced in the supplemented groups compared with the UBD. In addition reduced levels of Acetatifactor were detected in the UBD compared with the CON group and the supplements managed to improve this reduction (Figure 8B). Supplementation of SP, PE and AB to the UBD decreased the abundance of Butyricimonas and Rhodospirillales/other in a significant way compared to the UBD without supplements, and no statistical difference compared to the CON group was observed (Figure 8B). We also saw reduction in Parabacteroides and Parasutterella bacteria compared to the UBD while administration of the different supplements was performed. The abundance of these two bacteria was even higher compared to the CON group with PE and AB mushrooms addition to the diet (Figure 8B). All four supplements elevated the levels of Acetatifactor compared to the UBD alone, and MIX supplementation was the only one that cause to a significant difference compared to the CON group as well (Figure 8B). Tendency to higher abundance of Escherichia /Shigella in the UBD compare to the CON group was found and supplementation with SP, PE and MIX lowered this effect in a significant way (Figure 8B). These results suggest that some of the supplements were able to equalize the abundance of different bacteria to CON level, while others had a partial effect.

4. Discussion

In this research we establish a model that simulates the massive changes in the general population food habits and mimics the diet prevalent in both developed and developing countries. We shed a light on the importance of balanced nutrition, and revealed that consumption of unbalanced diet lead to severe growth retardation, disrupted bone structural parameters, decreased bone mineral density (BMD) and reduced strength and rigidity. Low doses of botanical dietary supplements were able to rescue the damage caused by this diet; leading to partial or full rescue to the control levels.
To elucidate the foundation for the beneficial effect the supplements have on the rats’ growth and skeletal quality, several explanations were tested. First, the nutritional composition of the supplements was examined to see if it compensates as for the deficiencies in the UBD. Despite the abundant ingredients in algae and the mushrooms, their addition to the diet did not quantitatively compensate for the huge deficiencies. Due to the minute amount added to the UBD diet (as recommended by manufacturers) the supplements failed to increase the amount of essential amino acids and minerals in the diet. Even when the increased food consumption of the supplemented groups was taken into account the amounts of essential nutrients did not reach the minimal requirements according to the NRC [39]. The dietary supplements do not improve the rats’ growth and skeletal quality by reversing the shortages of the unbalanced diet; suggesting they function as growth enhancers rather than nutritional improvers.
Next we examined whether the supplements effect might be related to their impact on the GM. Rise attention had been recently given to the effects microbial communities have on the human body and its overall health; and to the negative health consequences associated with microbial imbalance [13]. This imbalance could be attributed to the gain or loss of community members, known as dysbiosis. Dysbiosis, is commonly associated with impaired gut barrier function, increased pathogen translocation and inflammatory cell activation, and has been directly linked to several pathological conditions such as: inflammatory bowel disease (IBD), irritable bowel syndrome (IBS), diabetes, obesity and cancer [42,43]. High fat diets, sugar rich-diet and high red meat consumption are examples of diets already found to induce dysbiosis [44]. Bone metabolism is tightly regulated by the immune system, both locally and systemically. Macrophage, dendritic cells and T-helper cell (Th) subpopulations, produce pro-inflammatory cytokines such as TNFα, TGFβ, Type I interferon (IFNγ), IL-6, 17 and 22, known to regulate bone growth, and promotes resorption over formation [45].
Moreover, metabolic damage was found in dysbiosis due to reduced vitamins biosynthesis, impaired nutrient absorption and enzymatic breakdown of compounds as polysaccharides and polyphenols. The colon receives indigestible carbohydrates and proteins that represent 10–30% of the total ingested energy. Without the activity of the colonic microbiota these nutrients would generally be eliminated via the stool [46]. Short chain fatty acids (SCFA) are the most investigated guts products, that have been linked to T cells regulation [47], improved calcium absorption [48] and IGF-1 secretion, a dominant bone growth factor [49]. Indole, fermentation product of aromatic amino acid tryptophan, stimulates the production of mucin and the proliferation of goblet cells, to maintain intestinal impermeability [50]. It is well established that vitamins from B group and vitamin K can be synthesized by the GM and be used by the host [51]. Thus, decreased bacterial diversity, as found in the UBD group may act negatively on the skeleton.
When focusing on the specific changes occurring in the bacteria’s levels following administration of an unbalanced diet and the supplements, we found several bacteria associated with inflammation to be higher in the UBD group. For instance, Mucispirillum which is associated with active colitis [52]. Actinobacteria and Proteobacteria phylum are associated with IBD. Escherichia/Shigella, genus from the Proteobacteria phylum is one of the main bacteria associated with Crohn’s disease [43] and Parasutterella, a relatively new genus also part of the Proteobacteria phylum with advanced IBD [53]. Higher levels of these bacteria were found in the UBD group and reduction occurred in the group consuming the supplements. Up to 50% of IBD patients experience at least one extra-intestinal manifestation; among them is the involvement of reduced bone density [54]. It appears that alerted bone profile in IBD can be caused by number of reasons: inflammatory process in the gut with high levels of cytokines, nutritional alterations and nutrients deficiency due to inadequate diet intake and/or malabsorption [55].
Higher levels of Parabacteroides species, that were elevated by consumption of UBD in our experiments, were found in osteoporosis and osteopenia patients [56]. Osteopenia and osteoporosis are two bone diseases often associated with ageing. Although we used young rats in our experiments, their bone profile and skeletal characterization were comparable to those suffering from osteoporosis; reduced skeletal strength and BMD along with microstructural defects. Interestingly, SP, PE and MIX supplementation were able to decrease the elevated levels of Parabacteroides in the rats’ caecum.
Decreased levels of Lachnospiraceae_NK4A136_group and Acetatifactor were detected in the UBD group. Previous studies have indicated that Lachnospiraceae _NK4A136_group may have an anti-inflammatory effect [57]. It is positively correlated with acetic and butyric acid, two forms of SCFA, SOD, GSH-Px and GSH [58]. Acetatifactor metabolite products can activate the G protein-coupled bile acids (BAs) receptor (TGR5) signaling, which has a few beneficial functions in the skeletal system [59]. First, TGR5 can promote the secretion of glucagon-like peptide-1 (GLP-1) [59,60], that was recently associated with bone tissue. GLP-1 receptor (GLP-1 R) has been found on immature osteoblastic cell lines, on mouse osteoblast cells and in bone marrow stem cells [61]. It was found to increase the number of osteoblasts on the surface of trabecular bone and up regulate the expression of Runx2, ALP, collagen type 1, and osteocalcin [62]. GLP-1 may also promote bone formation through the reduction of sclerostin expression, which results in increased Wnt signaling pathway and bone formation [61,62]. Suppression of bone resorption by GLP-1 is associated with the expression of GLP1R in thyroid C cells where it can promote the secretion of calcitonin, an inhibitor of bone resorption by osteoclasts [63]. TGR5 signaling can also increase the activity of type 2 deiodinase (D2) which converts thyroxine (T4) into the active form Triiodothyronine (T3) [59]. T3 acts in chondrocytes and osteoblasts to regulate intramembranous and edonchondral ossification and control the rate of linear growth and bone maturation. Studies suggest that T3 regulates a number of key growth factor signaling pathways including IGF-1, FGF, PTHrP, Ihh and Wnt to influence skeletal growth [64]. In all four supplemented groups, improvement in Acetatifactor abundance was observed. Hormonal testes we preformed reveled significant decrease in T4 in the UBD group, which improved thanks to the supplements. No differences were found in TSH, T3 and IGF-1 (data not shown).
Lastly, we evaluated the active ingredients or fermentation products found in the supplements. Spirulina contains several active ingredients as phycocyanin, β-carotene, γ-linolenic fatty acid and omega-3 fatty acids; all of which show antioxidant, immunomodulatory and anti-inflammatory activities [24,65]. Spirulina also contains polysaccharides, resistant starch and oligosaccharides which known as dietary fibers. Humans lack the ability to digest β 1→4 linkages so the undigested materials reach the large intestine and undergo fermentation to SCFA, which has a variety of health beneficial as described earlier [65].
The bioactive substances found in mushrooms are mainly divided into: β-glucans, types of polysaccharides, and secondary metabolites as polyphenols [20]. Current literature suggests that long-term consumption of diets rich in polyphenols protects against osteoporosis and play a role in bone metabolism. They found to be protective against oxidative stress and capable of down regulating osteoclast differentiation and activity by inhibition of RANKL and COX-2 production and stimulation of OPG [66]. Das, Mukherjee, and Mitra et al. showed that black tea extract prevents bone loss induced by ovariectomy (OVX) in rats. The OVX group had significantly higher serum ALP, TRAP and lower bone density, which all were improved with the supplementation [67]. Bone regeneration effects of β-glucans as well as its anti-osteoporotic properties have been studied using in vitro and animal experiments. Increased BMD and calcium bioavailability were observed in rat model for osteoporosis. Human-based clinical trials are limited and reveled beneficial effects when added to healthy menopausal women: significant changes in serum bone-specific ALP, serum osteocalcin, urinary cross-linked N-telopeptide of type-1 collagen, urinary cross-linked C-telopeptide of type-1 collagen, calcium, and phosphorus levels [68,69].
In this study, we created for the first time a model that simulates the massive changes in the general population food habits: higher consumption of fat and insufficient amount of protein, vitamins and minerals. We investigated the diet effect as intact and not as separate components; the primary advantage of this approach is the translational potential from model organisms to humans.
Dramatic effect to the bone was observed along with slowdown in the rats’ growth rate. Supplementation of the deficient diet with botanical supplements showed surprising results: significant improvement to the morphological and mechanical bone parameters and growth patterns similar to those of the control group. Since we show that the mushrooms and spirulina were unable to correct the deficiencies of the diet, we suggest that the major positive effect of the supplements were attributed to their effect on the microbiome or thanks to their active ingredients and fermentation products.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/nu13124316/s1. Table S1: Statistical differences of body weight presented in graph X. Table S2: Statistical differences of body length presented in graph X. Figure S1: Energy intake and weight gain per week. Table S3: Amino acids content in the three supplements. Table S4: % Amino acids per kg diet from the NRC. Table S5: Minerals content in the three supplements. Table S6: % Minerals per kg diet from the NRC.

Author Contributions

Conceptualization, R.S. (Roni Sides) and E.M.-O.; Methodology, R.S. (Roni Sides), A.S., S.G.-F., J.Z., R.K.-A., R.R., O.B., S.P. and E.M.-O.; Formal analysis, R.S. (Roni Sides) and M.D.; Resources, E.M.-O.; writing—original draft preparation, R.S. (Roni Sides) and E.M.-O.; writing—review and editing, R.S. (Roni Sides), E.M.-O. and R.S. (Ron Shahar); supervision, E.M.-O., R.S. (Ron Shahar) and M.D.; project administration, S.P.; funding acquisition, E.M.-O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Ministry of agriculture, “Tevat Noah- Nitzan” grant number 12-05-0027.

Institutional Review Board Statement

All procedures were approved by the Hebrew University Animal Care Committee (permit number AG-18-15530-2).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. UNICEF; WHO. The International Bank for Reconstruction and Development/The World Bank. In Levels and Trends in Child Malnutrition: Key Findings of the 2019 Edition of the Joint Child Malnutrition Estimates; WHO: Geneva, Switzerland, 2019. [Google Scholar]
  2. Gödecke, T.; Stein, A.J.; Qaim, M. The global burden of chronic and hidden hunger: Trends and determinants. Glob. Food Sec. 2018, 17, 21–29. [Google Scholar] [CrossRef]
  3. Lucas, T.; Horton, R. The 21st-century great food transformation. Lancet 2019, 393, 386–387. [Google Scholar] [CrossRef]
  4. Branca, F.; Demaio, A.; Udomkesmalee, E.; Baker, P.; Aguayo, V.M.; Barquera, S.; Dain, K.; Keir, L.; Lartey, A.; Mugambi, G.; et al. A new nutrition manifesto for a new nutrition reality. Lancet 2020, 395, 8–10. [Google Scholar] [CrossRef]
  5. Lustig, R.H. Ultraprocessed Food: Addictive, Toxic, and Ready for Regulation. Nutrients 2020, 12, 3401. [Google Scholar] [CrossRef]
  6. Moubarac, J.-C. Ultra-Processed Foods in Canada: Consumption, Impact on Diet Quality and Policy Implications; University of Montreal: Montreal, QC, Canada, 2017. [Google Scholar] [CrossRef]
  7. Aslam, M.N.; Jepsen, K.J.; Khoury, B.; Graf, K.H.; Varani, J. Bone structure and function in male C57BL/6 mice: Effects of a high-fat Western-style diet with or without trace minerals. Bone Rep. 2016, 5, 141–149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Zaretsky, J.; Griess-fishheimer, S.; Carmi, A.; Shmul, T.T.; Ofer, L.; Sinai, T.; Penn, S. Ultra-processed food targets bone quality via endochondral ossi fi cation. Bone Res. 2021, 9, 14. [Google Scholar] [CrossRef] [PubMed]
  9. Bonjour, J.-P.; Theintz, G.; Law, F.; Slosman, D.; Rizzoli, R. Peak bone mass. Osteoporos. Int. 1994, 4, S7–S13. [Google Scholar] [CrossRef] [PubMed]
  10. Lozupone, C.A.; Stombaugh, J.I.; Gordon, J.I.; Jansson, J.K.; Knight, R. Diversity, stability and resilience of the human gut microbiota. Nature 2012, 489, 220–230. [Google Scholar] [CrossRef] [Green Version]
  11. Fan, Y.; Pedersen, O. Gut microbiota in human metabolic health and disease. Nat. Rev. Microbiol. 2021, 19, 55–71. [Google Scholar] [CrossRef] [PubMed]
  12. Ohlsson, C.; Sjögren, K. Effects of the gut microbiota on bone mass. Trends Endocrinol. Metab. 2015, 26, 69–74. [Google Scholar] [CrossRef] [PubMed]
  13. Christakos, S.; Veldurthy, V.; Patel, N.; Wei, R.; Ramsey, W.; Isales, C.M.; Lavoie, B.; Chishti, S.; Parameswaran, N.; Britton, R.A.; et al. Understanding the Gut-Bone Signaling Axis; Advances in Experimental Medicine and Biology Series; McCabe, L.R., Parameswaran, N., Eds.; Springer International Publishing: Cham, Switzerland, 2017; Volume 1033, ISBN 978-3-319-66651-8. [Google Scholar]
  14. The Council for Responsible Nutrition. Who Takes Dietary Supplements? And Why? The Council for Responsible Nutrition: Washington, DC, USA, 2019. [Google Scholar]
  15. Kamiński, M.; Kręgielska-Narożna, M.; Bogdański, P. Determination of the Popularity of Dietary Supplements Using Google Search Rankings. Nutrients 2020, 12, 908. [Google Scholar] [CrossRef] [Green Version]
  16. Matos, J.; Cardoso, C.; Bandarra, N.M.; Afonso, C. Microalgae as healthy ingredients for functional food: A review. Food Funct. 2017, 8, 2672–2685. [Google Scholar] [CrossRef] [PubMed]
  17. de Morais, M.G.; da Silva Vaz, B.; de Morais, E.G.; Costa, J.A.V. Biologically Active Metabolites Synthesized by Microalgae. Biomed Res. Int. 2015, 2015, 835761. [Google Scholar] [CrossRef] [Green Version]
  18. Hoseini, S.M.; Khosravi-Darani, K.; Mozafari, M.R. Nutritional and Medical Applications of Spirulina Microalgae. Mini-Rev. Med. Chem. 2013, 13, 1231–1237. [Google Scholar] [CrossRef] [PubMed]
  19. Reis, F.S.; Barros, L.; Martins, A.; Ferreira, I.C.F.R. Chemical composition and nutritional value of the most widely appreciated cultivated mushrooms: An inter-species comparative study. Food Chem. Toxicol. 2012, 50, 191–197. [Google Scholar] [CrossRef] [Green Version]
  20. Valverde, M.E.; Hernández-Pérez, T.; Paredes-López, O. Edible Mushrooms: Improving Human Health and Promoting Quality Life. Int. J. Microbiol. 2015, 2015, 376387. [Google Scholar] [CrossRef] [PubMed]
  21. Mattila, P.; Könkö, K.; Eurola, M.; Pihlava, J.-M.; Astola, J.; Vahteristo, L.; Hietaniemi, V.; Kumpulainen, J.; Valtonen, M.; Piironen, V. Contents of Vitamins, Mineral Elements, and Some Phenolic Compounds in Cultivated Mushrooms. J. Agric. Food Chem. 2001, 49, 2343–2348. [Google Scholar] [CrossRef]
  22. Iñiguez, M.; Pérez-Matute, P.; Villanueva-Millán, M.J.; Recio-Fernández, E.; Roncero-Ramos, I.; Pérez-Clavijo, M.; Oteo, J.-A. Agaricus bisporus supplementation reduces high-fat diet-induced body weight gain and fatty liver development. J. Physiol. Biochem. 2018, 74, 635–646. [Google Scholar] [CrossRef]
  23. Jeong, S.C.; Jeong, Y.T.; Yang, B.K.; Islam, R.; Koyyalamudi, S.R.; Pang, G.; Cho, K.Y.; Song, C.H. White button mushroom (Agaricus bisporus) lowers blood glucose and cholesterol levels in diabetic and hypercholesterolemic rats. Nutr. Res. 2010, 30, 49–56. [Google Scholar] [CrossRef] [PubMed]
  24. Wu, Q.; Liu, L.; Miron, A.; Klímová, B.; Wan, D.; Kuča, K. The antioxidant, immunomodulatory, and anti-inflammatory activities of Spirulina: An overview. Arch. Toxicol. 2016, 90, 1817–1840. [Google Scholar] [CrossRef]
  25. Kulshreshtha, A.; Zacharia, A.J.; Jarouliya, U.; Bhadauriya, P.; Prasad, G.B.K.S.; Bisen, P.S. Spirulina in health care management. Curr. Pharm. Biotechnol. 2008, 9, 400–405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Reeves, P.G. Committee Report AIN-93 Purified Diets for Laboratory Rodents: Final Report of the American Institute of Nutrition Ad Hoc Writing Committee on the Reformulation of the AIN-76A Rodent Diet. J. Nutr. 1939, 123, 1939–1951. [Google Scholar] [CrossRef] [PubMed]
  27. Koren, N.; Simsa-maziel, S.; Shahar, R.; Schwartz, B.; Monsonego-ornan, E. ScienceDirect Exposure to omega-3 fatty acids at early age accelerate bone growth and improve bone quality. J. Nutr. Biochem. 2014, 25, 623–633. [Google Scholar] [CrossRef] [PubMed]
  28. Simsa-maziel, S.; Zaretsky, J.; Reich, A.; Koren, Y.; Shahar, R.; Monsonego-ornan, E. IL-1RI participates in normal growth plate development and bone modeling. Am. J. Physiol. Endocrinol. Metab. 2013, 305, 15–21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Idelevich, A.; Kerschnitzki, M.; Shahar, R.; Monsonego-Ornan, E. 1,25(OH)2D3 Alters Growth Plate Maturation and Bone Architecture in Young Rats with Normal Renal Function. PLoS ONE 2011, 6, e20772. [Google Scholar] [CrossRef]
  30. Bouxsein, M.L.; Boyd, S.K.; Christiansen, B.A.; Guldberg, R.E.; Jepsen, K.J.; Müller, R. Guidelines for assessment of bone microstructure in rodents using micro-computed tomography. J. Bone Miner. Res. 2010, 25, 1468–1486. [Google Scholar] [CrossRef] [PubMed]
  31. Rozner, R.; Vernikov, J.; Griess-Fishheimer, S.; Travinsky, T.; Penn, S.; Schwartz, B.; Mesilati-Stahy, R.; Argov-Argaman, N.; Shahar, R.; Monsonego-Ornan, E. The Role of Omega-3 Polyunsaturated Fatty Acids from Different Sources in Bone Development. Nutrients 2020, 12, 3494. [Google Scholar] [CrossRef] [PubMed]
  32. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. Ser. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  33. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Kronenberg, H.M. Developmental regulation of the growth plate. Nature 2003, 423, 332–336. [Google Scholar] [CrossRef] [PubMed]
  35. Rizzoli, R. Nutritional aspects of bone health. Best Pract. Res. Clin. Endocrinol. Metab. 2014, 28, 795–808. [Google Scholar] [CrossRef] [PubMed]
  36. Cole, J.H.; van der Meulen, M.C.H. Whole Bone Mechanics and Bone Quality. Clin. Orthop. Relat. Res. 2011, 469, 2139–2149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Maynard, L.A. The atwater system of calculating the caloric value of diets. J. Nutr. 1944, 28, 443–452. [Google Scholar] [CrossRef]
  38. Le, A.; Ng, A.; Kwan, T.; Cusmano-Ozog, K.; Cowan, T.M. A rapid, sensitive method for quantitative analysis of underivatized amino acids by liquid chromatography–tandem mass spectrometry (LC–MS/MS). J. Chromatogr. B 2014, 944, 166–174. [Google Scholar] [CrossRef]
  39. Subcommittee on Laboratory Animal Nutrition; Committee on Animal Nutrition; Board on Agriculture; National Research Council. Nutrient Requirements of Laboratory Animals, 4th ed.; National Academies Press: Washington, DC, USA, 1995; ISBN 978-0-309-05126-2. [Google Scholar]
  40. Gonzalez, M.H.; Souza, G.B.; Oliveira, R.V.; Forato, L.A.; Nóbrega, J.A.; Nogueira, A.R.A. Microwave-assisted digestion procedures for biological samples with diluted nitric acid: Identification of reaction products. Talanta 2009, 79, 396–401. [Google Scholar] [CrossRef] [PubMed]
  41. Lu, J.; Argov-Argaman, N.; Anggrek, J.; Boeren, S.; van Hooijdonk, T.; Vervoort, J.; Hettinga, K.A. The protein and lipid composition of the membrane of milk fat globules depends on their size. J. Dairy Sci. 2016, 99, 4726–4738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Clemente, J.C.; Ursell, L.K.; Parfrey, L.W.; Knight, R. The Impact of the Gut Microbiota on Human Health: An Integrative View. Cell 2012, 148, 1258–1270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Belizário, J.E.; Faintuch, J. Metabolic Interaction in Infection; Springer International Publishing: Cham, Switzerland, 2018; Volume 109, ISBN 978-3-319-74931-0. [Google Scholar]
  44. Yang, Q.; Liang, Q.; Balakrishnan, B.; Belobrajdic, D.P.; Feng, Q.J.; Zhang, W. Zhang Role of Dietary Nutrients in the Modulation of Gut Microbiota: A Narrative Review. Nutrients 2020, 12, 381. [Google Scholar] [CrossRef] [Green Version]
  45. Watrous, D.A.; Andrews, B.S. The metabolism and immunology of bone. Semin. Arthritis Rheum. 1989, 19, 45–65. [Google Scholar] [CrossRef]
  46. Krajmalnik-Brown, R.; Ilhan, Z.-E.; Kang, D.-W.; DiBaise, J.K. Effects of Gut Microbes on Nutrient Absorption and Energy Regulation. Nutr. Clin. Pract. 2012, 27, 201–214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Smith, P.M.; Howitt, M.R.; Panikov, N.; Michaud, M.; Gallini, C.A.; Bohlooly-Y, M.; Glickman, J.N.; Garrett, W.S. The Microbial Metabolites, Short-Chain Fatty Acids, Regulate Colonic Treg Cell Homeostasis. Science 2013, 341, 569–573. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Wallace, T.C.; Marzorati, M.; Spence, L.; Weaver, C.M.; Williamson, P.S. New Frontiers in Fibers: Innovative and Emerging Research on the Gut Microbiome and Bone Health. J. Am. Coll. Nutr. 2017, 36, 218–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Yan, J.; Herzog, J.W.; Tsang, K.; Brennan, C.A.; Bower, M.A.; Garrett, W.S.; Sartor, B.R.; Aliprantis, A.O.; Charles, J.F. Gut microbiota induce IGF-1 and promote bone formation and growth. Proc. Natl. Acad. Sci. USA 2016, 113, E7554–E7563. [Google Scholar] [CrossRef] [Green Version]
  50. Postler, T.S.; Ghosh, S. Understanding the Holobiont: How Microbial Metabolites Affect Human Health and Shape the Immune System. Cell Metab. 2017, 26, 110–130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Rowland, I.; Gibson, G.; Heinken, A.; Scott, K.; Swann, J.; Thiele, I.; Tuohy, K. Gut microbiota functions: Metabolism of nutrients and other food components. Eur. J. Nutr. 2018, 57, 1–24. [Google Scholar] [CrossRef] [Green Version]
  52. Rooks, M.G.; Veiga, P.; Wardwell-Scott, L.H.; Tickle, T.; Segata, N.; Michaud, M.; Gallini, C.A.; Beal, C.; van Hylckama-Vlieg, J.E.; Ballal, S.A.; et al. Gut microbiome composition and function in experimental colitis during active disease and treatment-induced remission. ISME J. 2014, 8, 1403–1417. [Google Scholar] [CrossRef] [PubMed]
  53. Chiodini, R.J.; Dowd, S.E.; Chamberlin, W.M.; Galandiuk, S.; Davis, B.; Glassing, A. Microbial population differentials between mucosal and submucosal intestinal tissues in advanced Crohn’s disease of the ileum. PLoS ONE 2015, 10, e0134382. [Google Scholar] [CrossRef]
  54. Sgambato, D.; Gimigliano, F.; De Musis, C.; Moretti, A.; Toro, G.; Ferrante, E.; Miranda, A.; De Mauro, D.; Romano, L.; Iolascon, G.; et al. Bone alterations in inflammatory bowel diseases. World J. Clin. Cases 2019, 7, 1908–1925. [Google Scholar] [CrossRef] [PubMed]
  55. Agrawal, M.; Arora, S.; Li, J.; Rahmani, R.; Sun, L.; Steinlauf, A.F.; Mechanick, J.I.; Zaidi, M. Bone, Inflammation, and Inflammatory Bowel Disease. Curr. Osteoporos. Rep. 2011, 9, 251–257. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, J.; Wang, Y.; Gao, W.; Wang, B.; Zhao, H.; Zeng, Y.; Ji, Y.; Hao, D. Diversity analysis of gut microbiota in osteoporosis and osteopenia patients. PeerJ 2017, 5, e3450. [Google Scholar] [CrossRef] [Green Version]
  57. Li, J.; Wu, T.; Li, N.; Wang, X.; Chen, G.; Lyu, X. Bilberry anthocyanin extract promotes intestinal barrier function and inhibits digestive enzyme activity by regulating the gut microbiota in aging rats. Food Funct. 2019, 10, 333–343. [Google Scholar] [CrossRef] [PubMed]
  58. Li, H.; Liu, F.; Lu, J.; Shi, J.; Guan, J.; Yan, F.; Li, B.; Huo, G. Probiotic Mixture of Lactobacillus plantarum Strains Improves Lipid Metabolism and Gut Microbiota Structure in High Fat Diet-Fed Mice. Front. Microbiol. 2020, 11, 512. [Google Scholar] [CrossRef] [Green Version]
  59. Jia, E.; Liu, Z.; Pan, M.; Lu, J.; Ge, Q. Regulation of bile acid metabolism-related signaling pathways by gut microbiota in diseases. J. Zhejiang Univ. B 2019, 20, 781–792. [Google Scholar] [CrossRef]
  60. Thomas, C.; Gioiello, A.; Noriega, L.; Strehle, A.; Oury, J.; Rizzo, G.; Macchiarulo, A.; Yamamoto, H.; Mataki, C.; Pruzanski, M.; et al. TGR5-Mediated Bile Acid Sensing Controls Glucose Homeostasis. Cell Metab. 2009, 10, 167–177. [Google Scholar] [CrossRef] [Green Version]
  61. Schiellerup, S.P.; Skov-Jeppesen, K.; Windeløv, J.A.; Svane, M.S.; Holst, J.J.; Hartmann, B.; Rosenkilde, M.M. Gut Hormones and Their Effect on Bone Metabolism. Potential Drug Therapies in Future Osteoporosis Treatment. Front. Endocrinol. 2019, 10, 75. [Google Scholar] [CrossRef]
  62. Zhao, C.; Liang, J.; Yang, Y.; Yu, M.; Qu, X. The Impact of Glucagon-Like Peptide-1 on Bone Metabolism and Its Possible Mechanisms. Front. Endocrinol. 2017, 8, 98. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Yavropoulou, M.; Yovos, J. Incretins and bone: Evolving concepts in nutrient-dependent regulation of bone turnover. Hormones 2013, 12, 214–223. [Google Scholar] [CrossRef] [Green Version]
  64. Van Der Eerden, B.C.J.; Karperien, M.; Wit, J.M. Systemic and Local Regulation of the Growth Plate. Endocr. Rev. 2003, 24, 782–801. [Google Scholar] [CrossRef] [Green Version]
  65. Wells, M.L.; Potin, P.; Craigie, J.S.; Raven, J.A.; Merchant, S.S.; Helliwell, K.E.; Smith, A.G.; Camire, M.E.; Brawley, S.H. Algae as nutritional and functional food sources: Revisiting our understanding. J. Appl. Phycol. 2017, 29, 949–982. [Google Scholar] [CrossRef] [PubMed]
  66. Luka, Đ.; Fu, A.; Muhvi, D.; Giacometti, J. The role of polyphenols on bone metabolism in osteoporosis. Food Res. Int. 2015, 77, 290–298. [Google Scholar] [CrossRef]
  67. Das, A.S.; Mukherjee, M.; Mitra, C. Evidence for a prospective anti-osteoporosis effect of black tea (Camellia Sinensis) extract in a bilaterally ovariectomized rat model. Asia Pac. J. Clin. Nutr. 2004, 13, 210–216. [Google Scholar]
  68. Choi, J.S.; Park, M.Y.; Kim, J.D.; Cho, H.R.; Choi, I.S.; Kim, J.W. Safety and efficacy of polycalcium for improving biomarkers of bone metabolism: A 4-week open-label clinical study. J. Med. Food 2013, 16, 263–267. [Google Scholar] [CrossRef] [Green Version]
  69. Muhammad, K.; Bashir, I.; Choi, J. Clinical and Physiological Perspectives of β-Glucans: The Past, Present, and Future. Int. J. Mol. Sci. 2017, 18, 1906. [Google Scholar] [CrossRef]
Figure 1. Growth patterns and food consumption. (A) Body weight during the entire experiment (g). (B) Body length during the entire experiment (cm). (C) Energy intake (kcal/rat/day). (D) Body length in the end of the experiment (cm). (E) Femur length measured by micro-CT (mm). (F) 3rd–5th lumbar vertebra length measured by micro-CT (mm). Values are expressed as mean ± SD of n = 8 rats/group, different superscript letters are significantly different (p < 0.05) by one-way ANOVA followed by Tukey’s test.
Figure 1. Growth patterns and food consumption. (A) Body weight during the entire experiment (g). (B) Body length during the entire experiment (cm). (C) Energy intake (kcal/rat/day). (D) Body length in the end of the experiment (cm). (E) Femur length measured by micro-CT (mm). (F) 3rd–5th lumbar vertebra length measured by micro-CT (mm). Values are expressed as mean ± SD of n = 8 rats/group, different superscript letters are significantly different (p < 0.05) by one-way ANOVA followed by Tukey’s test.
Nutrients 13 04316 g001
Figure 2. Histological evaluation of the tibial growth plates (GPs) of 9 weeks old rats. Transverse tissue sections of 5 μm were prepared by microtome. GP width and number of cells has been measured in cell A software at 10 selected locations. (A) Safranin-O staining of the tibial GP. (B) GP width (µm); Number of cells. Values are expressed as mean ± SD of n = 5 rats/group, different superscript letters are significantly different (p < 0.05) by one-way ANOVA followed by Tukey’s test.
Figure 2. Histological evaluation of the tibial growth plates (GPs) of 9 weeks old rats. Transverse tissue sections of 5 μm were prepared by microtome. GP width and number of cells has been measured in cell A software at 10 selected locations. (A) Safranin-O staining of the tibial GP. (B) GP width (µm); Number of cells. Values are expressed as mean ± SD of n = 5 rats/group, different superscript letters are significantly different (p < 0.05) by one-way ANOVA followed by Tukey’s test.
Nutrients 13 04316 g002
Figure 3. Three-dimensional reconstructions of femora scaled for mineral density using AMIRA software. Sagittal sections of the distal femora bone; Transverse sections of cortical area at mid-diaphysis; Transverse sections of trabecular bones at the distal epiphysis above the growth plate.
Figure 3. Three-dimensional reconstructions of femora scaled for mineral density using AMIRA software. Sagittal sections of the distal femora bone; Transverse sections of cortical area at mid-diaphysis; Transverse sections of trabecular bones at the distal epiphysis above the growth plate.
Nutrients 13 04316 g003
Figure 4. Macronutrient’s composition. Caloric distribution percentage of macronutrients in the three supplements.
Figure 4. Macronutrient’s composition. Caloric distribution percentage of macronutrients in the three supplements.
Nutrients 13 04316 g004
Figure 5. Essential amino acids and minerals consumption. (A) Grams of essential amino acids consumed by the rats during the entire experiment. (B) Grams of essential minerals for skeletal development consumed by the rats during the entire experiment. The NRC-CON group is based on the CON food intake and the NRC requires recommendations.
Figure 5. Essential amino acids and minerals consumption. (A) Grams of essential amino acids consumed by the rats during the entire experiment. (B) Grams of essential minerals for skeletal development consumed by the rats during the entire experiment. The NRC-CON group is based on the CON food intake and the NRC requires recommendations.
Nutrients 13 04316 g005
Figure 6. Fatty acid profiles. Lipids were analyzed by gas chromatography–mass spectroscopy.
Figure 6. Fatty acid profiles. Lipids were analyzed by gas chromatography–mass spectroscopy.
Nutrients 13 04316 g006
Figure 7. Bacterial diversity between and within samples. (A) Alpha diversity was calculated using Shannon index. (B) Principal coordinates analysis was calculated based on bray-curtis distances matrix. p-value was calculated using Pairwise comparisons Wilcox. Test.
Figure 7. Bacterial diversity between and within samples. (A) Alpha diversity was calculated using Shannon index. (B) Principal coordinates analysis was calculated based on bray-curtis distances matrix. p-value was calculated using Pairwise comparisons Wilcox. Test.
Nutrients 13 04316 g007
Figure 8. Change in the bacterial profile. (A) Results were analyzed using Deseq2 and assessed using heatmap. (B) Differential abundance analysis was carried out using Deseq2 on the counts products. p value was adjusted for multiple testing by Benjamin-Hochberg procedure (Padj < 0.05). (*) Significant difference compared with the CON group. (#) Significant difference compared with the UBD group.
Figure 8. Change in the bacterial profile. (A) Results were analyzed using Deseq2 and assessed using heatmap. (B) Differential abundance analysis was carried out using Deseq2 on the counts products. p value was adjusted for multiple testing by Benjamin-Hochberg procedure (Padj < 0.05). (*) Significant difference compared with the CON group. (#) Significant difference compared with the UBD group.
Nutrients 13 04316 g008
Table 1. Macronutrient, vitamin and mineral composition of the experimental diets. Values are presented as gram per dry weight.
Table 1. Macronutrient, vitamin and mineral composition of the experimental diets. Values are presented as gram per dry weight.
Control DietUnbalanced Diet
Ingredientg\kg Dietg\kg Diet
Cornstarch397438.7
Casein (≥85% protein) 200105
Dextrinized cornstarch (90–94% tetrasaccharides)132145.7
Sucrose 100110
Soybean oil70118.8
Fiber5054
Mineral mix (AIN-93G-MX)3518.9
Vitamin mix (AIN-93-VX)105.4
L-Cystine31.6
Choline bitartrate (41.1% choline)2.51.34
Tert-butylhydroquinone0.0140.015
Table 2. Caloric composition of the experimental diets. Values are presented as calories per weight with added water.
Table 2. Caloric composition of the experimental diets. Values are presented as calories per weight with added water.
Control DietUnbalanced Diet
IngredientKcal\kg Diet% KcalKcal\kg Diet% Kcal
Cornstarch140263.51559.465
Dextrinized cornstarch (90–94% tetrasaccharides)465.5518
Sucrose352.7392
Casein (≥85% protein)70520.537410
L-Cystine10.65.6
Soybean oil555.516949.925
Total calories3491.3 3799
Kcal/gr3.49 3.799
Table 3. Supplements added to the diet. Values are presented as gram per dry weight and as Kcal per dry weight of the supplements added to kg diet.
Table 3. Supplements added to the diet. Values are presented as gram per dry weight and as Kcal per dry weight of the supplements added to kg diet.
Ingredientg\kg DietKcal\kg Diet
Spirulina algae519
Pleurotus eryngii 3.515.765
Agaricus bisporus3.513.618
Table 4. Morphometric and mechanic characteristics.
Table 4. Morphometric and mechanic characteristics.
CONUBDUBD + SPUBD + PEUBD + ABUBD + MIX
Femora trabecular bone microarchitecture
BV/TV (%)38.222 ± 2.81 a23.261 ± 2.79 c30.44 ± 3.49 b30.263 ± 2.99 b33.214 ± 1.29 b33.697 ± 3.19 b
Tb.Th (mm)0.119 ± 0.009 a0.105 ± 0.005 b0.118 ± 0.007 a0.117 ± 0.006 a0.117 ± 0.04 a0.121 ± 0.008 a
Tb.N (1/mm)3.203 ± 0.15 a2.205 ± 0.19 c2.583 ± 0.2 b2.592 ± 0.19 b2.841 ± 0.16 b2.771 ± 0.16 b
Tb.Sp (mm)0.264 ± 0.06 c 0.572 ± 0.12 a 0.466 ± 0.06 ab 0.45 ± 0.05 b0.398 ± 0.07 b 0.378 ± 0.03 b
5th vertebra trabecular bone microarchitecture
BV/TV (%)36.797 ± 1.4 a30.568 ± 1.24 c32.841 ± 1.77 bc33.146 ± 2.41 bc35.011 ± 2.15 ab31.881 ± 2.57 c
Tb.Th (mm)0.128 ± 0.005 a0.111 ± 0.003 d0.116 ± 0.003 c0.115 ± 0.003 cd0.121 ± 0.003 b0.115 ± 0.003 cd
Tb.N (1/mm)2.871 ± 0.09 a2.762 ± 0.1 a2.84 ± 0.12 a2.869 ± 0.17 a2.901 ± 0.15 a2.775 ± 0.19 a
Tb.Sp (mm)0.297 ± 0.017 a0.328 ± 0.017 a0.312 ± 0.012 a0.307 ± 0.025 a 0.313 ± 0.03 a 0.322 ± 0.037 a
Femora cortical bone microarchitecture
Ct.Ar/Tt.Ar (%)48.925 ± 1.69 a40.378 ± 1.87 c43.534 ± 1.06 b42.945 ± 1.69 b44.6 ± 1.38 b44.565 ± 1.61 b
Cs.Th (mm) 0.466 ± 0.02 a0.343 ± 0.02 c0.384 ± 0.01 b0.385 ± 0.02 b0.387 ± 0.01 b0.392 ± 0.01 b
Ma.Ar (mm2)5.324 ± 0.4 a5.108 ± 0.59 a5.081 ± 0.25 a5.299 ± 0.35 a4.828 ± 0.3 a4.981 ± 0.39 a
BMD (g/cm3)1.26 ± 0.03 a1.184 ± 0.01 c1.205 ± 0.02 bc1.199 ± 0.02 bc1.229 ± 0.03 ab1.23 ± 0.04 ab
Femora bone mechanical properties
Stiffness (N/mm)268.72 ± 45.5 a195.6 ± 43 b225.32 ± 22.9 ab197.34 ± 44 b228.49 ± 33.9 ab202.57 ± 39.1 b
Yield point (N)45.55 ± 3.9 a26.31 ± 4.2 c36.01 ± 3.7 b34.64 ± 2.7 b37.8 ± 2.6 b37.28 ± 3 b
Fracture load (N)74.13 ± 7.5 a42.69 ± 11.4 b52.06 ± 6.2 b41.81 ± 9.3 b51.25 ± 6.5 b50.21 ± 11.2 b
Max load (N)90.47 ± 7.7 a63.04 ± 9.3 c75.32 ± 6.7 b68.81 ± 5.7 bc70.01 ± 5.9 bc73.26 ± 6.2 bc
E to F (N × mm)83.01 ± 16.5 a54.93 ± 10.6 c71.8 ± 11.8 abc72.04 ± 11.6 abc63.2 ± 8.7 bc74.84 ± 17.5 ab
Bones were scanned by Micro-CT to determine geometric parameters. After reconstructing 2D and 3D analyses were performed. Trabecular bone parameters: Bone volume over total volume, BV/TV (%); Trabecular thickness, Tb.Th (mm); Trabecular number, Tb.N (1/mm); Trabecular separation, Tb.Sp (mm). Cortical bone parameters: Cortical area fraction, Ct.Ar/Tt.Ar (%); Cortical thickness, Cs.Th (mm); Medullary area, Ma.Ar (mm2); Bone mineral density, BMD (mm3). Mechanical properties were evaluated using three-point bending experiment. Biomechanical parameters obtained from load–displacement curve: Stiffness (N/mm); Yield point (N); Fracture load (N); Max load (N); Energy to fracture, E to F (N × mm). Values are expressed as mean ± SD of n = 8 rats/group, different superscript letters are significantly different (p < 0.05) by one-way ANOVA followed by Tukey’s test.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sides, R.; Griess-Fishheimer, S.; Zaretsky, J.; Shitrit, A.; Kalev-Altman, R.; Rozner, R.; Beresh, O.; Dumont, M.; Penn, S.; Shahar, R.; et al. The Use of Mushrooms and Spirulina Algae as Supplements to Prevent Growth Inhibition in a Pre-Clinical Model for an Unbalanced Diet. Nutrients 2021, 13, 4316. https://doi.org/10.3390/nu13124316

AMA Style

Sides R, Griess-Fishheimer S, Zaretsky J, Shitrit A, Kalev-Altman R, Rozner R, Beresh O, Dumont M, Penn S, Shahar R, et al. The Use of Mushrooms and Spirulina Algae as Supplements to Prevent Growth Inhibition in a Pre-Clinical Model for an Unbalanced Diet. Nutrients. 2021; 13(12):4316. https://doi.org/10.3390/nu13124316

Chicago/Turabian Style

Sides, Roni, Shelley Griess-Fishheimer, Janna Zaretsky, Astar Shitrit, Rotem Kalev-Altman, Reut Rozner, Olga Beresh, Maïtena Dumont, Svetlana Penn, Ron Shahar, and et al. 2021. "The Use of Mushrooms and Spirulina Algae as Supplements to Prevent Growth Inhibition in a Pre-Clinical Model for an Unbalanced Diet" Nutrients 13, no. 12: 4316. https://doi.org/10.3390/nu13124316

APA Style

Sides, R., Griess-Fishheimer, S., Zaretsky, J., Shitrit, A., Kalev-Altman, R., Rozner, R., Beresh, O., Dumont, M., Penn, S., Shahar, R., & Monsonego-Ornan, E. (2021). The Use of Mushrooms and Spirulina Algae as Supplements to Prevent Growth Inhibition in a Pre-Clinical Model for an Unbalanced Diet. Nutrients, 13(12), 4316. https://doi.org/10.3390/nu13124316

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

Article Metrics

Back to TopTop