Hepatic glycogen storage diseases are associated to microbial dysbiosis

Karina Colonetti; Bruna Bento dos Santos; Tatiéle Nalin; Carolina Fischinger Moura de Souza; Eric W. Triplett; Priscila Thiago Dobbler; Ida Vanessa Doederlein Schwartz; Luiz Fernando Wurdig Roesch

doi:10.1371/journal.pone.0214582

Abstract

Introduction

The gut microbiome has been related to several features present in Glycogen Storage Diseases (GSD) patients including obesity, inflammatory bowel disease (IBD) and liver disease.

Objectives

The primary objective of this study was to investigate associations between GSD and the gut microbiota.

Methods

Twenty-four GSD patients on treatment with uncooked cornstarch (UCCS), and 16 healthy controls had their faecal microbiota evaluated through 16S rRNA gene sequencing. Patients and controls were ≥3 years of age and not on antibiotics. Faecal pH, calprotectin, mean daily nutrient intake and current medications were recorded and correlated with gut microbiome.

Results

Patients’ group presented higher intake of UCCS, higher prevalence of IBD (n = 04/24) and obesity/overweight (n = 18/24) compared to controls (n = 0 and 06/16, respectively). Both groups differed regarding diet (in patients, the calories’ source was mainly the UCSS, and the intake of fat, calcium, sodium, and vitamins was lower than in controls), use of angiotensin-converting enzyme inhibitors (patients = 11, controls = 0; p-value = 0.001) multivitamins (patients = 22, controls = 01; p-value = 0.001), and mean faecal pH (patients = 6.23; controls = 7.41; p = 0.001). The GSD microbiome was characterized by low diversity and distinct microbial structure. The operational taxonomic unit (OTU) abundance was significantly influenced by faecal pH (r = 0.77; p = 6.8e-09), total carbohydrate (r = -0.6; p = 4.8e-05) and sugar (r = 0.057; p = 0.00013) intakes.

Conclusions

GSD patients presented intestinal dysbiosis, showing low faecal microbial diversity in comparison with healthy controls. Those findings might be due to the disease per se, and/or to the different diets, use of UCSS and of medicines, and obesity rate found in patients. Although the main driver of these differences is unknown, this study might help to understand how the nutritional management affects GSD patients.

Citation: Colonetti K, Bento dos Santos B, Nalin T, Moura de Souza CF, Triplett EW, Dobbler PT, et al. (2019) Hepatic glycogen storage diseases are associated to microbial dysbiosis. PLoS ONE 14(4): e0214582. https://doi.org/10.1371/journal.pone.0214582

Editor: Mathias Chamaillard, "INSERM", FRANCE

Received: August 23, 2018; Accepted: March 17, 2019; Published: April 2, 2019

Copyright: © 2019 Colonetti et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: Raw sequences were deposited in the Sequence Read Archive (SRA), accession SRP093885. Records are accessible at https://www.ncbi.nlm.nih.gov/sra/SRP093885. Run numbers SRR7464574 to SRR7464613.

Funding: The present study was funded by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Fundação de Amparo à Pesquisa do Estado do Rio Grande do Sul (FAPERGS)- EDITAL PRONEX FAPERGS/CNPq 12/2014, processo 16/2551-0000492-7; Fundo de Incentivo à Pesquisa e Eventos do Hospital de Clínicas de Porto Alegre (FIPE-HCPA)- n. 150218 and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Hepatic Glycogen Storage Diseases (GSD) are genetic disorders caused by deficient activity of one of the enzymes involved in the glycogenolysis pathway. The global incidence is estimated at 1 case per 20,000–43,000 live births. The most common types of GSD are GSD I, GSD III and GSD Ixα [1].

In GSD I, glucose-6-phosphate cannot be dephosphorylated to free glucose. There are two major subtypes of GSDI: Ia (~80%), caused by mutations in the G6PC gene, and GSD Ib (~20%), caused by mutations in the SLC37A4 gene. The proteins produced from G6PC (catalytic activity) and SLC37A4 (transporter) work together [2]. GSD Ia involves glycogenolysis and gluconeogenesis, and the clinical manifestations are increased weight, hepatomegaly, failure to thrive, fasting hypoglycaemia, high lactate, hyperuricemia, nephromegaly and hyperlipidaemia [3]. In addition to the features presented in GSD Ia, GSD Ib also presents with susceptibility to recurrent infections, impaired neutrophil and monocyte function, and inflammatory bowel disease (Crohn’s-like IBD) [1].

Mutations in the AGL gene cause GSD type III, in which the defective glycogen debranching enzyme blocks glycogenolysis, stopping the conversion of glycogen to glucose-1-phosphate [4]. At the same time, gluconeogenesis is enhanced to help maintain endogenous glucose production. Hepatomegaly in type III GSD generally improves with age, but affected individuals may develop chronic liver disease (cirrhosis) and liver failure later in life [5].

GSD IX is caused by the inability of phosphorylase b kinase (PHKA) to break down the glycogen in liver and/or muscle cells. Type IXα glycogenosis is an X-linked disease caused by mutations in the alpha subunit of PHKA. The signs and symptoms typically begin in early childhood, but GSD IX is usually milder than the other types [6].

The treatment for the aforementioned types of GSD involves nutritional adjustments primarily, with the periodic and frequent administration of large amounts of uncooked cornstarch (UCCS) and restriction of simple carbohydrates [7] to maintain normoglycaemia and avoid glycogen storage. Usually, higher and frequent doses of UCCS are prescribed for type Ia patients and lower doses for type IX patients. The dose is adjusted according to weight and metabolic demand [8]. GSD III and IX patients may require a hyperproteic diet with fewer restrictions for simple sugars. Sometimes additional medications may be necessary.

During the last decades, our understanding of the human being has changed. We know now that the eukaryote cells encoded by our genome are not the only component of our body. Symbiont prokaryotic cells inhabiting many cavities of our body provide metabolic functions far beyond the scope of our own physiological capabilities [9]. These cells play an important role in health and disease states [10]. The gut microbes are the most studied human associated microbial communities and consists of trillions of microbes and millions of functional genes [11]. Healthy humans present a remarkable microbial diversity but with similar functions indicating that different microbial communities are associated with a healthy microbiome [12]. The gut microbiome can be influenced by diet, lifestyle, drugs and genetics of the host [13], and has been related to several features present in GSD patients including obesity, IBD and liver disease [14]. This work aimed to investigate possible associations between GSD and the gut microbiota.

Methods

This study was a cross-sectional, observational convenience sampling study, which included 24 GSD patients (Ia = 15, Ib = 5, III = 1, IXα = 3) and 16 healthy controls. All patients were recruited from the outpatient clinics of the Medical Genetics Service at Hospital de Clínicas de Porto Alegre (MGS-HCPA), Brazil from Jan/2016 to May/2017. As inclusion criteria, the subjects (patients and controls) were ≥ 3 years old and not on antibiotics. The GSD patients also were required to: a) have a genetic diagnosis of GSD and b) be on treatment with UCCS. The healthy controls were recruited by invitation as they came to routine appointments at Santa Cecília Basic Health Unit, Porto Alegre, Brazil. All subjects received a kit and printed instructions for stool collection, storage, and transport. They were also provided with printed instructions to record three days of dietary information. Each participant collected their own frozen fecal sample and three-day dietary record and submitted them to an outpatient clinic during their next routine check-up. Upon returning to the clinic, each participant answered a brief questionnaire about personal features including weight and height, eating habits, intestinal habits, medicines of recent and/or continuous usage and lifestyle. The study protocol was approved by the Ethics Committee of Hospital de Clínicas de Porto Alegre (HCPA). All participants and/or legal guardians signed an informed consent.

