Introduction

The incidence of type 2 diabetes (T2D) is increasing and 592 million people worldwide are projected to be affected by the disease by 2035, which makes T2D a major public-health challenge1. The pivotal feature of T2D is hyperglycaemia caused by hepatic- and peripheral insulin resistance, accompanied by progressing beta-cell dysfunction. However, etiologically and pathophysiologically T2D is a heterogeneous disorder with several genetic and environmental factors affecting the risk of developing the diabetic phenotype. Among environmental factors lifestyle-related variables such as physical inactivity and energy-dense dieting have long been recognized as important determinants of T2D development, but other potentially modifiable risk factors have recently been identified, including gut microbiota dysbiosis2,3 and dietary acid load4.

Experimental lowering of blood pH within the normal range has been shown to reduce insulin sensitivity5. A diet rich in acidogenic foods can induce whole-body acid-base imbalance in the form of low-grade metabolic acidosis6 and markers of diet-induced acidosis have been associated with insulin resistance7,8. Interestingly, in a recent study of 66,485 French middle-aged women dietary acid load predicted incident T2D over a period of 14 years independent of other known risk factors4; a result that was independently corroborated in a study of 120′053 Americans9.

While early reports of dysbiosis in T2D patients2,3 have likely been confounded by pharmaceutical therapy10, a recent study of 277 Danes reported specific functional shifts in the faecal microbiome of insulin resistant, non-diabetic individuals11. Through integrated analyses of metagenomics and metabolomics, combined with mechanistic studies in mice, Prevotella copri and Bacteroides vulgatus were identified as bacterial drivers of insulin resistance and aggravated glucose intolerance, thereby substantiating the etiological role of the gut microbiota in T2D pathogenesis. Although it has been demonstrated that dietary nutrient composition is a strong determinant of intestinal ecology12, with dietary changes being mirrored by compositional changes within days13, it remains elusive whether the microbial features associated with insulin resistance and T2D in humans are induced by specific dietary patterns, including acidogenic diets.

Alkaline drinking water has been proposed by proponents of alternative medicine as a remedy to counteract the effects of an acidogenic diet. While many claims of health promoting effects of alkaline water are scientifically unsubstantiated, changes in drinking water pH have been reported to affect both gut microbiota composition and host metabolism as evidenced by two recent studies in diabetes-prone non-obese diabetic (NOD) mice14,15. A study by Wolf and colleagues showed that neutral water (pH ≈ 7) compared to acidic water (pH ≈ 3) increased the incidence of diabetes14. Acidic water also decreased the abundance of Firmicutes and increased the abundance of Actinobacteria and Proteobacteria. Contrastingly, a study by Sofi and colleagues found that changing the drinking water from acidic (pH ≈ 3) to neutral (pH ≈ 7) decreased diabetes incidence and rate of disease progression, accompanied by a reduction in the abundance of the genus Bacteroides and some species of Prevotella (P. oris, P. melaninogenica, P. loescheii, and P. copri), whereas other species of Prevotella (P. multiformis and P. shahii) were increased along with species of Parabacteroides (P. goldsteinii and P. distasonis)15. Differences in baseline gut microbiota, along with difference in other environmental circumstances, have been proposed as an explanation for the conflicting results regarding diabetes incidence and rate of progression16,17. Importantly, both studies showed substantial differences in the gut microbiota of mice drinking acidic and neutral water, demonstrating that drinking water pH can have a profound impact on the gut microbial community. NOD mice are the most commonly used rodent model of autoimmune diabetes18. Interestingly, a recent study has shown that NOD mice develop insulin resistance, with hepatic insulin resistance being present already in the pre-diabetic phase and muscle insulin sensitivity being impaired as early as three days after onset of diabetes19. Taken together these observations suggest that the effects of changes in drinking water pH reported by Wolf et al. and Sofi et al. may in part be influenced by mechanisms affecting insulin resistance in addition to augmented autoimmune insulitis.

In the present study, we investigated possible compositional changes in the gut microbial community and glucose regulation of young healthy adults following intake of alkaline vs. neutral drinking water.

Methods

Subjects and study design

