ORIGINAL RESEARCH article
Sec. Animal Nutrition and Metabolism
Volume 8 - 2021 | https://doi.org/10.3389/fvets.2021.654223
The Effects of 7 Days of Feeding Pulse-Based Diets on Digestibility, Glycemic Response and Taurine Levels in Domestic Dogs
- 1Department of Veterinary Biomedical Sciences, Saskatoon, SK, Canada
- 2Department of Food and Bioproduct Sciences, University of Saskatchewan, Saskatoon, SK, Canada
Grain-based carbohydrate sources such as rice comprise 30–50% of commercial pet foods. Some pet foods however have removed the use of grains and have instead incorporated pulses, such as peas and lentils, resulting in grain-free diets. The hypothesis was dog diets with higher levels of dietary fiber will produce a low glycemic response due to decreased rates of digestion and lowered bioavailability of all macronutrients and increased fecal bile salt excretion. This in turn was hypothesized to produce lower plasma concentrations of cysteine, methionine and taurine after 7 days of feeding each test diet in dogs. Six diets were formulated at an inclusion level of 20% available carbohydrate, using white rice flour (grain) or whole pulse flours from smooth pea, fava bean, red lentil or 2 different wrinkled pea varieties (CDC 4,140–4 or Amigold) and fed to beagles in a randomized, cross-over, blinded design. After 7 days feeding each diet, fasting blood glucose was the lowest in the lentil (3.5 ± 0.1 mmol/L) and wrinkled pea (4,140–4; 3.6 ± 0.1 mmol/L) diet periods, while peak glucose levels was lowest after feeding the lentil diet (4.4 ± 0.1 mmol/L) compared to the rice diet. Total tract apparent digestibility of all macronutrients as well as taurine differed among diets yet plasma taurine was not outside normal range. Decreased macronutrient and amino acid digestibility was associated with increasing amylose and dietary fiber content but the specific causative agent could not be determined from this study. Surprisingly, digestibility decreases were not due to increased bile salt loss in the feces since increasing dietary fiber content led to decreased fecal bile salt levels. In conclusion, although pulse-based canine diets have beneficial low glycemic properties, after only 7 days, these pulse-based diets decrease macronutrient and amino acid digestibility. This is likely related at least in part to the lower animal protein content, but on a long-term basis could put domestic dogs at risk for low taurine and dilated cardiomyopathy.
The global pet food industry has been steadily growing and is projected to reach a value of $91 billion USD by 2022 (1). In a majority of homes with pets, pet owners feed them commercially prepared diets as they are affordable and nutritionally complete with the belief that they promote animal health (2). Pet diets in North America are formulated to meet nutritional requirements based on the standards set by the Association of American Feed Control Officials (AAFCO).
Pet owners have an active role in choosing diets to feed their pets and generally make their decisions based on: (i) true knowledge of what comprises a healthy diet, (ii) perception of nutritional requirements, (iii) human diet trends and (iv) overall opinions of the pet food industry. Importantly, owners will often change their pet's diet to match the diet consumed by the human owner (3). Pet foods are marketed strategically to owners, with trendy claims such as “organic,” “natural” and “grain-free” often found with premium pet foods (4). Grain-free diets exclude the use of grains such as wheat, corn, or rice flours and instead incorporate pulses such as peas, lentils and fava beans as the major carbohydrate source (5). While pulses are common dietary ingredients for both humans and animals, they are considered to be highly nutritious primarily due to their high levels of protein, in addition to carbohydrate, fiber, vitamins and minerals (6). Pulse crops are slowly digested due to relatively high amylose, resistant starch and dietary fiber content (7, 8). This characteristic of pulses can be utilized to optimize satiety through a lowered glycemic response and glycemic index (6), a feature that is also in dogs [(9, 10), Briens et al., unpublished].
Postprandial blood glucose responses can vary based on the type of carbohydrate and the rate of digestion (11). In dogs, diets with different levels of starch cause variations in postprandial glucose responses (12). In human studies, foods with a low glycemic index, such as legumes or pulses, are slowly digested, which provides a slow and sustained release of glucose into the blood stream (11). Alternatively, foods that are quickly digested, such as white bread, have a high glycemic index, which result in rapid and high postprandial blood glucose levels (11). Similarly, in both dogs and cats, legumes result in a slow or negligible postprandial glucose response while grains result in a faster, higher postprandial blood glucose response [(9, 13), Briens et al., unpublished]. Thus, pulse-containing diets promote glucose control, insulin sensitivity, satiety, weight control and longer-term health in humans and dogs (10, 14).
While high-protein pulses can decrease the cost to produce a high protein pet food due to a lower need for animal protein, a downside to pulses is that plant protein lacks taurine. Moreover, compared to cereal grains, pulses contain limited amounts of the non-essential amino acid cysteine and reduced amounts of the essential amino acids, methionine and tryptophan (6). Cysteine and methionine are used by the dog liver and the central nervous system to synthesize taurine via the transsulfuration pathway. Thus, taurine is not considered to be an essential amino acid in dogs. While dilated cardiomyopathy is common in dogs (15), some cases are associated with low plasma taurine and can be reversed with taurine supplementation (16). Canine dilated cardiomyopathy was reported by the US Food and Drug Administration (FDA), based initially on 9 case studies, to be associated with feeding grain-free diets in July 2018. This led to a decrease in grain-free dog food sales, decreased use of pulses in pet food and losses to the pulse-growing agriculture sector. Confusion among veterinarians and pet owners as to whether grain-free diets are healthy for dogs was further exacerbated by a recent acknowledgment by the FDA in November 2020 that causes of dilated cardiomyopathy in dogs may be more complicated than just a single ingredient such as pulses and is instead likely multi-factorial. To begin to address some of these questions, this study aimed to explore whether grain-free diets lead to taurine deficiency and if this is associated with simultaneous low cysteine and methionine levels. Moreover, mechanisms by which pulses could deplete taurine, cysteine or methionine require experimental confirmation.
One possible explanation is dietary fiber which was linked to decreased taurine levels in dogs (17, 18). Fiber increases the excretion of fecal bile acids, and since taurocholate is the most common bile salt in dogs, this leads to increased fecal loss of taurine (19, 20). In addition to taurine loss, high dietary fiber decreases protein digestibility, resulting in decreased cysteine and methionine (precursors of taurine) bioavailability. Again, while the slow carbohydrate digestibility of pulses has some beneficial health effects such as low glycemic index, high dietary fiber of pulses may exacerbate the already low taurine, cysteine and methionine found in grain-free diets (20, 21). What is unclear is what component of dietary fiber is responsible for the increased fecal bile salt loss, decreased protein digestibility and decreased sulfur amino acid bioavailability. As a first exploration, this study aims to use diets with increasing levels of dietary fiber and amylose content using different pulses as well as to explore its role in these processes in dogs.
The hypothesis was that dog diets with higher levels of dietary fiber will produce a low glycemic response due to decreased rates of digestion and lowered bioavailability of all macronutrients and increased fecal bile salt excretion. This in turn is hypothesized to produce lower plasma concentrations of cysteine, methionine and taurine after 7 days of feeding each test diet in dogs. In order to investigate these hypotheses, whole and complete diets formulated to include 20% available carbohydrate using a grain (white rice flour) compared to whole pulse flours from smooth pea (CDC Inca), fava bean (CDC Snowdrop), red lentil (CDC Maxim) or 2 different wrinkled pea varieties (CDC 4,140–4 or Amigold) were fed to beagles in a randomized, cross-over, blinded design. After 7 days of feeding each diet, macronutrient digestibility and glycemic responses were examined along with plasma concentrations of cysteine, methionine, and taurine as well as fecal bile salt concentrations in beagle dogs.
All procedures and handling of the dogs were conducted following protocols approved by the University of Saskatchewan's Animal Research Ethics Board according to guidelines that were established by the Canadian Council on Animal Care (Animal Utilization Protocol #20190055).
Adult Beagle dogs (n = 8; 4 spayed females, 4 neutered males; 8.87 ± 0.90 kg, 2–4 years old) were obtained from certified scientific breeders (Marshal Bioresources, North Rose, NY, USA and King Fisher International, Stouffville, Ontario, Canada). Beagles were housed at the Animal Care Unit (ACU) in the Western College of Veterinary Medicine at the University of Saskatchewan, Saskatoon, SK, Canada. Beagles were group-housed during the day in a large enclosure to allow for daily socialization but were individually kenneled during feedings and overnight. Dogs were walked and socialized on a daily basis. The Beagles were also provided regular health examinations, deworming and routine vaccinations from certified veterinarians to ensure optimal health.
The test diets included one control (rice; a grain-containing diet) and five pulse-based diets [all grain-free diets: smooth pea (CDC Inca), wrinkled pea (4,140–4 variety), wrinkled pea (Amigold variety), red lentil (CDC Maxim) and fava bean (CDC Snowdrop)]. All diets were formulated at 20% available carbohydrate using locally obtained flours (flour proximate analyses shown in Supplementary Table 1). A non-digestible Celite marker was also incorporated at 1% for determination of total tract apparent digestibility and measured as acid-insoluble ash in proximate analyses of diet and feces (22). Diets were formulated using the software Creative Concept 5 (Creative Formulation Concepts, Pierz, MN, USA) and were structured to meet the nutritional requirements for canine adult maintenance (see Table 1 for formulations). These requirements were based on AAFCO and the National Research Council recommendations. Feed ingredients were sourced from local and commercial sources as needed and diets were extruded into a dry kibble format using a laboratory-scale, co-rotating, twin-screw extruder (Baker Perkins Ltd, Peterborough, UK) at the University of Winnipeg (Food Science Laboratory, Winnipeg, MB, Canada). All diets were extruded under the same conditions as described in Supplementary Table 2. Diets were then vacuum coated with fat at the Canadian Feed Research Centre (North Battleford, SK, Canada). Samples of all diets were then sent to analytical laboratories for proximate and amino acid analysis (Central Testing Laboratory Ltd., Winnipeg, MB, Canada) according to AOAC standards (23). Dry matter was determined by oven-drying the sample and crude protein determined using the Kjeldahl method while non-fiber carbohydrate and fat were determined through acid-hydrolysis solvent extraction. Metabolizable energy (ME) content of diets was determined through calculation (see footnote to Table 2 for equation). Fiber analyses were performed according to the AOAC 2011.25 method (Eurofins, Toronto, ON, Canada).
Table 1. Formulation of test diets with increasing fiber and amylose content listed from left to right. All diets were formulated to include 20% available carbohydrate.
Dogs were fed twice daily, weighed weekly and body condition scored using a 9-point scale (24). During the pre-trial phase of at least a month, each dog was fed a standard commercial diet (Purina Proplan, Mississauga, ON Canada) and individual food portion/maintenance energy required to maintain ideal weight (body condition score of 4–5 on a 9-point scale) determined. Once on trial, isocaloric test diet portions to that determined for each dog in the pre-trial phase was calculated and used throughout the feeding trial without any further adjustment to portion size.
Glycemic Index & Digestibility Testing
To establish both glycemic responses and starch digestibility of diets as well as the effects on circulating amino acid and taurine levels, Beagles were fed each test diet for 7 days. This was done in a randomized, cross-over, blinded, repeated measures design study with a 3-day washout period on the commercial diet between each test diet and another 7-day feeding period followed by 3-day washout repeated until all diets had been tested in each dog. Total tract apparent digestibility was determined in feces collected on the sixth and seventh day of each feeding period. After collection, feces were kept at−20°C until they were dried at 65°C for 72 h. Fecal samples were then sent to an analytical laboratory to assess nutrient excretion (Central Testing Laboratory Ltd., Winnipeg, MB, Canada). In addition, another portion of dried feces was used to conduct total bile acid analyses using a commercial kit according to manufacturer instructions (Total Bile Acid Assay Kit, Cell Biolabs Inc. San Diego, California, USA) which uses a colorimetric enzyme driven reaction in which bile acids are incubated in the presence of 3-alpha hydroxysteroiddehydrogenase and thio-NADH.
Total tract apparent digestibility was calculated using the formula:
On day 7 of feeding, dogs were fasted overnight and 8.0 mL of whole blood was collected the next day from the jugular vein. While still fasted, 5.0 mL of whole blood was collected into EDTA tubes and centrifuged at 2,200 RPM to obtain plasma. Samples were stored at−80°C until assayed for plasma cysteine (measured as the total of the cysteine dimer, cystine, plus the deprotonated form of cysteine called half-cystine), methionine, and taurine at a contract analytical laboratory (Amino Acid Laboratory, University of California Davis, Davis, CA, USA). Beagles were then catheterized in the cephalic vein and fasting blood glucose was determined using an Ultra2 glucometer (OneTouch, LifeScan Canada ULC, Malvern, PA, USA). Depending on the diet fed during the week, dogs were either fed glucose (oral glucose tolerance test) providing 1 g/kg body weight of glucose after consumption of the commercial diet or fed a portion of the test diet fed that week to provide 1 g/kg available carbohydrate (1 g/kg divided by % available carbohydrate) according to established protocols in dogs in our group (9). The amount of available carbohydrate in each diet was determined using a commercially obtained test kit (Available Carbohydrate Assay Kit, Megazyme, Bray, Ireland, see Table 2). Available carbohydrate was defined as total digestible starch (TDS) plus maltodextrins, sucrose, D-glucose, D-fructose and lactose. The available carbohydrate method measures glucose liberated in vitro after 4 h incubation (AOAC Method 2020.07). Blood glucose levels were monitored using a glucometer over a 5-h time period (time 0, 5, 10, 15, 20, 30, 45, 60, 90, 120, 150, 180, 210, 240, 270, 300) according to methods established in this group (10).
Sulfur Containing Amino Acid Plasma Analysis
Plasma sample analyses were conducted at a contract analytical laboratory (Amino Acid Laboratory, University of California Davis, Davis, CA, USA). Modified AOAC Official Method 994.12 alternative III was performed according to GLP (taurine and methionine recovery rates were 97–102%). Variances between duplicates were <5%. Cystine results were corrected by multiplying factor of 2 (recovery rate is about 50%).
Data Handling and Statistics
All data were tested for normality and outliers using the Kolomogorov-Smirnov test, Q-Q plots and box plots. Depending on the normality of the data either a repeated measure, 1-way ANOVA or a repeated measure, 1-way ANOVA on ranked data was then conducted followed by post-hoc Tukey's tests if significance was achieved. Differences were considered statistically significant at p ≤ 0.05. Values obtained during the oral glucose tolerance test were not used for statistical analysis and were provided strictly for reference or in glycemic index calculations. Principal components analysis was also performed on the six dietary treatment groups (34 variables studied). Factors were reduced to two components and variables studied were reduced to 24 where ≥83% of the variance across the data was explained. Analyses were performed using SigmaPlot 12.0 and Systat 12.0 (Systat Software Inc. San Jose, CA, USA).
Proximate and Amino Acid Analysis of Diets
As expected, crude fiber increased in the diets from the lowest rice control diet to the highest wrinkled pea diets (Table 2). Crude fat of the diets ranged from 16.01% to 19.99%, but was unrelated to dietary fiber (Table 2). Similarly, metabolizable energy ranged from 3,806 to 4,121 kcal/kg. In the test diets, crude protein (% dry matter) ranged from 26.13% to 33.69% which was well-above the AAFCO minimum of 18% dietary protein (25). However, dietary protein content slightly decreased with increasing fiber content, as higher levels of pulse ingredients were incorporated (Table 2). All test diets met or exceeded the AAFCO requirements of 0.33% methionine and 0.65% cystine+methionine (25), except the wrinkled pea CDC 4,140–4 diet which had 0.5% cystine+methionine (Table 2). Dietary methionine was highest in the rice diet at 0.84% on a dry matter basis, with all pulse-based diets containing lower methionine to a low of 0.38% for the high fiber CDC 4,140–4 wrinkled pea diet (Table 2). Dietary cystine ranged from a high of 1.01% for the rice diet to a low of 0.38% for the CDC 4,140–4 wrinkled pea diet, similar to the methionine results (Table 2). The grain-containing rice diet also had the highest level of taurine at 0.14% (Table 2). In contrast, the pulse containing diets had lower, varying amounts of taurine which ranged from 0.07% to 0.12% as shown in Table 3.
Fiber, Amylose and Carbohydrate Content of the Diets
Diets were formulated to contain 20% available carbohydrate and results of the available carbohydrate kit for diets showed good agreement with the target (Table 3). Available carbohydrate content (defined as the amount of glucose liberated by amyloglucosidase+pancreatin in 4 h in vitro at 37°C) in the formulated test diets ranged from 22.4% to 25.2%, as shown in Table 3. This available carbohydrate value was then used to calculate what meal size needed to be fed to provide 1 g available carbohydrate per kilogram body weight of the dogs during glycemic testing (Table 3). The same test diets prior to fat coating were tested for in vitro starch digestibility using Englyst methodology in separate study by this group (27). Uncoated wrinkled pea diets after extrusion (both Amigold and 4,140–4) had higher resistant starch content (15.1–19.0% on a dry starch basis) compared to the other pulse-based or rice-based diets (2.7–5.5% on a dry starch basis). Also, both wrinkled pea diets had lower gelatinized starch after extrusion (11.9–12.3% dry matter) compared with those of round pea, lentil, fava bean, and rice diets (17.1–19.6%; 9).
Dietary fiber analyses of the fat coated test diets used in this feeding study demonstrated that they had varying levels of dietary fiber (Table 3). Low molecular weight dietary fiber varied from <0.6–4.9% for the rice and the wrinkled pea (Amigold) diets, respectively (Table 3). High molecular weight dietary fiber was subdivided into two categories: insoluble high molecular weight dietary fiber and soluble high molecular weight dietary fiber. The soluble high molecular weight dietary fiber varied from 0.9% to 1.9% (rice and Amigold wrinkled pea diets, respectively), while insoluble high molecular weight dietary fiber varied from 15.2% to 21.8% (rice and Amigold wrinkled pea diets, respectively; Table 3). Similarly, total dietary fiber ranged from 16.1% to 28.6% (rice and Amigold wrinkled pea diets, respectively; Table 3). Finally, the wrinkled pea (Amigold) pulse-containing diet had the highest amylose content (14.6%) while the rice-containing diet had the lowest level of amylose at 5.1% (Table 3).
Total tract apparent digestibility of crude protein varied across the diets (ANOVA, p < 0.001). Crude protein digestibility varied from 72.95% to 84.35%, the rice diet had the highest crude protein digestibility, and the wrinkled pea (Amigold) diet had the lowest (Table 4). For sulfur-containing amino acids, the total tract apparent digestibility of cystine, methionine, cystine+methionine and taurine varied among test diets (ANOVA, p < 0.001 for all but taurine where p = 0.034; Table 4). Variability was noted for cystine digestibility which ranged from 6.5% to 90% (Table 4). This variability was not due to diarrhea. Although fecal output and score (quality) were not quantitated in this study, qualitatively, no obvious changes were observed among diets. All other macronutrient (fat and starch) digestibility values were also different among diets (Table 4).
Table 4. Total tract apparent digestibility analyses of protein, fat, starch, cystine, methionine, cystine+methionine and taurine in the 6 different test diets formulated at 20% available carbohydrate with variable amounts of fiber, fed for 7 days.
Blood glucose increased from fasting after feeding glucose or a meal with 1 g available carbohydrate/kg body weight, then returned to baseline within 5 h in dogs (Figure 1). From these glycemic response figures, quantitative data and statistical analyses on peak, time to peak, area under the curve and glycemic index were calculated (Table 5). In dogs fed the commercial diet for 7 days, the fasting blood glucose level prior to the oral glucose tolerance test was 3.8 ± 0.2 mmol/L. After 7 days of feeding the test diets, fasting blood glucose had variation with diet in these same dogs (Figure 1; ANOVA, p < 0.001; Table 5). After ingestion of the glucose standard, glucose levels peaked at 52.5 ± 5.5 min with a blood glucose value of 6.3 ± 0.2 mmol/L (Figure 1, Table 5). Time to peak blood glucose was longer after a meal of a whole diet and the peak was lower when compared to the response to the glucose standard (Figure 1).
Figure 1. Time course of blood glucose response in fasted dogs after oral glucose challenge (1 g glucose/kg body weight) after 7 days of feeding a commercial husbandry diet (filled circles). Time course of blood glucose response in fasted dogs after consuming a meal of each test diet (1 g available carbohydrate/kg body weight) is also shown. Dogs were fed each test diet formulated at 20% available carbohydrate with increasing levels of fiber and amylose (lowest to highest from top to bottom in figure legend) for 7 days prior to an overnight fast, followed by the test meal the next day. Data is shown as Mean ± SEM; n = 8 dogs.
Table 5. Quantitative measures of glycemic response in fasted beagles after feeding a meal of different test diets.
Dogs fed the lentil-based diet and wrinkled pea (4,140–4) diet for 7 days had the lowest fasting blood glucose levels (3.5 ± 0.1 mmol/L) and (3.6 ± 0.1 mmol/L), respectively (Table 5). The peak blood glucose was also different among diets (ANOVA, p = 0.01), while the time to peak was not different (ANOVA, p = 0.20). Dogs fed the rice diet had the highest peak blood glucose at 5.0 ± 0.09 mmol/L, while the lowest peak blood glucose was observed with the lentil diet at 4.4 ± 0.1 mmol/L (Table 5). The area under the blood glucose response curve (AUC) was different among diets (ANOVA, p = 0.02; Table 5), ranging from 810.9 ± 15.8 mmol/L x min for the rice diet to a low of 726.5 ± 21.7 mmol/L x min for the lentil diet (Table 5). Glycemic index values of the diets followed the AUC trend, with differences among diets (ANOVA, p < 0.001; Table 5). The rice diet had the highest glycemic index of 95.7 ± 2.2, while the lowest glycemic index was observed in the lentil diet at 85.8 ± 2.8 (Table 5).
Sulfur Containing Amino Acids in Plasma
After 7 days of feeding each diet, fasting plasma taurine was different in the dogs after feeding different diets (ANOVA, p = 0.021; Table 6). Feeding rice and the lentil diets produced the highest plasma taurine levels in the dogs at 99 and 111 nmol/ml (equivalent to μmol/L), while feeding the fava bean diet produced the lowest plasma taurine at 73 nmol/ml (Table 6). Despite differences in dietary content of cystine and methionine (Table 2), after 7 days of feeding each diet, fasting levels of plasma methionine and half-cystine were not significantly different (Table 6). In contrast, plasma cysteine varied greatly among the dogs fed the different diets with the highest levels when the lentil and rice diets were fed.
Table 6. Plasma amino acid levels of taurine, half-cystine, cysteine and methionine observed in fasted dogs after 7 days of feeding each test diet.
Bile Acid Assay
Fecal bile acid content was significantly different among diets, with the rice diet having the highest value and the wrinkled pea 4,140–4 diet having the lowest (Figure 2).
Figure 2. Fecal bile acid content from dogs after 7 days of feeding each test diet. Diets are listed in order of increasing fiber and amylose content from left to right. Data is shown as Mean ± SEM; n = 8 dogs. Different letters indicate significant differences using Tukey's post-hoc analysis (P < 0.05) after 1-Way Repeated Measures ANOVA.
Principal Components Analysis
In order to assess the relationships among the variables tested, a Principal Components Analysis (PCA) was run initially with all variables measured in all 6 test diets. The top 24 variables that were correlated to the first two components were retained in the final PCA analysis. These 24 variables explained 83.13% of variance in this final analysis (see Table 7). A factor loadings plot of these variables (Figure 3A) clustered variables with the greatest tendency to be positively related to each other. A plot of weighted factor scores for each of these 24 variables for each diet resulted in a graph where the rice diet was clearly separated from the two wrinkled pea diets, while fava bean, smooth pea and lentil diets were intermediate (Figure 3B). This confirms the order we predicted based on fiber and amylose content. In support of this, the greatest separation among the diets was along the x-axis showing Factor 1 scores (Figure 3B) where high positive scores indicated diets with high dietary content of chicken meal, crude protein, sulfur amino acids (cystine, methionine and taurine) as well as high fecal bile salt content, available carbohydrate and fat digestibility, but low values for dietary content of amylose and all fiber fractions (Figure 3B). The diets did not separate as much on the y-axis except for the lentil diet, but the 95% confidence intervals were very large for this diet and it did overlap with all diets on this axis (Figure 3B). High positive y-axis values for Factor 2 were associated with high fasting blood glucose, long time to peak blood glucose, high glycemic index and high plasma methionine values, but low plasma cystine, half-cysteine and taurine values.
Table 7. Principal Components Analysis of the top 24 variables examined in this study that explained 83% of the variation among diets.
Figure 3. Plots from Principal Components Analysis of the top 24 variables that explained the most variability among diets. (A) Factor loadings plot for factor 1 and factor 2. (B) A plot of weighted factor scores for components 1 vs. components 2 for each test diet in this study from the Principal Components Analysis. The size of the shape indicates the 95% confidence intervals and distance between shapes indicates greater difference among diets. Higher, positive scores for Factor 1 were associated with high dietary chicken content, dietary crude protein, dietary sulfur amino acid content (SAA; cysteine, methionine and taurine), fecal bile salts, dietary available carbohydrate and fat digestibility, but low amylose and all fiber types. In contrast, a high positive score for Factor 2 was associated with high fasting blood glucose, long time to peak blood glucose, high glycemic index and high plasma methionine, but low plasma cysteine, plasma half-cysteine and plasma taurine. SOL_HMWDF, soluble high molecular weight dietary fiber; INSOL_HMWDF, insoluble high molecular weight dietary fiber; TDF, total dietary fiber; LMWDF, low molecular weight dietary fiber; TTP_GLUCOSE, time to peak glucose; FASTING_GLUC, fasting blood glucose; AUC_GLUCOSE, area under the curve for blood glucose response; GI, glycemic index; PLASMA_MET, plasma methionine; MET_DIET, dietary level of methionine; TAUR_DIET, dietary level of taurine; DIG_FAT, total tract apparent digestibility of fat; DIG_STARCH, total tract apparent digestibility of starch; FECAL BILE, fecal bile salts; DIG_PROTEIN, total tract apparent digestibility of protein; DIG_MET, total tract apparent digestibility of methionine; AVAIL_CHO, available carbohydrate in diet; PLAS_TAU, plasma taurine; PLAS_HALFCYS, plasma half-cysteine; PLAS_CYS, plasma cystine.
The most important findings of this study were that the pulse-based, grain-free diets produced a lowered glycemic response, which could be utilized to promote increased satiety and decreased risk of diabetes mellitus. This postprandial glycemic response was negatively associated with plasma taurine and cysteine levels, but positively associated with plasma methionine levels in dogs after 7-days of feeding each test diet. Surprisingly, fiber content (all fractions) and dietary amylose content were not strongly related to plasma sulfur amino acid levels. Instead, fiber and amylose were negatively correlated to digestibility of all macronutrients, including sulfur containing amino acids. Despite increasing levels of dietary fiber as amylose content increased among diets, excretion of fecal bile acids in this study unexpectedly decreased. In addition, on a short-term feeding basis (7 days), grain-free diets did not cause a detrimental impact on the plasma taurine status in dogs, despite decreased taurine digestibility. This study showed promising effects of grain-free, pulse-based diets that could be utilized for the improvement of health in dogs.
Test Diet Properties and Relationship to Fiber Content
All test diets were formulated to be as similar as possible, while aiming to achieve ~20% available carbohydrate. This produced desired variations among the diets regarding their dietary fiber and amylose content. Diets with increasing levels of dietary fiber and amylose were among the diets that incorporated the largest amounts of pulse flour which provided plant protein and subsequently needed decreased amounts of animal protein (chicken meal) to be isonitrogenous. In addition, decreasing the amount of animal protein among diets produced decreasing levels of taurine in diets. Other studies confirmed that seafood and poultry products contained high concentrations of taurine, while plant products do not (28). Dietary fiber also did not appear to have any impacts on the other amino acids studied within the different diets such as cystine, methionine and cysteine + methionine. All pulse diets had lower levels of these amino acids in comparison to the rice diet due to pulses possessing limited content of these amino acids (6, 29). This could explain why there were variations of dietary cystine, methionine and cystine+methionine in the test diets.
In this study increasing amylose content among the different diets was associated with increasing levels of total dietary fiber. This was due to pulses with higher amylose levels containing larger amount of resistant starch, which is the fraction of starch that reaches the colon intact, as it is not digested in the small intestine (30). As reviewed by (30), resistant starch such as amylose is also viewed as dietary fiber, partially explaining the link between increased dietary fiber and amylose among diets. Increased amounts of amylose and dietary fiber influence nutrient digestibility and postprandial glucose responses (13). Thus, based on those diets we produced, our results support our hypothesis that the diets with higher fiber and lower animal-source protein exert negative effects on macronutrient and amino acid digestibility and produce lower glycemic responses. However, after 7-days of feeding diets, all dog plasma levels of cysteine, methionine and taurine were within normal range. Longer term studies are needed to confirm this.
Effect of Different Pulse-Based Diets on Glycemic Response
The diets for this study were chosen to demonstrate how increasing dietary fiber produced positive attributes on glycemic responses but might have negative impact on other nutrients in dogs. In humans, increased consumption of dietary fiber and amylose decreased postprandial blood glucose (31–34). High amylose pulse crops when used as an ingredient in dog diets minimized the risk of obesity and diabetes mellitus (13) and pulse-based diets have low glycemic properties in dogs [(9, 13, 35), Briens et al., unpublished]. However, a study systematically using a gradient of amylose content in dog diets had not yet been performed. The current study demonstrates that while amylose and total dietary fiber content is clearly an important factor in decreasing glycemic response in dogs, the spurious low results with the lentil diet point likely contributed to the lack of correlation of all the glycemic endpoints to dietary fiber and amylose content.
In addition to dietary fiber content determining glycemic response, other factors include digestion rate, amount of diet ingested, processing factors and dietary composition (13). In this study, all diets had 22–25% total digestible starch. However, since wrinkled pea flour contains much lower starch levels (~34%), the high amylose wrinkled pea diets contained up to 65% flour which contributed a large amount of fiber compared to the low amylose diets. In human studies, resistant starch is an indication of high levels of amylose and contributes to decreased starch digestibility and increased levels of total dietary fiber (36). In similar studies, changes in dietary fiber impacts glycemic response by slowing down the rate of passage of feed and the rate of hydrolysis on polysaccharides in starch (12, 37). Postprandial glucose responses are further impacted by dietary fiber as it is believed to prolong glucose absorption, thus reducing variations in glucose responses (13). Overall, glycemic response in this study was less impacted by the varying levels of amylose and dietary fiber within the different test diets and instead related to some other unidentified factor that was associated with high plasma taurine and cysteine. Further studies are needed to confirm this finding and explore how this happens.
Effect of Amylose vs. Fiber on Macronutrient and Amino Acid Digestibility
The results of the current study which found decreasing digestibility of all macronutrients (protein, fat and starch) with increasing dietary amylose and total dietary fiber agrees well-with a study conducted by (38). Digestibility of crude protein decreased as total dietary fiber consumption increased in dogs (38). This could be explained possibly by endogenous factors within the pulse flours having intrinsic interactions to form structures with starch, such as amylose, to limit the digestibility of protein (39). Similar to what is seen with protein digestion, lipids also interact with starches, creating a single-helical structure with amylose molecules and limits the enzymatic digestibility of starch (40). These amylose-lipid complexes that are resistant to starch digestion are formed during the exposure to elevated temperatures, which occur during the process of extrusion (41, 42). Due to the possibility of these processes, lipid digestibility in this study decreased as amylose levels and total dietary fiber increased. For example, a study in growing pigs demonstrated that increasing levels of fiber decreased the total tract apparent digestibility of both crude protein and fat (43).
In addition to amylose and dietary fiber impacting the digestibility of macronutrients, they also impact the digestibility of amino acids. In this study increasing levels of amylose and dietary fiber were negatively correlated to digestibility of cystine, methionine and taurine. Diets containing high amounts of dietary fiber not only lead to a greater possibility of sulfur amino acid excretion, but also greater microbial overgrowth and taurine assimilation by the microflora (18). This would be detected as increased total tract apparent digestibility, so is inconsistent with our observations. What this study could not determine is whether a particular fiber fraction or amylose content were driving changes in sulfur amino acid digestibility since all of fiber fractions and amylose were all equally negatively associated with digestibility.
Effect of Different Pulse-Based Diets on Plasma Levels of Sulfur-Containing Amino Acids
Unlike cats, dogs are able to synthesize taurine from the sulfur-containing amino acids cysteine and methionine (15). Thus, taurine is not considered an essential dietary nutrient for dogs, while methionine and methionine+cysteine have dietary minima in dogs (25). In this study, plasma taurine levels ranged from 73 to 111 nmol/mL in fasted dogs after 7 days of feeding test diets which falls within previously reported taurine reference ranges of 63–194 nmol/mL (16, 44–46). Other studies disagree on whether or not grain-free diets contribute to decreased plasma taurine levels in dogs. Some studies determined that dogs consuming grain-free diets have an increased prevalence of taurine deficiencies (16, 44), while others have noted no change or improvements in taurine status in dogs when consuming grain-free diets (47, 48). Despite lower dietary taurine content with increasing dietary fiber and amylose in the current study, plasma taurine remained within normal range. This could be due to the short-term nature of the current study (7 days per diet). Another important factor linked to a lack of consistent change in plasma taurine may be that taurine levels in target tissues such as the heart are more relevant and these tissue levels could be depleted before plasma levels of these free amino acids fall (17). Moreover, the beagle breed is not predisposed to either taurine deficiencies or dilated cardiomyopathy (19). Taken together, this study demonstrated that short-term consumption of both grain-containing and grain-free diets had no major effect on plasma taurine levels in dogs.
In this study, plasma levels of half-cystine ranged from 16 to 19 nmol/ml while plasma levels of methionine varied from 50 to 56 nmol/mL. Levels for both of these amino acids are lower than, but close to, the value reported in a dog study with 46 and 57 nmol/mL, respectively (46). Plasma cysteine values had much greater and unexpected variation from 132 to 2,732 nmol/mL among diets in the current study. Cysteine numbers are suspect for two reasons. First, cysteine is unstable after sample collection and rapidly forms disulfide bonds with itself to form the dimer called cystine or with other plasma proteins which are subsequently removed before analysis (personal communication from Amino Acid Laboratory, University of California Davis, Davis, CA, USA). Second, in this experiment, the cysteine data was more likely overestimated due to interference during HPLC analysis from L-alpha-aminoadipic acid that co-eluted with cysteine (49). Coupled with the fact that dietary cystine levels did not vary that much, it seems likely that dietary fiber has no effect on plasma cysteine, methionine or taurine levels, at least after 7 days of feeding test diets.
Effect of Dietary Fiber and Amylose on Fecal Bile Acids
Contrary to the original hypothesis, fecal excretion of total bile acids decreased as dietary fiber and amylose increased in this study in beagles after 7 days of feeding test diets. The findings of this study disagree with reports that dietary fiber can bind bile acids within the intestinal lumen, leading to increased fecal excretion of bile acids (47, 50, 51). However, the results of this study agree with another report that grain-free diets did not lead to an increased excretion of bile acids in dogs (19). Soluble dietary fiber abundant in pulses, was proposed to lower taurine availability in companion animals (18). One of the major roles of taurine in dogs is conjugation with bile acids to form the predominant bile salt, taurocholate (17, 18). While soluble fibers bind bile acids to prevent their reabsorption through the entero-hepatic circulation, thereby lowering lipids (a beneficial health effect), this effect could also deplete taurine via taurocholate loss in the feces, leading to taurine wasting (17). Alternatively, legumes high in dietary fiber can also act as prebiotics for the gut microbiota, changing the microbial population composition and overgrowth, which could enhance taurine degradation in the small intestine before it can be absorbed (45). Future studies need to explore whether feeding periods >7 days cause increased bile salt loss, whether specific classes of bile acids are affected or whether the effects on intestinal microbiota lead to taurine depletion in dogs and if so, which dietary fiber components, if any, contribute to this loss.
Study Strengths and Limitations
A strength of this study was the detailed glycemic response data for each diet. Another strength of this study was the use of multiple pulses to study if dietary fiber and amylose are impacting variables in different pulses or if there is another component that should be studied. One limitation of the study design was the sample size of eight beagles. However, we combated this limitation by using a cross-over, repeated-measures design, allowing all dogs to rotate through all diets and be studied. Another limitation of this study could be insufficient duration of feeding each diet for plasma sulfur containing amino acids to change. Ongoing studies in this group are exploring effects of longer-term feeding periods but establishing what happens in shorter time frames is also important information. For digestibility measurements, a limitation was the use of total tract apparent digestibility method which was necessary since our studies require minimally invasive, non-lethal techniques to allow our dogs to be adopted into homes once retired. True digestibility can only be assessed through either lethal sampling to remove digesta along the length of the intestinal tract or through ileal sampling that requires surgical creation of a permanently disfiguring ileal cannula.
Pulses are beneficial at producing a low glycemic response in dogs and higher amylose pulses such as the wrinkled pea (4,140–4) have superior low glycemic properties in dogs. However, even pulses such as red lentil with relatively low amylose and dietary fiber could also be processed to produce a low glycemic diet. The trade-off for beneficial low glycemic properties of high dietary fiber and high amylose pulses is decreased macronutrient and amino acid digestibility. However, this study did not find high fiber or amylose to be associated with increased fecal bile acid secretion and instead observed a decrease. In addition, due to the limitations of diet formulation in this study, high-fiber pulse diets contained less animal-source protein and higher plant-based proteins compared with the rice diet, which could be another important factor leading to decreased digestibility of certain nutrients. However, despite decreased nutrient digestibility, plasma levels of sulfur-containing amino acids, including taurine remained within normal range after 7 days of feeding test diets, suggesting at least in the short term, that the benefits may outweigh any negative nutritional effects of high fiber canine diets.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
The animal study was reviewed and approved by University of Saskatchewan Animal Research Ethics Board.
CQ: formulated and created test diets, designed and performed experiments, and analyzed data and co-wrote the paper. YR: material preparation, starch analyses, data curation, validation, and review & editing. TM: performed experiments and review & editing. YA: Conceptualization, funding acquisition, investigation, project administration, resources, and review & editing. LW: Supervised animal research, experimental design, funding acquisition, investigation, resources, project administration, review & editing, and co-wrote the paper. All authors contributed to the article and approved the submitted version.
Funding for this study was provided by the Alberta Pulse Growers, Manitoba Pulse and Soybean Growers, Ontario Bean Growers, Saskatchewan Pulse Growers, Pulse Canada and Natural Sciences and Engineering Research Council (NSERC). In-kind donation of some feed ingredients were provided by Horizon Pet Foods (Rosthern, SK Canada). Fava bean, pea, lentil and wrinkled pea flours were kindly provided by T Wakertin (Plant Sciences, University of Saskatchewan, Saskatoon, SK Canada) and by M Tulbek of Alliance Grain Traders (Saskatoon, SK Canada).
Conflict of Interest
The authors declare that this study received funding or in-kind support from the Alberta Pulse Growers, Manitoba Pulse and Soybean Growers, Ontario Bean Growers, Saskatchewan Pulse Growers, Pulse Canada, Natural Sciences and Engineering Research Council (NSERC), Alliance Grain Traders (Saskatoon, SK Canada), Horizon Pet Foods (Rosthern, SK Canada) and Dr. Tom Warkentin (University of Saskatchewan, SK). The funders were not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.
Thank you to Luciana Reis, Bright Boafo Boamah, Elise Bokshowan, Marina Subramaniam and the Animal Care Unit Staff at the University of Saskatchewan for assisting in animal care and restraint for this research project.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2021.654223/full#supplementary-material
2. Dzanis DA. The association of american feed control officials dog and cat food nutrient profiles: substantiation of nutritional adequacy of complete and balanced pet foods in the United States. J Nutr. (1994) 124(suppl. 12):2535S. doi: 10.1093/jn/124.suppl_12.2535S
3. Suarez L, Peña C, Carretón E, Juste MC, Bautista-Castaño I, Montoya-Alonso JA. Preferences of owners of overweight dogs when buying commercial pet food. J Animal Physiol Animal Nutr. (2012) 96:655–9. doi: 10.1111/j.1439-0396.2011.01193.x
4. Pirsich W, von Hardenberg LM, Theuvsen L. The pet food industry: an innovative distribution channel for marketing feed products from welfare friendly production to consumers? Int J Food System Dynamics. (2017) 8:250–61. doi: 10.18461/ijfsd.v8i3.836
5. Prantil LR, Heinze CR, Freeman LM. Comparison of carbohydrate content between grain-containing and grain-free dry cat diets and between reported and calculated carbohydrate values. J Feline Med Surg. (2018) 20:349–55. doi: 10.1177/1098612X17710842
8. Ren Y, Setia R, Warkentin TD, Ai Y. Functionality and starch digestibility of wrinkled and round pea flours of two different particle sizes. Food Chem. (2021) 336:127711. doi: 10.1016/j.foodchem.2020.127711
9. Adolphe JL, Drew MD, Huang Q, Silver TI, Weber LP. Postprandial impairment of flow-mediated dilation and elevated methylglyoxal after simple but not complex carbohydrate consumption in dogs. Nutr Res. (2012) 32:278–84. doi: 10.1016/j.nutres.2012.03.002
10. Adolphe JL, Silver TI, Childs H, Drew MD, Weber LP. Short-term obesity results in detrimental metabolic and cardiovascular changes that may not be reversed with weight loss in an obese dog model. Bri J Nutr. (2014) 112:647–56. doi: 10.1017/S0007114514001214
12. Nguyen P, Dumon H, Biourge V, Pouteau E. Glycemic and insulinemic responses after ingestion of commercial foods in healthy dogs: influence of food composition. J Nutr. (1998) 128:2654S. doi: 10.1093/jn/128.12.2654S
13. Carciofi AC, Takakura FS, De-Oliveira LD, Teshima E, Jeremias JT, Brunetto MA, et al. Effects of six carbohydrate sources on dog diet digestibility and post-prandial glucose and insulin response. J Animal Physiol Animal Nutr. (2008) 92:326–36. doi: 10.1111/j.1439-0396.2007.00794.x
14. Englyst KN, Englyst HN, Hudson GJ, Cole TJ, Cummings JH. Rapidly available glucose in foods: an in vitro measurement that reflects the glycemic response. Am J Clin Nutr. (1999) 69:448–54. doi: 10.1093/ajcn/69.3.448
15. Backus RC, Ko KS, Fascetti AJ, Kittleson MD, MacDonald KA, Maggs DJ, et al. Low plasma taurine concentration in Newfoundland dogs is associated with low plasma methionine and cyst (e) ine concentrations and low taurine synthesis. J Nutr. (2006) 136:2525–33. doi: 10.1093/jn/136.10.2525
16. Kaplan JL, Stern JA, Fascetti AJ, Larsen JA, Skolnik H, Peddle GD, et al. Taurine deficiency and dilated cardiomyopathy in golden retrievers fed commercial diets. PLoS ONE. (2018) 13:e0209112. doi: 10.1371/journal.pone.0209112
17. Sanderson SL, Gross KL, Ogburn PN, Calvert C, Jacobs G, Lowry SR, et al. Effects of dietary fat and L-carnitine on plasma and whole blood taurine concentrations and cardiac function in healthy dogs fed protein-restricted diets. Am J Vet Res. (2001) 62:1616–23. doi: 10.2460/ajvr.2001.62.1616
18. Mansilla WD, Marinangeli CP, Ekenstedt KJ, Larsen JA, Aldrich G, Columbus DA, et al. Special topic: The association between pulse ingredients and canine dilated cardiomyopathy: addressing the knowledge gaps before establishing causation. J Animal Sci. (2019) 97:983–97. doi: 10.1093/jas/sky488
19. Pezzali JG, Acuff HL, Henry W, Alexander C, Swanson KS, Aldrich CG. Effects of different carbohydrate sources on taurine status in healthy beagle dogs. J Animal Sci. (2020) 98:skaa010. doi: 10.1093/jas/skaa010
25. Association of American Feed Controls Organization AAFCO. AAFCO Methods for Substantiating Nutritional Adequacy of Dog and Cat Foods. AAFCO Model Pet Food and Specialty Pet Food Regulations PF2, 4, 7, 8, 9 and/or 10 (2013). Association of American Feed Control Officials. Official Publication. Champaign, IL.
27. Ren Y, Quilliam C, Weber LP, Warkentin TD, Tulbek MC, Ai Yl. Development of Low-Glycemic Pet Foods From Normal and High-Amylose Pulse Varieties. Cereals & Grains 20 Online Meeting by Cereals & Cereal Grains Association. (2020). Available online at: https://aaccnet.confex.com/aaccnet/2020/meetingapp.cgi/Paper/5606. (accessed March 16, 2021).
28. Spitze AR, Wong DL, Rogers QR, Fascetti AJ. Taurine concentrations in animal feed ingredients; cooking influences taurine content. J Animal Physiol Animal Nutr. 87:251–62. doi: 10.1046/j.1439-0396.2003.00434.x
30. Newberry M, Berbezy P, Belobrajdic D, Chapron S, Tabouillot P, Regina A, et al. High-amylose wheat foods: a new opportunity to meet dietary fiber targets for health. Cereal Foods World. (2018) 63:188–93. doi: 10.1094/CFW-63-5-0188
33. Behall KM, Scholfield DJ, Yuhaniak I, Canary J. Diets containing high amylose vs amylopectin starch: effects on metabolic variables in human subjects. Am J Clin Nutr. (1989) 49:337–44. doi: 10.1093/ajcn/49.2.337
35. Adolphe JL, Drew MD, Silver TI, Fouhse J, Childs H, Weber LP. Effect of an extruded pea or rice diet on postprandial insulin and cardiovascular responses in dogs. J Animal Physiol Animal Nutr. (2015) 99:767–76. doi: 10.1111/jpn.12275
36. Li H, Gidley MJ, Dhital S. High-amylose starches to bridge the “fiber gap”: development, structure, and nutritional functionality. Compr Rev Food Sci Food Safety. (2019) 18:362–79. doi: 10.1111/1541-4337.12416
38. Beloshapka AN, Cross TWL, Swanson KS. Graded dietary resistant starch concentrations on apparent total tract macronutrient digestibility and fecal fermentative end-products and microbial populations of healthy adult dogs. J Animal Sci. (2020). 99(1):1–9. doi: 10.1093/jas/skaa409
42. Lin S, Hsieh F, Huff HE. Effects of lipids and processing conditions on degree of starch gelatinization of extruded dry pet food. LWT-Food Sci Technol. (1997) 30:754–61. doi: 10.1006/fstl.1997.0271
43. Zhang W, Li D, Liu L, Zang J, Duan Q, Yang W, et al. The effects of dietary fiber level on nutrient digestibility in growing pigs. J Animal Sci Biotechnol. (2013) 4:17. doi: 10.1186/2049-1891-4-17
44. Ontiveros ES, Whelchel BD, Yu J, Kaplan JL, Sharpe AN, Fousse SL, et al. Development of plasma and whole blood taurine reference ranges and identification of dietary features associated with taurine deficiency and dilated cardiomyopathy in golden retrievers: A prospective, observational study. PLos ONE. (2020) 15:e0233206. doi: 10.1371/journal.pone.0233206
45. Tôrres CL, Backus RC, Fascetti AJ, Rogers QR. Taurine status in normal dogs fed a commercial diet associated with taurine deficiency and dilated cardiomyopathy. J Animal Physiol Animal Nutr. (2003) 87:359–72. doi: 10.1046/j.1439-0396.2003.00446.x
46. Delaney SJ, Kass PH, Rogers QR, Fascetti AJ. Plasma and whole blood taurine in normal dogs of varying size fed commercially prepared food. J Animal Physiol Animal Nutr. (2003) 87:236–44. doi: 10.1046/j.1439-0396.2003.00433.x
47. Donadelli RA, Pezzali JG, Oba PM, Swanson KS, Coon C, Varney J, et al. A commercial grain-free diet does not decrease plasma amino acids and taurine status but increases bile acid excretion when fed to labrador retrievers. Trans Animal Sci. (2020) 4:txaa141. doi: 10.1093/tas/txaa141
48. Adin D, DeFrancesco TC, Keene B, Tou S, Meurs K, Atkins C, et al. Echocardiographic phenotype of canine dilated cardiomyopathy differs based on diet type. J Vet Cardiol. (2019) 21:1–9. doi: 10.1016/j.jvc.2018.11.002
50. Garcia-Diez F, Garcia-Mediavilla V, Bayon JE, Gonzalez-Gallego J. Pectin feeding influences fecal bile acid excretion, hepatic bile acid and cholesterol synthesis and serum cholesterol in rats. J Nutr. (1996) 126:1766–71. doi: 10.1093/jn/126.7.1766
Keywords: digestibility, taurine, glycemic response, pulses, canine, fiber
Citation: Quilliam C, Ren Y, Morris T, Ai Y and Weber LP (2021) The Effects of 7 Days of Feeding Pulse-Based Diets on Digestibility, Glycemic Response and Taurine Levels in Domestic Dogs. Front. Vet. Sci. 8:654223. doi: 10.3389/fvets.2021.654223
Received: 15 January 2021; Accepted: 06 April 2021;
Published: 05 May 2021.
Edited by:Joseph Wakshlag, Cornell University, United States
Reviewed by:Sergio S. González-Muñoz, Colegio de Postgraduados (COLPOS), Mexico
Charles Gregory Aldrich, Kansas State University, United States
Copyright © 2021 Quilliam, Ren, Morris, Ai and Weber. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lynn P. Weber, email@example.com