As a routine, GSD patients seen at the MGS-HCPA who are on UCCS therapy also receive a multivitamin prescription. Despite optimum dietary treatment other drugs could also be prescribed, mainly for type I patients, such as allopurinol, to prevent gout and urate nephropathy; angiotensin converting enzyme inhibitors, to slow-down or prevent further deterioration of renal function; citrate, to preventing or ameliorating urolithiasis and nephrocalcinosis, in addition to correcting lactacidaemia; statins to treat hypercholesterolaemia [15]; and mainly for Ib patients, G-CSF to treat neutropenia, neutrophil dysfunction and IBD; and the intestinal anti-inflammatory mesalazine (5-amino-salicylic acid), also to treat IBD [16].

Nutritional assessment, clinical data and statistical analysis

Macro and micronutrients intake by the subjects were estimated from the three-day food records through the Nutribase software (NB16Cloud, CyberSoft, Inc., Phoenix, AZ, USA). The daily nutrient intake of each participant was the sum of the nutrients of each food item. The average of the three-day intake was used for further analysis. Multivitamin consumption and other medications were not included in the nutritional assessment but were considered as variables that potentially were modifying the gut microbial composition, so they were tested by Permutational Multivariate Analysis of Variance. Clinical data, such as IBD and other relevant conditions, were accessed from medical records. BMI-for-age and Z-scores were calculated within the World Health Organization (WHO) AnthroPlus software suite. A qualitative classification for this data followed the WHO criteria [17].

Statistical analysis among the groups was performed using PASW Statistics for Windows software (Vs18.0, 2009, SPSS Inc., Chicago, USA). Numerical variables were compared using the Mann-Whitney U test. Categorical variables were compared using X², Fisher’s exact test or Continuity Correction, when necessary (with statistical significant determined by the threshold p ≤ 0.05). Statistical analyses with the microbiome feature are described below.

Bacterial DNA extraction, 16S rRNA gene amplifications and sequencing

The bacterial DNA was isolated from 0.3 mg of frozen faecal sample with QIAamp DNA Stool Mini Kit (Qiagen, Valencia, CA, USA) (Qiagen) according to manufacturer instructions and stored at -20°C until use. The NanoVue system (GE Healthcare, Chicago, IL, USAGE Healthcare) was used to assess the quality of extractions for downstream applications. For the sequencing step, the library was prepared following the procedures described by Barboza et al. [18]. Briefly, region V4 of 16S rRNA gene was amplified with the barcoded bacterial/archaeal primers 515F and 806R [19] in a reaction containing 2U of Platinum Taq DNA High Fidelity Polymerase (Invitrogen, Carlsbad, CA, USA), 4 μL 10X High Fidelity PCR Buffer, 2 mM MgSO4, 0.2 mM dNTPs, 0.1 μM of both the 806R barcoded primer and the 515F primer, 25μg of Ultrapure BSA (Invitrogen, Carlsbad, CA, USA) and approximately 50 ng of DNA template in a final volume of 25 μL. After an initial denaturation step of 5 min at 95°C, 30 cycles of 94°C for 45 s, 56°C for 45 s and 72°C for 1 minute were performed, followed by a final extension step of 10 min at 72°C. After visualization on agarose gel 1.5%, the PCR products were purified with the Agencourt AMPure XP Reagent (Beckman Coulter, Brea, CA, USA) and the final concentration of the PCR product was quantified with the Qubit Fluorometer kit (Invitrogen, Carlsbad, CA, USA) following the manufacturer's recommendations. Finally, the reactions were combined in equimolar concentrations to create a mixture composed of 16S gene amplified fragments of each sample. This composite sample was used for library preparation with the Ion OneTouch 2 System using the Ion PGM Template OT2 400 Kit (Thermo Fisher Scientific, Waltham, MA, USA). Sequencing was performed with Ion PGM Sequencing 400 on the Ion PGM System using Ion 318 Chip v2.

16S profiling data analysis

The Fastq files exported from the Ion PGM System were analysed with the BMP Operating System (BMPOS) [20] according to the recommendations of the Brazilian Microbiome Project [21]. Briefly, an Operational Taxonomic Unit (OTU) table was built using reads truncated at 200 bp and quality filtered with a maximum expected error of 0.5. After removing singletons, the sequences were clustered into OTUs at cutoff of 97% similarity, and chimeras were checked and removed to obtain representative sequences for each microbial phylotype. Taxonomic classification was carried out in QIIME version 1.9.1 [22] based on the UCLUST method against the SILVA ribosomal RNA gene database version v132 [23] with a confidence threshold of 80%. Downstream analyses were carried out with dataset rarefied to the minimum library size [24,25] in the R environment [26] using the phyloseq package [27] and vegan package [28]. The online software Microbiome Analyst [29] was used to further detect microbial biomarkers associated with GSD patients. After Cumulative Sum Scaling (CSS) normalization [30], the dataset was analysed by the non-parametric factorial Kruskal-Wallis (KW) sum-rank test followed by Linear Discriminant Analysis [31]. To make sure the biomarkers observed were not only driven by IBD-like patients, we performed one analysis using the full dataset and another analysis excluding all four IBD-like patients and matched controls.

Faecal calprotectin assay and pH measurement

Frozen faecal samples of patients and controls were thawed and aliquoted at room temperature (20°C) to perform the pH measures and calprotectin assay. To determine the faecal pH, the samples were diluted 1:10 (w/v) in distilled water. After homogenization and incubation for 5 min at room temperature, the faecal pH was measured by an electronic pH-meter (K39-1014B, KASVI, PR, Brazil) three minutes after complete electrode immersion.

The faecal calprotectin was quantified from 100 mg of faecal sample with the RIDASCREEN Calprotectin test (R-Biopharm AG) according to the manufacturer’s instructions. Calprotectin is a calcium-/zinc-binding protein, highly stable and resistant to degradation by intestinal contents (pancreatic secretions, proteases, and bacterial degradation). It is mainly produced by neutrophils in inflammation and has been amply confirmed in intestinal inflammatory diseases [32]. Calprotectin was evaluated to verify gut inflammation across groups and its influence over the number of OTUs. Due to the small sample size of GSD III and IXα, just the subtypes Ia and Ib (groups containing >15% of total sample) were compared. Results for GSD Ia and GSD Ib patients were presented as median (Q1-Q3) and as min-max to GSD III and IXα. To test the correlation among calprotectin and OUT richness, patients who were on mesalazine were excluded from analysis.

Results

The characteristics of the patients and controls are summarized in Table 1. The nutrient intake varied significantly between groups (S1 Table); the largest variation observed was the higher total carbohydrate and calorie intakes in the GSD group due to UCCS usage. The amount of protein consumed (g) and the number of calories derived from proteins did not differ between patients and controls. However, the percentage of total caloric intake from proteins was lower in patients. Patients ingested less fat (g and Kcal/day) and had a lower percentage of fat in the diet. Regarding micronutrients, patients’ diet was poor in calcium and sodium, and in vitamins B3, H, D and E in comparison to the control group’s diet.

Download:

Table 1. Sample characterization, analysis of potential confounding variables and their effect on microbial communities.

https://doi.org/10.1371/journal.pone.0214582.t001

The intakes of macro and micronutrients were similar among all the GSD types, with some kcal variation from carbohydrate intake due the difference in UCCS consumption among groups (S2 Table).

Overall 16S rRNA sequencing results, sequence quality control and control for confounding variables