Volunteers were recruited by advertisement via online resources and by poster-advertisement at the University of Copenhagen campuses. Healthy, non-smoking men, aged 18 to 35 years, with a body mass index between 20.0 to 27.0 kg·m−2 were eligible for inclusion. Volunteers were ineligible if they had previously had abdominal surgery, had not been weight stabile for at least two months, had received antibiotic treatment within 2 months, or had a habitual consumption of caloric or artificially sweetened beverages in excess of 1 L/week or an average intake of alcohol in excess of 7 units/week. Participants in need of medical treatment during the study period were excluded.

A total of 30 volunteers were included in a non-blinded, randomized cross-over study with two 2-week intervention periods with an interposed washout period of at least 3 weeks. Participants were examined before and after each intervention period (Fig. 1). Treatment sequence was assigned at random by use of a computer algorithm. Generation of allocation sequence, enrolment of participants and assignment to interventions was conducted by THH and MTT. During intervention periods participants were instructed to quench thirst by drinking the experimental water. Participants were required to ingest a minimum of two litres distributed throughout the day, and record the intake on a daily basis. Intake of regular tap water during the intervention periods was prohibited. Intake of warm water-based beverage (e.g. tea and coffee) was allowed, but had to be based on the experimental water to the greatest extent possible. A limited amount of alcohol (not exceeding 7 drinks/week) was accepted. Participants were provided with a portable, air-tight container for the purpose of quenching thirst outside the home. In the 3-week washout period between the two intervention periods participants were instructed to return to their normal drinking habits. Compliance during intervention periods was ensured by regular e-mail and text message reminders.

Figure 1
figure 1

Flowchart.

The study was approved by the Regional Committee on Health Research Ethics for the Capital Region of Denmark (Protocol# H-15011825) and conducted according to the Helsinki Declaration. Written informed consent was obtained from all participants. The study was registered retrospectively at clinicaltrials.gov (NCT02917616; https://clinicaltrials.gov/ct2/show/NCT02917616) on September 26, 2016.

Production of experimental water

Each participant was provided with a validated commercial device (PUREPRO USA CORP, IL USA) capable of producing water with specific pH by means of ion exchange and electrolysis, to be used during each intervention period. The devices were installed on the tap in the participants’ homes, which were all located in the Copenhagen area. For each participant, the same device was used during both intervention periods. During each installation, the device was tested and duplicate pH measures were made on the setting used for the intervention. For the intervention with neutral drinking water the mean pH was 7.19 (range 7.01–7.38) and for the intervention with alkaline drinking water mean pH was 9.06 (range 8.64–9.44), corresponding to a mean difference in pH of 1.87 (range 1.55–2.33) between interventions.

Examinations

Participants were examined in the morning following a 10-hour overnight fast, before and after each intervention period. Examinations took place at the Novo Nordisk Foundation Center for Basic Metabolic Research, Unit for Physiological Studies located at Frederiksberg Hospital, Frederiksberg, Denmark.

Gastrointestinal symptoms

Gastrointestinal symptoms were assessed by a digital symptom questionnaire using a visual analogue scale, containing registration of overall gastrointestinal symptoms, bloating, abdominal pain, constipation, diarrhoea, flatulence, metallic taste, stool consistency, nausea, tiredness and bowel habits.

Anthropometrics

Participants were weighed on an electronic scale (TANITA BC-420MA, Tanita Corporation of America, IL USA) without shoes, dressed in light clothing or underwear. Height was measured to the nearest 0.5 cm without shoes using a wall-mounted stadiometer (ADE MZ10023, ADE, Germany). Waist and hip circumference was measured to the nearest cm in erect position midway between the iliac crest and the lower costal margin and at the level of the pubic symphysis, respectively. Sagittal abdominal diameter was measured in the supine position at the high point of the iliac crest during exhalation. BMI was calculated by dividing the weight (kg) by the square of the height (m) and body fat percentage was assessed with bioelectric impedance analysis (TANITA BC-420MA). Blood pressure was assessed in the supine position following a 5-minute rest as the average of 3 repeated measurements made by an automatic sphygmomanometer (A&D Medical, Japan).

Urine

Second void urine samples were collected in the fasting state and pH was measured immediately after collection using a calibrated pH meter (HI 99192, Hanna Instruments, RI USA).

Blood

Venous blood was collected in the fasting state and analysed for concentrations of plasma glucose, serum insulin, and plasma high-sensitivity C-reactive protein (hsCRP). Plasma glucose was analysed using the hexokinase method (GLUC3, Roche/Hitachi Cobas e501), serum insulin was analysed using an enzyme-linked chemiluminescent immunoassay (Insulin, Roche/Hitachi Cobas e401), and hsCRP was analysed using a particle-enhanced turbidimetric immunoassay (CRPHS, Roche/Hitachi Cobas e701). Following fasting blood sampling the participants were asked to consume 75 g of glucose dissolved in 250 mL of cold water over 5 minutes. Blood was sampled at 0, 30, 60, 90 and 120 minutes. Blood samples obtained during the OGTT were analysed for plasma glucose and serum insulin, as described.

Faeces

Faecal samples were collected before and after each intervention following standardized operating procedures, including home sampling with immediate freezing at −18 °C and transfer on dry ice for final storage at −80 °C within 48 hours.

DNA extraction, 16S rRNA library preparation and sequencing

Genomic DNA was isolated from 200 mg of faeces using the NucleoSpinSoil kit (Macherey-Nagel GmbH & Co. KG, Germany) following the manufacturer’s instruction. For the cell lysis buffer SL2 + Enhancer buffer SX were used, the subsequent vortex step was replaced with repeated bead beating. DNA yield, purity and integrity were assessed using a Qubit 2.0 fluorometer, a NanoDrop 2000 spectrometer (Thermo Fisher Scientific Inc., MA USA) and agarose gel electrophoresis. Library preparation with polymerase chain reaction (PCR) amplification was performed using 20 ng bacterial DNA, 0.2 μM of each barcoded forward and reverse primer, and HotMasterMix (5Prime) solution in a total volume of 25 μL. To target the variable region 4 (V4) of the 16S rRNA gene a forward primer 515 F (5′ AATGATACGGCGACCACCGAGATCTACAC <i5> TATGGTAATTGTGTGCCAGCMGCCGCGGTAA 3′) and a reverse primer 806 R (5′AAGCAGAAGACGGCATACGAGAT <i7> AGTCAGTCAGCCGGACTACHVGGGTWTCTAAT 3′) were used; each primer consists of the appropriate Illumina adapter, an 8-nucleotide index sequence i5 and i7, a 10-nucleotide pad sequence, a 2-nucleotide linker, and the gene-specific primer20,21. The PCR reaction was performed in a Thermocycler (Eppendorf AG, Germany), using the following parameters: 3 minutes at 94 °C, followed by 28 cycles of 20 seconds at 94 °C, 30 seconds at 55 °C and 54 seconds at 72 °C. The samples were purified with a magnetic-bead based clean-up and size selection kit (Macherey-Nagel GmbH & Co. KG, Germany). Amplicons were visualized by gel electrophoresis and quantified by a Qubit 2.0 fluorometer. A master DNA pool was generated from the purified products in equimolar ratios. The DNA was sequenced using an Illumina MiSeq platform (MiSeq Reagent Kits v2, 500 cycles), generating a total of 10,862,648 (range 33,320–424,278) paired-end reads, which were merged using FLASH22, generating contigs comprising 250 ± 6 base pairs. Expected error filtering (Emax = 0.5) in USEARCH23 was used to exclude low quality contigs. Using Quantitative Insights Into Microbial Ecology (QIIME) v.1.520 the remaining high-quality contigs (3,644,107; range 10,906-153,524) were de-multiplexed and assigned to operational taxonomic units (OTU) by a minimum 97% sequence similarity against Greengenes v.13.824 using closed reference OTU picking.

Statistical analyses

All statistical tests were performed using R v3.3.1 (www.r-project.org). No a priori sample size calculation was made. The primary outcome was difference in Shannon diversity of species level OTUs after the alkaline intervention relative to the neutral water intervention. Secondary outcomes included between- and within-intervention changes in richness, community structure and membership, and taxonomic composition. Whereas nominal P-values ≤0.05 are reported, analyses were adjusted for multiplicity ad modum Benjamini-Hochberg25 and a false discovery rate (FDR) ≤5% was observed for significance.

Intervention effect on physiological outcomes

Between-intervention effect (pH9 vs. pH7) on physiological outcomes was analysed by analysis of covariance (ANCOVA) in a linear mixed model framework with subject as a random factor. Models were fitted using restricted maximum likelihood and include a treatment × visit interaction with a level for each combination of treatment (pH7 or pH9) and visit (1 through 4), allowing a) participants to have individual baselines in the first intervention period and b) for the baseline in the second period to depend on the preceding treatment, thereby taking into account inadequate randomization and potential carry-over effects, respectively. Treatment effect was estimated by a post-hoc t-test as the average difference between pre- and post-intervention visits within a period. For outcomes measured repeatedly at each visit (i.e. glucose and insulin during OGTT) a three-way time × treatment × visit interaction was included, instead of the two-way treatment × visit interaction, and a Gaussian spatial correlation structure was applied, assuming repeated measurements to be serially correlated within a visit. The global treatment effect on plasma glucose and serum insulin during an OGTT was assessed using a χ2 test. Model assumptions were checked visually by inspection of residual plots (homoscedasticity) and normal probability plots (Gaussian distribution). Logarithmic transformation of the dependent variable was applied where appropriate. Areas under the curve (AUC) were calculated using the trapezoidal method and incremental values represent the expression above fasting values. Indices of beta-cell function and insulin resistance were calculated as described in Table S1. Difference in compliance between intervention periods among study completers was tested using a Wilcoxon signed-rank test.