After quality filtering of the 16S rRNA reads, a total of 1,786,582 high-quality sequences longer than 200 bp were obtained. To analyse whether the number of sequences from each sample was representative of the underlying bacterial community, sequence coverage was calculated (S3 Table). An average of 44,664 sequences per sample was obtained with average sequence coverage of 0.99 at the 3% dissimilarity level. This sequencing depth was sufficient to obtain excellent representation of the microbial community in these samples.

Results for suspected confounding variables that potentially were modifying the gut microbial composition are presented at Table 1 and S1 Table. The gut microbial communities were not affected by sex, age, nor the nutritional status of the subjects tested. Faecal pH was lower in patients (6.23) than in controls (7.41), and this variable affected the presence/absence and abundance of the gut microbes, with a reduced OTU count in lower pH. Only 18% of controls (n = 3) and 41% of patients (n = 10) used antibiotics within the 6 months prior to data collection. The use of antibiotics within the 6 months prior to sampling did not affect the presence/absence of microbes (p = 0.252) nor microbial relative abundance (p = 0.179) in these samples.

Hepatic GSD is associated with an abnormal gut microbial community

The analysis of overall microbial community structure revealed significant differences between patients and controls (Fig 1). According to the PERMANOVA, the microbial community structure between patients and controls differed by the presence and absence of taxa (r² = 0.182; p = 0.003) and by their relative abundances (r² = 0.166; p = 0.001). The analysis indicated that the relative abundance of taxa contributed 16% of the variation in the microbial community between patients and controls while the presence/absence of specific taxa contributed 18% to that variation.

Download:

Fig 1.

Principal coordinates analysis (PCoA) based on Bray Curtis distance matrix (A) and Euclidean distance matrix (B) show the separation of gut microbiomes between GSD patients and controls. Each point represents a microbial community from one subject; colours indicate the treatment.

https://doi.org/10.1371/journal.pone.0214582.g001

Microbial diversity as measured by richness of OTUs and by the Shannon diversity index also differed significantly (p < 0.01) between patients and controls (Fig 2). On average, control stool samples possessed 184 OTUs while the patients had only 100 OTUs. The average Shannon diversity index was 3.49 and 2.48 in controls and patients, respectively. Together, these beta and alpha diversity analyses indicated that the GSD gut microbiome is characterized by low diversity and distinct microbial structures.

Download:

Fig 2. Alpha diversity measurements of microbial communities in the GSD patients and control groups.

Each panel represents one alpha diversity measure: Richness = total number of OTUs observed, Shannon = microbial index of diversity. Boxes span the first to third quartiles; the horizontal line within the boxes represents the median. Whiskers extending vertically from the boxes indicate variability outside the upper and lower quartiles. *** indicates a statistical difference between treatments at cutoff p ≤ 0.001.

https://doi.org/10.1371/journal.pone.0214582.g002

Defining the main taxa associated with the gut microbiota of patients and controls

Specific microbial phylotypes present within the gut community might drive the main differences observed in GSD patients. To find those microbes, biomarker screening analysis was performed at different taxonomic levels. A total of 14 phyla were detected within these samples. However, more than half of the community was dominated by only three phyla: Bateroidetes (58% in controls; 47% in patients), Firmicutes (34% in controls; 39% in patients) and Proteobacteria (5.8% in controls; 10% in patients) (Fig 3). All of the other phyla had very low relative abundances. LEfSe analysis identified three microbial phyla as biomarkers with Actinobacteria and Proteobacteria overrepresented in patients while Euryarchaeota was underrepresented. In particular, Proteobacteria presented a very high LDA score (more than 3.9 orders of magnitude), reflecting a marked increase in relative abundance in patients and consistently low abundance in controls. Firmitutes had a marginally-significant difference between patients and controls (p = 0.043 and LDA score = 4.53 but FDR = 0.07).

Download:

Fig 3. The average relative abundance of phyla found in GSD patients and healthy controls.

Phyla followed by an asterisk (*) are different, both in terms of statistics and biological consistency, between patients and controls at p and FDR ≤ 0.05: Euryarchaeota (LDA score = 1.75), Actinobacteria (LDA score = 3.06) and Proteobacteria (LDA score = 3.94). Firmicutes was marginally significantly different with p = 0.064, LDA score = 4.52 and FDR = 0.112.

https://doi.org/10.1371/journal.pone.0214582.g003

At the genus level, nineteen microbial biomarkers were different, both in terms of statistics and biological consistency, between patients and controls (Table 2). Those genera were higher in controls. In patients, those genera were in low abundance and in some cases totally absent. The lack of those microbes might be reflected in the alpha and beta diversity results as mentioned previously (Figs 1 and 2). Besides, Lactobacillus and Escherichia/Shigella were found to be dominant in patients with a very high LDA score (4.36 and 3.89, respectively), highlighting the biological importance of those microbes in GSD. To remove any biases caused by patients with IBD-like symptoms (n = 4), all IBD-like patients and their respective controls were removed from the dataset and a new biomarker analysis was performed (Table 2). Similar trends as observed within the full dataset were still present in this reduced dataset. However, the Lactobacillus genus, found previously in higher abundance in patients was not observed within the dataset without IBD-like patients. On the other hand, Escherichia/Shigella was still found to be more abundant in patients than in controls (LDA score = 3.85).

Download:

Table 2. Microbial biomarkers differentiating patients with hepatic glycogenosis diseases and healthy controls.

https://doi.org/10.1371/journal.pone.0214582.t002

Correlations between the gut microbiota, diet, faecal pH and gut inflammation

Spearman correlations were calculated between the microbiome, diet, faecal pH and calprotectin (Fig 4).

Download:

Fig 4. Correlations between the microbiota and diet, faecal pH, and gut inflammation.

https://doi.org/10.1371/journal.pone.0214582.g004

The faecal pH values varied between patients and controls (Table 1), and this was important for shaping their respective differences in gut microbiomes. Differences were determined with the Euclidian distance matrix (for presence/absence of taxa) and the Bray Curtis distance matrix (for relative microbial abundance). Faecal pH was correlated with the total number of microbial OTUs such that higher faecal pH seemed to support more OTUs.

Microbial richness correlated negatively with total carbohydrate but positively with simple carbohydrates (sugar). Calprotectin seemed to have no influence over the microbiome in terms of the number of OTUs (Fig 4). In addition, there was no correlation between this inflammatory marker and gut microbial richness.

Discussion

This is the first study about the fecal microbiota of GSD patients. In hepatic GSD, high and periodic amounts of UCCS plus dietetic restriction of fast-digestion carbohydrates are the main way to treat the genetic impairment in the glycogenolytic pathway. Our data suggest that the overload of UCCS can lead to low fecal pH by favouring some bacterial genera capable of utilizing complex carbohydrates in detriment of others. The low fecal pH, in turn, also acts as an environmental selection factor to the bacteria in the lumen. Dysbiosis has been associated with IBD and obesity. IBD includes inflammatory bowel diseases of unknown aetiology and has two main forms: ulcerative colitis and Crohn’s disease (CD). CD is a chronic disease that can affect any region in the digestive tract but is more likely to involve the small and large intestines and the perianal region [33]. Enteropathy is related to type I patients, and despite GSD Ib patients are classically described as prone to IBD-Crohn’s-like due the impaired neutrophil activity, this does not explain why patients with GSD Ia also displayed serologic markers altered for IBD, even if asymptomatic [34]. Its not clear if UCCS is the cause of obesity in GSD patients [35], but the microbiome might be associated with it. Here we discuss why the changes in microbiota could be considered as a factor of influence in the phenotype of these patients and why the UCCS usage, even though not exclusively, is an important factor that contribute to that.