Intervention effect on gut microbiota composition and diversity

Downstream analyses of 16S sequencing data were performed in R using the phyloseq package v1.16.226. Prior to alpha- and beta-diversity analyses samples were rarefied to an equal sequencing depth of 10,306 reads. For analysis of between-intervention effect at OTU level, a set of core OTUs with a mean relative abundance ≥0.05% and present in at least one sample in ≥50% of participants were selected. Similarly, genus level analyses included all genera present at least once in ≥50% of participants. Between-intervention effect on richness, alpha-diversity indices and relative abundance of individual OTUs and genera was tested using linear mixed models (as for physiological outcomes), following logarithmic transformation. When testing the treatment effect on a given taxa, a pseudo count equal to the lowest detected relative abundance of that taxa was added prior to logarithmic transformation. Participants in which a given taxon was not present in any of the samples were excluded from the analyses of that taxon in order to minimize the effect of zero-inflation. Between-intervention effect on community structure was assessed by principal coordinate (PCoA) ordination using Bray-Curtis, unweighted- and weighted UniFrac metrics, and analysed using permutational multivariate analysis of variance (PERMANOVA) as implemented in the vegan R package v2.4.0. Analyses of within-intervention effect was analysed in a mixed linear model framework contrasting pre- and post-intervention measurements with a random effect of subject and adjusted for intervention sequence. For analyses considering interventions as technical replicates (as if participants received the same intervention twice) we applied a mixed linear model contrasting pre- and post-intervention measurements adjusted for intervention sequence, with a random effect of subject within period.

Results

The study was carried out from November 2015 until April 2016. A total of 30 male participants were enrolled; however, one participant withdrew consent prior to randomization. Twenty-nine participants completed the first intervention period, and 24 completed both periods; 3 were excluded due to intercurrent illness, and 2 withdrew from the study due to personal circumstances (Fig. 1). Data from all available time points are included in the analyses.

Gut microbiota

We did not observe an effect of changing drinking water pH on overall diversity as represented by Shannon’s index which was the primary outcome. Similarly, no effect was observed for OTU based richness, estimated richness (Chao1) or Simpson’s reciprocal index (Fig. S1) when comparing the alkaline and neutral water interventions. However, in post hoc analyses comparing pre- and post-intervention samples for each intervention we observed a 15.8% (95% CI: 6.9–25.5%; P = 0.0003) increase in observed richness and a 15.4% (95% CI: 5.3–26.5%; P = 0.002) increase in estimated (Chao1) richness after the neutral-water intervention which remained significant after correction for multiple testing (Q = 0.003 and 0.008, respectively). Likewise, there was a 9.5% (95% CI: 0.25–19.6%; P = 0.04) increase in overall diversity represented by the Shannon index, although the effect was not statistically significant after correction for multiplicity. We did not observe any change in richness or overall diversity during the alkaline-water intervention.

Principal coordinate ordination of community structure and community membership as assessed by Bray-Curtis dissimilarity and UniFrac (unweighted and weighted) distances revealed highly individual responses to both the alkaline and neutral water interventions (Fig. 2). The effect of alkaline water relative to neutral water on community structure was analysed by comparing post-intervention samples from the two periods (Fig. 2G–I). Using PERMANOVA analysis of Bray-Curtis dissimilarity and UniFrac distances, we did not observe any difference in effect between the two interventions. Within the individual intervention periods there was an effect (P = 0.03) of the neutral water on community membership as assessed by unweighted UniFrac (Fig. 2D), but it did not remain statistically significant when adjusted for multiplicity (Q = 0.21). We did not observe any effect on community structure or community membership during the alkaline intervention (Fig. 2A–C).