Since the introduction of UCCS treatment for GSD, the focus changed from mortality to morbidity and control of long-term complications [36], such as metabolic syndrome and related symptoms [37,38]. GSD type I patients are prone to obesity, and it is suspected that UCCS contributes to the aforementioned features [35,39]. GSD I patients also are subject of heavier doses of UCCS and more restrict diet in comparison with types III and IX [35]. Regarding antibiotics, although its usage clearly drives changes in the gut microbial community, subjects who were treated with antibiotics within 6 months prior to data collection, but not during the study itself, were not affected by the previous antibiotic usage.

We found that the phyla Actinobacteria and Proteobacteria were overrepresented in patients while the Euryarchaeota was underrepresented. The microbiome of GSD patients present low diversity and was highly dominated by Escherichia/Shigella.

One possible driver of the differences in gut microbiomes between patients and controls is UCCS overload, which creates an acidic environment [34,40]. In the human body, acids are generated by regular metabolic activities and through the daily intake of food [41]. Fecal pH was lower in patients than controls and stool acidification might lead to an alteration in the relative abundances of fermenting bacteria, decreasing the conversion of unabsorbable starches to short chain fatty acids (SCFAs) [34].

SCFAs, including butyrate, are compounds made by bacteria in the gut that affect several physiologic functions and serve anti-inflammatory roles [42]. Fecal pH was associated with beta diversity and bacterial families belonging to the Clostridia class, an important producer of butyrate in the gut. Several genera of SCFA-producing bacteria—Coprococcus, Blautia, Anaerostipes, Odoribacter and Faecalibacterium—were decreased in patients. Those genera were also identified in paediatric patients with Crohn’s Disease [43]. Besides, Coprococcus and Faecalibacterium were found to have significantly low abundance in patients with nonalcoholic fatty liver disease, independently of body mass index and insulin resistance [43].

The bacterial species residing within the mucous layer of the colon may influence whether host cellular homeostasis is maintained or inflammatory mechanisms are triggered. A mutualistic relationship between the colonic microbiota, their metabolic products and the host immune system is likely involved [44]. The phylum Proteobacteria was more abundant in patients than in controls while the phylum Euryarchaeota was less abundant. Proteobacteria is a gram-negative phylum with an outer membrane mainly composed of lipopolysaccharides (LPS), which are known to sustain systemic levels of low-grade inflammation [45]. Higher levels of Proteobacteria can be considered a strong marker of dysbiosis [46]. This phylum is prevalent in patients with liver cirrhosis [47]. Several serological markers for IBD were altered in GSD-Ia patients [34], and GSD Ib patients are prone to IBD CD-like. Despite the fact that calprotectin seemed not to influence the number of OTUs gut inflammation (calprotectin >50μg/g) was verified in several patients. GSD type Ib patients have shown a concentration of calprotectin ≤50μg/g, except for one patient, who had an active IBD diagnosed in the same week. This might be due to a remission state and the use of anti-inflammatory mesalazine by these patients.

In general, dysbiosis can be categorized as a) loss of beneficial organisms, b) excessive growth of potentially harmful organisms and c) loss of overall microbial diversity. These three categories often occur at the same time [48]. Dysbiosis has been implicated in a wide range of diseases, including IBD, liver disease and obesity, that are secondary manifestations in GSD patients [49]. The reason for dysbiosis remains unclear, but the overload of UCCS contributes to those characteristics. The food intake records showed a difference in the intake of calories, mainly due to the administration of UCCS in patients, as well as a difference in microbial signature that is known to be related to obesity. It is not known whether these microbiome changes are a cause or a consequence of the pathophysiologies. However, correcting the dysbiosis can improve health in some patients [50–52]. Dysbiosis can also provide biomarkers for disease detection and management [53].

Conclusion

In this study, we reported significant alterations in the intestinal environments of GSD patients versus healthy controls. Microbiota can be affected by abiotic and biotic factors, namely pH and inflammation, and the differences in these factors between patients and controls might be linked to both genetic disease and UCCS consumption. Several bacterial taxa were different in GSD patients than in controls, and those groups are consistent with the secondary phenotypic manifestations of GSD. The microbiome patterns of these patients may reinforce immune-metabolic pathways that already are altered by genetic impairment, and may also be a factor in the differential individual response to treatment. Patients may gain health and quality of life from the restoration of gut microbial diversity that has been diminished by high UCCS intake. Future research therefore should investigate ways to manipulate the gut microbiome and clarify the possible effects of doing so.

Supporting information

S1 Table. Differences in nutrient mean daily intake between healthy controls and GSD patients and their effect on microbial communities.

*Absolute number means that the estimative of ingestion was constant for all the subjects of the group.¹Mann-Whitney U test. ²Bray-Curtis. Significant (p<0.05) events are highlighted in bold.

https://doi.org/10.1371/journal.pone.0214582.s001

(PDF)

S2 Table. Summary of the finding of the GSD patients (n = 24).

OTU: operational taxonomic unit; UCCS: uncooked cornstarch; ACE: Angiotensin-converting-enzyme inhibitor (enalapril maleate); G-CSF: G-colony stimulating factor. ¹ Numeric variables were reported as Median (Q1-Q3) for GSD Ia and Ib and as Min- Max for GSD III and IXα. Qualitative variables were reported as absolute numbers. P-value was accessed to differences between the groups Ia and Ib. ^†Calprotectin and number of OTUs for patients on and without mesalazine were reported as Min-Max.

https://doi.org/10.1371/journal.pone.0214582.s002

(PDF)

S3 Table. Overall description of the 16S rRNA sequencing results among subjects.

https://doi.org/10.1371/journal.pone.0214582.s003

(PDF)

Acknowledgments

We thank all subjects for allowing us to collect samples and information for the study. PRONEX-FAPERGS, FIPE-HCPA, CAPES, CNPq

References