Figure 2
figure 2

Principal coordinate plots of faecal microbiota community structure. Principal coordinate plots (axes 1 and 2) of treatment effect on community structure based on unweighted UniFrac (A,D,G), weighted UniFrac (B,E,H), and Bray-Curtis distances (C,F,I). Samples from the alkaline (A–C) and neutral water intervention (D–F), as well as post-intervention samples from both periods (GI) are depicted separately. Samples from the same participant are connected by solid lines. P-values are from permutational multivariate analysis of variance.

Baseline composition of the core microbiota (Figs S2 and S3) was dominated by members of the two major phyla Bacteroidetes (49.7%) and Firmicutes (45.4%) with Prevotellaceae (24.9%), Lachnospiraceae (20.7%), Bacteroidaceae (18.6%) and Ruminococcaceae (15.6%) as the predominant families, albeit with substantial inter-individual variation. We did not observe any significant (FDR ≤5%) differential effect of the two interventions at OTU level; however, at nominal significance 10 OTUs were less abundant and one OTU was more abundant after the alkaline intervention than after the neutral intervention (Fig. 3). Aggregated at genus level results were similar with no FDR significant effects of alkaline- vs. neutral water (Table S2).

Figure 3
figure 3

Treatment effect on OTU composition. Change (%) in geometric means of relative abundance and corresponding P-values derived from linear mixed models of treatment effect i.e. alkaline versus neutral water. Prevalence indicates number of participants in which a given OTU is present in at least one sample. Abundance indicates mean relative abundance of a given OTU in baseline samples. Taxonomy of OTUs [Greengenes ID] is given at the lowest classified rank. Only OTUs with P ≤ 0.05 are annotated.

Looking at change in relative abundance of individual OTUs within each intervention period, we observed an effect of neutral water (P ≤ 0.05) on 40 out of 201 core OTUs, six of which were significant at the pre-specified FDR of 5% (Fig. 4A): three assigned to genus Bacteroides, one assigned to species Prevotella copri, one assigned to family Ruminococcaceae, and one assigned to order Clostridiales. We also observed an effect of the alkaline water intervention on the relative abundance of 11 core OTUs; but none were significant when adjusted for multiplicity (Fig. 4B). Interestingly, the increase in relative abundance of three OTUs assigned to genus Bacteroides, one assigned to genus Blautia and one to family Ruminococcaceae reached nominal significance during both periods, indicating a potential common effect of the two interventions (Fig. S4). Furthermore, when considering paired samples from both interventions as technical replicates, we observed a significant increase (FDR <5%) in relative abundance of nine OTUs (Fig. S5) assigned to order Clostridiales (2), family Ruminococcaceae (1), genus Bacteroides (5) and species Prevotella copri (1). We did not observe any genera of which the relative abundance changed significantly within either intervention period (Table S3).

Figure 4
figure 4

Intervention effect on OTU composition. Change in relative abundance (% increase/decrease in geometric mean) of core OTUs during the neutral water (A) and alkaline (B) water interventions. Red dotted line indicates the cut-off corresponding to a Q-value of 0.05, and the blue dotted line indicates a P-values of 0.05. Prevalence indicates number of participants in which a given OTU is present in at least one sample. Abundance indicates mean relative abundance (‰) of a given OTU in baseline samples. Taxonomy of OTUs [Greengenes ID] is given at the lowest classified rank. o, order. f, family. g, genus.

Glucose regulation, low-grade inflammation and urine pH

All participants were glucose tolerant and had normal fasting glucose at the baseline examination (Table 1). A global test, considering all sampled time-points during an OGTT, did not show a significant between-intervention difference in neither plasma glucose (P = 0.13) nor serum insulin (P = 0.53) concentration after the two interventions (Fig. S6A,B). Similarly, when considering stimulated glucose and insulin only (i.e. 30 through 120 minutes) there was no significant effect of alkaline water compared to neutral (P = 0.07 and 0.39 for glucose and insulin, respectively). Considering individual time-points separately, our analyses showed that plasma glucose was 9.3% (95CI: 16.5–2.6%) higher at 30 minutes after the alkaline water intervention compared to the neutral water intervention (Fig. 5). Correspondingly, incremental AUC of glucose from 0 to 30 minutes, and from 0 to 120 minutes, were higher (9.3% (16.5–2.6%; P = 0.006) and 41.9% (97.2–2.2%; P = 0.04), respectively), whereas AUCs unadjusted for fasting glucose were not. Importantly, none of the observed effects on plasma glucose withstood correction for multiplicity at a predefined FDR of 5%. There was no detectable effect of alkaline compared to neutral water on serum insulin concentration, nor did we observe an effect on indices of beta-cell function or insulin resistance (Table 2). Similarly, there was no discernible effect of alkaline compared to neutral water on hsCRP (P = 0.89). There was no measurable effect of alkaline water compared to neutral on fasting urine pH (−0.9% (95CI: –5.1–3.5%; P = 0.69)).

Table 1 Baseline Characteristics.
Figure 5
figure 5

Treatment effect on glucose and insulin during an oral glucose tolerance test. Curves are mean (±SEM) plasma glucose (A) and serum insulin (B) sampled at 5 time-points during an oral glucose tolerance test. P-value for difference in plasma glucose at 30 min from linear mixed model adjusted for age and body-mass index.

Table 2 Effect of alkaline versus neutral water on glycaemic variables.

Gastro-intestinal symptoms

We did not observe any difference in frequency of bowel movements, nor in severity of gastro-intestinal symptoms when comparing alkaline and neutral water interventions (Fig. S7A,B). However, when comparing pre- and post-intervention records, disregarding the type of intervention, VAS scores of abdominal discomfort and degree of constipation were increased during intervention by 71.6% (95CI: 10.4–166.6%; P = 0.02) and 42.0% (95CI: 3.5–94.8%; P = 0.03), respectively (Fig. S7C); albeit the effects did not withstand correction for multiplicity of testing (Q = 0.09 in both cases).

Compliance

Daily records of water consumption demonstrated that participants on average met the required intake of 2 L per day on 12 out of 14 days (range 8–14 days) during each intervention period. There was no indication of difference in compliance during the two interventions, among those who completed both intervention periods (P = 0.57, Wilcoxon signed-rank test).

Discussion

In the present study, the objective was to investigate whether changes in drinking water pH has a short-term impact on the intestinal microbial community, as well as glucose regulation of the host. During two periods of two weeks each, participant drank a minimum of 2 L of either alkaline (pH9) or neutral (pH7) water daily. When comparing the two interventions we did not observe any significant effect on neither gut microbiota, glucose regulation, nor low-grade inflammation.

In the present study we included healthy male volunteers in an attempt to reduce potential confounding effects related to sex, demographics and health status. Consequently, we cannot rule out that a change in drinking water pH might affect the composition of the gut microbiota in the case of pre-existing dysbiosis. Furthermore, it is possible that an effect on glucose metabolism might be more pronounced in insulin resistant individuals or individuals with reduced beta-cell function.

The composition of the gut microbiota at baseline was comparable to what has previously been reported in healthy individuals, both in terms of average abundances of the major taxa, as well as inter-individual variation27. As pointed out by Wolf et al.14 and Sofi et al.15 the contrasting results observed in NOD mice are possibly ascribable to baseline variation in the gut microbiota between facilities. Similar limitations apply when studying free-living humans, in which inter-individual compositional variation at baseline and intra-individual temporal variation is presumably much larger than in laboratory mice.