1. Özen H. Glycogen storage diseases: New perspectives. World J Gastroenterol. 2007;13: 2541. pmid:17552001
- View Article
- PubMed/NCBI
- Google Scholar
2. Chou JY, Jun HS, Mansfield BC. Glycogen storage disease type I and G6Pase-β deficiency: etiology and therapy. Nat Rev Endocrinol. 2010;6: 676–688. pmid:20975743
- View Article
- PubMed/NCBI
- Google Scholar
3. Bali DS, Chen YT, Austin S, Goldstein JL, Adam MP, Ardinger HH, et al. Glycogen Storage Disease Type I. 2006 Apr 19 [Updated 2016 Aug 25]. In: Adam MP, Ardinger HH, Pagon RA, et al., editors. GeneReviews® [Internet]. Seattle (WA): University of Washington, Seattle; 1993–2019. Available from: https://www.ncbi.nlm.nih.gov/books/NBK1312/.
4. Sentner CP, Hoogeveen IJ, Weinstein DA, Santer R, Murphy E, McKiernan PJ, et al. Glycogen storage disease type III: diagnosis, genotype, management, clinical course and outcome. J Inherit Metab Dis. 2016;39: 697–704. pmid:27106217
- View Article
- PubMed/NCBI
- Google Scholar
5. Kishnani PS, Austin SL, Arn P, Bali DS, Boney A, Case LE, et al. Glycogen Storage Disease Type III diagnosis and management guidelines. Genet Med. 2010;12: 446–463. pmid:20631546
- View Article
- PubMed/NCBI
- Google Scholar
6. Herbert M, Goldstein JL, Rehder C, Austin S, Kishnani PS, Bali DS. Phosphorylase Kinase Deficiency. 2011 May 31 [Updated 2018 Nov 1]. In: Adam MP, Ardinger HH, Pagon RA, et al., editors. GeneReviews® [Internet]. Seattle (WA): University of Washington, Seattle; 1993–2019. Available from: https://www.ncbi.nlm.nih.gov/books/NBK55061/.
7. Weinstein DA, Wolfsdorf JI. Effect of continuous glucose therapy with uncooked cornstarch on the long-term clinical course of type 1a glycogen storage disease. Eur J Pediatr. 2002;161 Suppl 1: S35–39.
- View Article
- Google Scholar
8. Kishnani PS, Austin SL, Abdenur JE, Arn P, Bali DS, Boney A, et al. Diagnosis and management of glycogen storage disease type I: a practice guideline of the American College of Medical Genetics and Genomics. Genet Med. 2014;16: e1–e1. pmid:25356975
- View Article
- PubMed/NCBI
- Google Scholar
9. Sleator RD. The human superorganism—of microbes and men. Med Hypotheses. 2010;74: 214–215. pmid:19836146
- View Article
- PubMed/NCBI
- Google Scholar
10. Zhang YJ, Li S, Gan RY, Zhou T, Xu DP, Li HB. Impacts of gut bacteria on human health and diseases. Int J Mol Sci. 2015;16: 7493–519. pmid:25849657
- View Article
- PubMed/NCBI
- Google Scholar
11. Friedrich MJ. Genomes of Microbes Inhabiting the Body Offer Clues to Human Health and Disease. JAMA. 2013;309: 1447. pmid:23571560
- View Article
- PubMed/NCBI
- Google Scholar
12. Consortium HMP. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486: 207–14. pmid:22699609
- View Article
- PubMed/NCBI
- Google Scholar
13. Hall AB, Tolonen AC, Xavier RJ. Human genetic variation and the gut microbiome in disease. Nat Rev Genet. 2017;18: 690–699. pmid:28824167
- View Article
- PubMed/NCBI
- Google Scholar
14. Marchesi JR, Adams DH, Fava F, Hermes GDA, Hirschfield GM, Hold G, et al. The gut microbiota and host health: a new clinical frontier. Gut. 2016;65: 330–339. pmid:26338727
- View Article
- PubMed/NCBI
- Google Scholar
15. Rake JP, Visser G, Labrune P, Leonard JV, Ullrich K, Smit GPA. Guidelines for management of glycogen storage disease type I—European Study on Glycogen Storage Disease Type I (ESGSD I). Eur J Pediatr. 2002;161 Suppl 1: S112–119.
- View Article
- Google Scholar
16. Visser G, Rake JP, Labrune P, Leonard JV, Moses S, Ullrich K, et al. Consensus guidelines for management of glycogen storage disease type 1b - European Study on Glycogen Storage Disease Type 1. Eur J Pediatr. 2002;161 Suppl 1: S120–123.
- View Article
- Google Scholar
17. WHO MULTICENTRE GROWTH REFERENCE STUDY GROUP, Onis M. WHO Child Growth Standards based on length/height, weight and age: WHO Child Growth Standards. Acta Paediatr. 2007;95: 76–85.
- View Article
- Google Scholar
18. Barboza ADM, Pylro VS, Jacques RJS, Gubiani PI. Seasonal dynamics alter taxonomical and functional microbial profiles in Pampa biome soils under natural grasslands. 2018; 18.
- View Article
- Google Scholar
19. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6: 1621–1624. pmid:22402401
- View Article
- PubMed/NCBI
- Google Scholar
20. Pylro VS, Morais DK, de Oliveira FS, dos Santos FG, Lemos LN, Oliveira G, et al. BMPOS: a flexible and user-friendly tool sets for microbiome studies. Microb Ecol. 2016;72: 443–447. pmid:27220974
- View Article
- PubMed/NCBI
- Google Scholar
21. Pylro VS, Roesch LFW, Morais DK, Clark IM, Hirsch PR, Tótola MR. Data analysis for 16S microbial profiling from different benchtop sequencing platforms. J Microbiol Methods. 2014;107: 30–37. pmid:25193439
- View Article
- PubMed/NCBI
- Google Scholar
22. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7: 335–336. pmid:20383131
- View Article
- PubMed/NCBI
- Google Scholar
23. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2012;41: D590–D596. pmid:23193283
- View Article
- PubMed/NCBI
- Google Scholar
24. Lemos LN, Fulthorpe RR, Triplett EW, Roesch LFW. Rethinking microbial diversity analysis in the high throughput sequencing era. J Microbiol Methods. 2011;86: 42–51. pmid:21457733
- View Article
- PubMed/NCBI
- Google Scholar
25. Lemos LN, Fulthorpe RR, Roesch LF. Low sequencing efforts bias analyses of shared taxa in microbial communities. Folia Microbiol (Praha). 2012;57: 409–413.
- View Article
- Google Scholar
26. R Development Core Team. R: A Language and Environment for Statistical Computing [Internet]. 2008. Available: http://www.R-project.org
- View Article
- Google Scholar
27. McMurdie PJ, Holmes S. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE. 2013;8: e61217. pmid:23630581
- View Article
- PubMed/NCBI
- Google Scholar
28. Oksanen J, Blanchet F G, Kindt R, Legendre P, Minchin P R, O’Hara R B, et al. Vegan: community ecology package. R package vegan, vers. 2.2–1. 2015.
- View Article
- Google Scholar
29. Dhariwal A, Chong J, Habib S, King IL, Agellon LB, Xia J. MicrobiomeAnalyst: a web-based tool for comprehensive statistical, visual and meta-analysis of microbiome data. Nucleic Acids Res. 2017;
- View Article
- Google Scholar
30. Paulson JN, Stine OC, Bravo HC, Pop M. Differential abundance analysis for microbial marker-gene surveys. Nat Methods. 2013;10: 1200. pmid:24076764
- View Article
- PubMed/NCBI
- Google Scholar
31. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12: R60. pmid:21702898
- View Article
- PubMed/NCBI
- Google Scholar
32. Bjarnason I. The Use of Fecal Calprotectin in Inflammatory Bowel Disease.: 4.
- View Article
- Google Scholar
33. Matsuoka K, Kobayashi T, Ueno F, Matsui T, Hirai F, Inoue N, et al. Evidence-based clinical practice guidelines for inflammatory bowel disease. J Gastroenterol. 2018;53: 305–353. pmid:29429045
- View Article
- PubMed/NCBI
- Google Scholar
34. Lawrence NT, Chengsupanimit T, Brown LM, Weinstein DA. High Incidence of Serologic Markers of Inflammatory Bowel Disease in Asymptomatic Patients with Glycogen Storage Disease Type Ia. In: Zschocke J, Baumgartner M, Morava E, Patterson M, Rahman S, Peters V, editors. JIMD Reports, Volume 24. Berlin, Heidelberg: Springer Berlin Heidelberg; 2015. pp. 123–128. https://doi.org/10.1007/8904_2015_452 pmid:26093626
35. dos Santos BB, Nalin T, Grokoski KC, Perry IDS, Refosco LF, Vairo FP, et al. Nutritional Status and Body Composition in Patients With Hepatic Glycogen Storage Diseases Treated With Uncooked Cornstarch—A Controlled Study. J Inborn Errors Metab Screen. 2017;5: 232640981773301.
- View Article
- Google Scholar
36. Moses SW. Historical highlights and unsolved problems in glycogen storage disease type 1. Eur J Pediatr. 2002;161 Suppl 1: S2–9.
- View Article
- Google Scholar
37. Melis D, Rossi A, Pivonello R, Salerno M, Balivo F, Spadarella S, et al. Glycogen storage disease type Ia (GSDIa) but not Glycogen storage disease type Ib (GSDIb) is associated to an increased risk of metabolic syndrome: possible role of microsomal glucose 6-phosphate accumulation. Orphanet J Rare Dis. 2015;10.
- View Article
- Google Scholar
38. Derks TGJ, van Rijn M. Lipids in hepatic glycogen storage diseases: pathophysiology, monitoring of dietary management and future directions. J Inherit Metab Dis. 2015;38: 537–543. pmid:25633903
- View Article
- PubMed/NCBI
- Google Scholar
39. Jelaska BK, Ostojić SB, Berović N, Kokić V. Continuous glucose monitoring in the treatment of obesity in patients with glycogen storage disease type Ia. Endocrinol Diabetes Metab Case Rep. 2014;2014.
- View Article
- Google Scholar
40. Kotarski SF, Waniska RD, Thurn KK. Starch hydrolysis by the ruminal microflora. J Nutr. 1992;122: 178–190. pmid:1370325
- View Article
- PubMed/NCBI
- Google Scholar
41. Bagchi B. Water in Biological and Chemical Processes: From Structure and Dynamics to Function [Internet]. Cambridge: Cambridge University Press; 2013. https://doi.org/10.1017/CBO9781139583947
42. Morrison DJ, Preston T. Formation of short chain fatty acids by the gut microbiota and their impact on human metabolism. Gut Microbes. 2016;7: 189–200. pmid:26963409
- View Article
- PubMed/NCBI
- Google Scholar
43. Da Silva HE, Teterina A, Comelli EM, Taibi A, Arendt BM, Fischer SE, et al. Nonalcoholic fatty liver disease is associated with dysbiosis independent of body mass index and insulin resistance. Sci Rep. 2018;8.
- View Article
- Google Scholar
44. Carding S, Verbeke K, Vipond DT, Corfe BM, Owen LJ. Dysbiosis of the gut microbiota in disease. Microb Ecol Health Dis. 2015;26.
- View Article
- Google Scholar
45. Guo H, Diao N, Yuan R, Chen K, Geng S, Li M, et al. Subclinical-Dose Endotoxin Sustains Low-Grade Inflammation and Exacerbates Steatohepatitis in High-Fat Diet–Fed Mice. J Immunol. 2016;196: 2300–2308. pmid:26810228
- View Article
- PubMed/NCBI
- Google Scholar
46. Shin N-R, Whon TW, Bae J-W. Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends Biotechnol. 2015;33: 496–503. pmid:26210164
- View Article
- PubMed/NCBI
- Google Scholar
47. Chen Y, Yang F, Lu H, Wang B, Chen Y, Lei D, et al. Characterization of fecal microbial communities in patients with liver cirrhosis. Hepatol Baltim Md. 2011;54: 562–572.
- View Article
- Google Scholar
48. DeGruttola AK, Low D, Mizoguchi A, Mizoguchi E. Current Understanding of Dysbiosis in Disease in Human and Animal Models: Inflamm Bowel Dis. 2016;22: 1137–1150. pmid:27070911
- View Article
- PubMed/NCBI
- Google Scholar
49. Rake JP, Visser G, Labrune P, Leonard JV, Ullrich K, Smit GPA. Glycogen storage disease type I: diagnosis, management, clinical course and outcome. Results of the European Study on Glycogen Storage Disease Type I (ESGSD I). Eur J Pediatr. 2002;161 Suppl 1: S20–34.
- View Article
- Google Scholar
50. Di Luccia B, Crescenzo R, Mazzoli A, Cigliano L, Venditti P, Walser J-C, et al. Rescue of Fructose-Induced Metabolic Syndrome by Antibiotics or Faecal Transplantation in a Rat Model of Obesity. Heimesaat MM, editor. PLOS ONE. 2015;10: e0134893. pmid:26244577
- View Article
- PubMed/NCBI
- Google Scholar
51. Vrieze A, Van Nood E, Holleman F, Salojarvi J, Kootte RS, Bartelsman JFWM, et al. Transfer of intestinal microbiota from lean donors increases insulin sensitivity in individuals with metabolic syndrome. Gastroenterology. 2012;143: 913–916.e7. pmid:22728514
- View Article
- PubMed/NCBI
- Google Scholar
52. Carnero-Gregorio M, Molares-Vila A, Corbalan-Rivas A, Villaverde-Taboada C, Rodriguez-Cerdeira C. Effect of VSL#3 Probiotic in a Patient with Glycogen Storage Disease Type Ia and Irritable Bowel Disease-like Disease. Probiotics Antimicrob Proteins. 2018;
- View Article
- Google Scholar
53. Pascal V, Pozuelo M, Borruel N, Casellas F, Campos D, Santiago A, et al. A microbial signature for Crohn’s disease. Gut. 2017;66: 813–822. pmid:28179361
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Özen H. Glycogen storage diseases: New perspectives. World J Gastroenterol. 2007;13: 2541. pmid:17552001
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Chou JY, Jun HS, Mansfield BC. Glycogen storage disease type I and G6Pase-β deficiency: etiology and therapy. Nat Rev Endocrinol. 2010;6: 676–688. pmid:20975743
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Bali DS, Chen YT, Austin S, Goldstein JL, Adam MP, Ardinger HH, et al. Glycogen Storage Disease Type I. 2006 Apr 19 [Updated 2016 Aug 25]. In: Adam MP, Ardinger HH, Pagon RA, et al., editors. GeneReviews® [Internet]. Seattle (WA): University of Washington, Seattle; 1993–2019. Available from: https://www.ncbi.nlm.nih.gov/books/NBK1312/.

[ref4] 4. Sentner CP, Hoogeveen IJ, Weinstein DA, Santer R, Murphy E, McKiernan PJ, et al. Glycogen storage disease type III: diagnosis, genotype, management, clinical course and outcome. J Inherit Metab Dis. 2016;39: 697–704. pmid:27106217
View Article
PubMed/NCBI
Google Scholar

[11] View Article

[12] PubMed/NCBI

[13] Google Scholar

[ref5] 5. Kishnani PS, Austin SL, Arn P, Bali DS, Boney A, Case LE, et al. Glycogen Storage Disease Type III diagnosis and management guidelines. Genet Med. 2010;12: 446–463. pmid:20631546
View Article
PubMed/NCBI
Google Scholar

[15] View Article

[16] PubMed/NCBI

[17] Google Scholar

[ref6] 6. Herbert M, Goldstein JL, Rehder C, Austin S, Kishnani PS, Bali DS. Phosphorylase Kinase Deficiency. 2011 May 31 [Updated 2018 Nov 1]. In: Adam MP, Ardinger HH, Pagon RA, et al., editors. GeneReviews® [Internet]. Seattle (WA): University of Washington, Seattle; 1993–2019. Available from: https://www.ncbi.nlm.nih.gov/books/NBK55061/.

[ref7] 7. Weinstein DA, Wolfsdorf JI. Effect of continuous glucose therapy with uncooked cornstarch on the long-term clinical course of type 1a glycogen storage disease. Eur J Pediatr. 2002;161 Suppl 1: S35–39.
View Article
Google Scholar

[20] View Article

[21] Google Scholar

[ref8] 8. Kishnani PS, Austin SL, Abdenur JE, Arn P, Bali DS, Boney A, et al. Diagnosis and management of glycogen storage disease type I: a practice guideline of the American College of Medical Genetics and Genomics. Genet Med. 2014;16: e1–e1. pmid:25356975
View Article
PubMed/NCBI
Google Scholar

[23] View Article

[24] PubMed/NCBI

[25] Google Scholar

[ref9] 9. Sleator RD. The human superorganism—of microbes and men. Med Hypotheses. 2010;74: 214–215. pmid:19836146
View Article
PubMed/NCBI
Google Scholar

[27] View Article

[28] PubMed/NCBI

[29] Google Scholar