A recent study by Murakami and colleagues investigated the effect of bicarbonate rich mineral water vs. tap water on glucose regulation and gut microbiota composition in 19 healthy volunteers (12 women)28. They reported a significant decrease in glycoalbumin, as well as an increase in the relative abundance of the families Christensenellaecae, Dehalobacteriaceae, Bacteroidaceae, Porphyromonadaceae, Rikenellaceae, Erysipelotrichaceae, Oxalobacteraceae and a decrease in Bifidobacteriaceae following consumption of 500 ml of bicarbonate rich water for 1 week. Of notice, the concentration of bicarbonate was 100-fold higher in the bicarbonate rich water compared to tap water, yet pH was lower (7.07 compared to 7.58 in tap water). This highlights important aspects regarding the potential effects of variation in drinking water pH. The pH as a measure of hydrogen ion activity of an aqueous solution says nothing about the buffering capacity of said solution. For instance, a solution of sodium bicarbonate is many times more efficient in buffering the acidic environment in the stomach than a solution of sodium hydroxide with the same pH. Given recent results showing an increase in families Streptococcaceae and Micrococcacaea associated with the use of proton pump inhibitors29,30, it is reasonable to assume that the potential effect of drinking water on the gut microbiota is linked to the buffering capacity. Mineral composition of the water may also have an effect on gut microbiota composition, irrespective of the acid-base properties. Divalent cations affect the stability of the outer membrane of Gram-negative bacteria31 and are involved in motility, transport, and cell differentiation processes32. In the present study, we used an electrochemical process to produce water with neutral and alkaline pH. Due to financial limitations, the mineral content of the experimental water consumed by individual participants was not assessed systematically.

As has recently been demonstrated using population-based datasets33, intervention studies aiming to identify a microbiome shift specific to a known association are likely to require considerable sample sizes. This is a first-in-man study and due to the unavailability of relevant effect size estimates no a priori sample size calculation could be made when designing the study. Thus, our findings are exploratory and it cannot be ruled out that a larger sample size would have shown a significant difference between the two interventions. For instance, post hoc analyses indicated that the present study had 19.5% statistical power to detect (paired t-test) a significant decrease in alpha diversity (Shannon’s index) at a two-tailed α = 0.05 given the observed effect size34. Consequently, a study designed to test the specific hypothesis of an effect of alkaline vs. neutral drinking water on alpha diversity would require a sample size of 180 individuals to have 80% statistical power at an α = 0.05.

Limitations of the present study include certain aspects of its design. Participants were instructed to maintain their usual diet pattern. Still, any unintended changes in diet during the study period could potentially influence both the gut microbiota and glucose regulation, obscuring any actual effect of changing drinking water pH. Participants and researchers were non-blinded to intervention sequence which could potentially result in biased estimates of treatment effects, especially of self-reported, subjective outcomes like gastro-intestinal symptom severity. We sought to ensure compliance by regular text message reminders and self-reporting of consumed water amounts, but no objective assessment of adherence to the interventions was made.

In post hoc analyses of changes in microbiota composition within each intervention arm we observed a significant increase in the relative abundance of six OTUs, as well as a significant increase in observed and estimated richness, during the neutral water intervention. No significant effects were observed during the alkaline interventions. Of the OTUs that were significant during the neutral water intervention five were among the OTUs with P-values below the 10th percentile during the alkaline intervention, two of which were nominally significant. Consequently, when analysing samples from the two interventions as technical replicates we observed nine OTUs that became more abundant when changing from regular tap water to experimental water. A possible explanation for this apparent effect could be that participants drank more during intervention periods than they usually did, because a minimum daily intake of 2 litres was required. Another possible explanation is that regular beverages were substituted by water during intervention periods. Seventy six percent of participants in the present study were habitual coffee drinkers. The effects of coffee (decaffeinated as well as caffeinated) on gastro-intestinal physiology include increased gastric acid secretion, gall-bladder contraction and rectosigmoid motility35. Parallel to these stimulatory effects coffee consumption has been shown to alter the gut microbiota composition in high-fat-fed rats36 and increase the abundance of Bifidobacterium spp. in humans37. Change in coffee consumption during interventions was not recorded and any effect of beverage substitution is speculative only. It is also possible that the observed compositional change is caused by a machine effect irrespective of the chosen pH; however, the exact nature of such an effect is unknown.

Conclusions

The present study did not show any differential effect of alkaline vs. neutral drinking water on the gut microbiota composition or glucose regulation and inflammatory state of healthy, male volunteers. However, larger studies are required in order to definitively rule out any effect of variation in drinking water pH on the gut microbiota or physiology of the host. Analyses of within intervention changes in microbiota composition indicated a possible effect of quantitative or qualitative changes in habitual drinking habits which requires further experimental elucidation.

Ethics approval and informed consent

The study was approved by the Regional Committee on Health Research Ethics for the Capital Region of Denmark (Protocol# H-15011825) and conducted according to the Helsinki Declaration. Written informed consent was obtained from all participants. The study was registered at www.clinicaltrials.gov (NCT02917616) on September 28, 2016.