[ref10] 10. Zhang YJ, Li S, Gan RY, Zhou T, Xu DP, Li HB. Impacts of gut bacteria on human health and diseases. Int J Mol Sci. 2015;16: 7493–519. pmid:25849657
View Article
PubMed/NCBI
Google Scholar

[31] View Article

[32] PubMed/NCBI

[33] Google Scholar

[ref11] 11. Friedrich MJ. Genomes of Microbes Inhabiting the Body Offer Clues to Human Health and Disease. JAMA. 2013;309: 1447. pmid:23571560
View Article
PubMed/NCBI
Google Scholar

[35] View Article

[36] PubMed/NCBI

[37] Google Scholar

[ref12] 12. Consortium HMP. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486: 207–14. pmid:22699609
View Article
PubMed/NCBI
Google Scholar

[39] View Article

[40] PubMed/NCBI

[41] Google Scholar

[ref13] 13. Hall AB, Tolonen AC, Xavier RJ. Human genetic variation and the gut microbiome in disease. Nat Rev Genet. 2017;18: 690–699. pmid:28824167
View Article
PubMed/NCBI
Google Scholar

[43] View Article

[44] PubMed/NCBI

[45] Google Scholar

[ref14] 14. Marchesi JR, Adams DH, Fava F, Hermes GDA, Hirschfield GM, Hold G, et al. The gut microbiota and host health: a new clinical frontier. Gut. 2016;65: 330–339. pmid:26338727
View Article
PubMed/NCBI
Google Scholar

[47] View Article

[48] PubMed/NCBI

[49] Google Scholar

[ref15] 15. Rake JP, Visser G, Labrune P, Leonard JV, Ullrich K, Smit GPA. Guidelines for management of glycogen storage disease type I—European Study on Glycogen Storage Disease Type I (ESGSD I). Eur J Pediatr. 2002;161 Suppl 1: S112–119.
View Article
Google Scholar

[51] View Article

[52] Google Scholar

[ref16] 16. Visser G, Rake JP, Labrune P, Leonard JV, Moses S, Ullrich K, et al. Consensus guidelines for management of glycogen storage disease type 1b - European Study on Glycogen Storage Disease Type 1. Eur J Pediatr. 2002;161 Suppl 1: S120–123.
View Article
Google Scholar

[54] View Article

[55] Google Scholar

[ref17] 17. WHO MULTICENTRE GROWTH REFERENCE STUDY GROUP, Onis M. WHO Child Growth Standards based on length/height, weight and age: WHO Child Growth Standards. Acta Paediatr. 2007;95: 76–85.
View Article
Google Scholar

[57] View Article

[58] Google Scholar

[ref18] 18. Barboza ADM, Pylro VS, Jacques RJS, Gubiani PI. Seasonal dynamics alter taxonomical and functional microbial profiles in Pampa biome soils under natural grasslands. 2018; 18.
View Article
Google Scholar

[60] View Article

[61] Google Scholar

[ref19] 19. Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6: 1621–1624. pmid:22402401
View Article
PubMed/NCBI
Google Scholar

[63] View Article

[64] PubMed/NCBI

[65] Google Scholar

[ref20] 20. Pylro VS, Morais DK, de Oliveira FS, dos Santos FG, Lemos LN, Oliveira G, et al. BMPOS: a flexible and user-friendly tool sets for microbiome studies. Microb Ecol. 2016;72: 443–447. pmid:27220974
View Article
PubMed/NCBI
Google Scholar

[67] View Article

[68] PubMed/NCBI

[69] Google Scholar

[ref21] 21. Pylro VS, Roesch LFW, Morais DK, Clark IM, Hirsch PR, Tótola MR. Data analysis for 16S microbial profiling from different benchtop sequencing platforms. J Microbiol Methods. 2014;107: 30–37. pmid:25193439
View Article
PubMed/NCBI
Google Scholar

[71] View Article

[72] PubMed/NCBI

[73] Google Scholar

[ref22] 22. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7: 335–336. pmid:20383131
View Article
PubMed/NCBI
Google Scholar

[75] View Article

[76] PubMed/NCBI

[77] Google Scholar

[ref23] 23. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2012;41: D590–D596. pmid:23193283
View Article
PubMed/NCBI
Google Scholar

[79] View Article

[80] PubMed/NCBI

[81] Google Scholar

[ref24] 24. Lemos LN, Fulthorpe RR, Triplett EW, Roesch LFW. Rethinking microbial diversity analysis in the high throughput sequencing era. J Microbiol Methods. 2011;86: 42–51. pmid:21457733
View Article
PubMed/NCBI
Google Scholar

[83] View Article

[84] PubMed/NCBI

[85] Google Scholar

[ref25] 25. Lemos LN, Fulthorpe RR, Roesch LF. Low sequencing efforts bias analyses of shared taxa in microbial communities. Folia Microbiol (Praha). 2012;57: 409–413.
View Article
Google Scholar

[87] View Article

[88] Google Scholar

[ref26] 26. R Development Core Team. R: A Language and Environment for Statistical Computing [Internet]. 2008. Available: http://www.R-project.org
View Article
Google Scholar

[90] View Article

[91] Google Scholar

[ref27] 27. McMurdie PJ, Holmes S. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE. 2013;8: e61217. pmid:23630581
View Article
PubMed/NCBI
Google Scholar

[93] View Article

[94] PubMed/NCBI

[95] Google Scholar

[ref28] 28. Oksanen J, Blanchet F G, Kindt R, Legendre P, Minchin P R, O’Hara R B, et al. Vegan: community ecology package. R package vegan, vers. 2.2–1. 2015.
View Article
Google Scholar

[97] View Article

[98] Google Scholar

[ref29] 29. Dhariwal A, Chong J, Habib S, King IL, Agellon LB, Xia J. MicrobiomeAnalyst: a web-based tool for comprehensive statistical, visual and meta-analysis of microbiome data. Nucleic Acids Res. 2017;
View Article
Google Scholar

[100] View Article

[101] Google Scholar

[ref30] 30. Paulson JN, Stine OC, Bravo HC, Pop M. Differential abundance analysis for microbial marker-gene surveys. Nat Methods. 2013;10: 1200. pmid:24076764
View Article
PubMed/NCBI
Google Scholar

[103] View Article

[104] PubMed/NCBI

[105] Google Scholar

[ref31] 31. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12: R60. pmid:21702898
View Article
PubMed/NCBI
Google Scholar

[107] View Article

[108] PubMed/NCBI

[109] Google Scholar

[ref32] 32. Bjarnason I. The Use of Fecal Calprotectin in Inflammatory Bowel Disease.: 4.
View Article
Google Scholar

[111] View Article

[112] Google Scholar

[ref33] 33. Matsuoka K, Kobayashi T, Ueno F, Matsui T, Hirai F, Inoue N, et al. Evidence-based clinical practice guidelines for inflammatory bowel disease. J Gastroenterol. 2018;53: 305–353. pmid:29429045
View Article
PubMed/NCBI
Google Scholar

[114] View Article

[115] PubMed/NCBI

[116] Google Scholar

[ref34] 34. Lawrence NT, Chengsupanimit T, Brown LM, Weinstein DA. High Incidence of Serologic Markers of Inflammatory Bowel Disease in Asymptomatic Patients with Glycogen Storage Disease Type Ia. In: Zschocke J, Baumgartner M, Morava E, Patterson M, Rahman S, Peters V, editors. JIMD Reports, Volume 24. Berlin, Heidelberg: Springer Berlin Heidelberg; 2015. pp. 123–128. https://doi.org/10.1007/8904_2015_452 pmid:26093626

[ref35] 35. dos Santos BB, Nalin T, Grokoski KC, Perry IDS, Refosco LF, Vairo FP, et al. Nutritional Status and Body Composition in Patients With Hepatic Glycogen Storage Diseases Treated With Uncooked Cornstarch—A Controlled Study. J Inborn Errors Metab Screen. 2017;5: 232640981773301.
View Article
Google Scholar

[119] View Article

[120] Google Scholar

[ref36] 36. Moses SW. Historical highlights and unsolved problems in glycogen storage disease type 1. Eur J Pediatr. 2002;161 Suppl 1: S2–9.
View Article
Google Scholar

[122] View Article

[123] Google Scholar

[ref37] 37. Melis D, Rossi A, Pivonello R, Salerno M, Balivo F, Spadarella S, et al. Glycogen storage disease type Ia (GSDIa) but not Glycogen storage disease type Ib (GSDIb) is associated to an increased risk of metabolic syndrome: possible role of microsomal glucose 6-phosphate accumulation. Orphanet J Rare Dis. 2015;10.
View Article
Google Scholar

[125] View Article

[126] Google Scholar

[ref38] 38. Derks TGJ, van Rijn M. Lipids in hepatic glycogen storage diseases: pathophysiology, monitoring of dietary management and future directions. J Inherit Metab Dis. 2015;38: 537–543. pmid:25633903
View Article
PubMed/NCBI
Google Scholar

[128] View Article

[129] PubMed/NCBI

[130] Google Scholar

[ref39] 39. Jelaska BK, Ostojić SB, Berović N, Kokić V. Continuous glucose monitoring in the treatment of obesity in patients with glycogen storage disease type Ia. Endocrinol Diabetes Metab Case Rep. 2014;2014.
View Article
Google Scholar

[132] View Article

[133] Google Scholar

[ref40] 40. Kotarski SF, Waniska RD, Thurn KK. Starch hydrolysis by the ruminal microflora. J Nutr. 1992;122: 178–190. pmid:1370325
View Article
PubMed/NCBI
Google Scholar

[135] View Article

[136] PubMed/NCBI

[137] Google Scholar

[ref41] 41. Bagchi B. Water in Biological and Chemical Processes: From Structure and Dynamics to Function [Internet]. Cambridge: Cambridge University Press; 2013. https://doi.org/10.1017/CBO9781139583947

[ref42] 42. Morrison DJ, Preston T. Formation of short chain fatty acids by the gut microbiota and their impact on human metabolism. Gut Microbes. 2016;7: 189–200. pmid:26963409
View Article
PubMed/NCBI
Google Scholar

[140] View Article

[141] PubMed/NCBI

[142] Google Scholar

[ref43] 43. Da Silva HE, Teterina A, Comelli EM, Taibi A, Arendt BM, Fischer SE, et al. Nonalcoholic fatty liver disease is associated with dysbiosis independent of body mass index and insulin resistance. Sci Rep. 2018;8.
View Article
Google Scholar

[144] View Article

[145] Google Scholar

[ref44] 44. Carding S, Verbeke K, Vipond DT, Corfe BM, Owen LJ. Dysbiosis of the gut microbiota in disease. Microb Ecol Health Dis. 2015;26.
View Article
Google Scholar

[147] View Article

[148] Google Scholar

[ref45] 45. Guo H, Diao N, Yuan R, Chen K, Geng S, Li M, et al. Subclinical-Dose Endotoxin Sustains Low-Grade Inflammation and Exacerbates Steatohepatitis in High-Fat Diet–Fed Mice. J Immunol. 2016;196: 2300–2308. pmid:26810228
View Article
PubMed/NCBI
Google Scholar

[150] View Article

[151] PubMed/NCBI

[152] Google Scholar

[ref46] 46. Shin N-R, Whon TW, Bae J-W. Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends Biotechnol. 2015;33: 496–503. pmid:26210164
View Article
PubMed/NCBI
Google Scholar

[154] View Article

[155] PubMed/NCBI

[156] Google Scholar

[ref47] 47. Chen Y, Yang F, Lu H, Wang B, Chen Y, Lei D, et al. Characterization of fecal microbial communities in patients with liver cirrhosis. Hepatol Baltim Md. 2011;54: 562–572.
View Article
Google Scholar

[158] View Article

[159] Google Scholar

[ref48] 48. DeGruttola AK, Low D, Mizoguchi A, Mizoguchi E. Current Understanding of Dysbiosis in Disease in Human and Animal Models: Inflamm Bowel Dis. 2016;22: 1137–1150. pmid:27070911
View Article
PubMed/NCBI
Google Scholar

[161] View Article

[162] PubMed/NCBI

[163] Google Scholar

[ref49] 49. Rake JP, Visser G, Labrune P, Leonard JV, Ullrich K, Smit GPA. Glycogen storage disease type I: diagnosis, management, clinical course and outcome. Results of the European Study on Glycogen Storage Disease Type I (ESGSD I). Eur J Pediatr. 2002;161 Suppl 1: S20–34.
View Article
Google Scholar

[165] View Article

[166] Google Scholar

[ref50] 50. Di Luccia B, Crescenzo R, Mazzoli A, Cigliano L, Venditti P, Walser J-C, et al. Rescue of Fructose-Induced Metabolic Syndrome by Antibiotics or Faecal Transplantation in a Rat Model of Obesity. Heimesaat MM, editor. PLOS ONE. 2015;10: e0134893. pmid:26244577
View Article
PubMed/NCBI
Google Scholar

[168] View Article

[169] PubMed/NCBI

[170] Google Scholar

[ref51] 51. Vrieze A, Van Nood E, Holleman F, Salojarvi J, Kootte RS, Bartelsman JFWM, et al. Transfer of intestinal microbiota from lean donors increases insulin sensitivity in individuals with metabolic syndrome. Gastroenterology. 2012;143: 913–916.e7. pmid:22728514
View Article
PubMed/NCBI
Google Scholar

[172] View Article

[173] PubMed/NCBI

[174] Google Scholar

[ref52] 52. Carnero-Gregorio M, Molares-Vila A, Corbalan-Rivas A, Villaverde-Taboada C, Rodriguez-Cerdeira C. Effect of VSL#3 Probiotic in a Patient with Glycogen Storage Disease Type Ia and Irritable Bowel Disease-like Disease. Probiotics Antimicrob Proteins. 2018;
View Article
Google Scholar

[176] View Article

[177] Google Scholar

[ref53] 53. Pascal V, Pozuelo M, Borruel N, Casellas F, Campos D, Santiago A, et al. A microbial signature for Crohn’s disease. Gut. 2017;66: 813–822. pmid:28179361
View Article
PubMed/NCBI
Google Scholar

[179] View Article

[180] PubMed/NCBI

[181] Google Scholar

Correction

Figures

Abstract

Introduction

Objectives

Methods

Results

Conclusions

Introduction

Methods

Nutritional assessment, clinical data and statistical analysis

Bacterial DNA extraction, 16S rRNA gene amplifications and sequencing

16S profiling data analysis

Faecal calprotectin assay and pH measurement

Results

Overall 16S rRNA sequencing results, sequence quality control and control for confounding variables

Hepatic GSD is associated with an abnormal gut microbial community

Defining the main taxa associated with the gut microbiota of patients and controls

Correlations between the gut microbiota, diet, faecal pH and gut inflammation

Discussion

Conclusion

Supporting information

S1 Table. Differences in nutrient mean daily intake between healthy controls and GSD patients and their effect on microbial communities.

S2 Table. Summary of the finding of the GSD patients (n = 24).

S3 Table. Overall description of the 16S rRNA sequencing results among subjects.

Acknowledgments

References