장 건강, 환경과 장내 유익균(내부발효) 대사 2024 네이처!!

[문형철]

1. nature  
2. nature microbiology  
3. articles  
4. article

Gut physiology and environment explain variations in human gut microbiome composition and metabolism

Download PDF

• Article
• Open access
• Published: 27 November 2024

Gut physiology and environment explain variations in human gut microbiome composition and metabolism

• Nicola Procházková, 
• Martin F. Laursen, 
• Giorgia La Barbera, 
• Eirini Tsekitsidi, 
• Malte S. Jørgensen, 
• Morten A. Rasmussen, 
• Jeroen Raes, 
• Tine R. Licht, 
• Lars O. Dragsted & 
• Henrik M. Roager 

Nature Microbiology volume 9, pages3210–3225 (2024)Cite this article

• 30k Accesses

• 30 Citations

• 379 Altmetric

• Metricsdetails

Abstract

The human gut microbiome is highly personal. However, the contribution of gut physiology and environment to variations in the gut microbiome remains understudied. Here we performed an observational trial using multi-omics to profile microbiome composition and metabolism in 61 healthy adults for 9 consecutive days. We assessed day-to-day changes in gut environmental factors and measured whole-gut and segmental intestinal transit time and pH using a wireless motility capsule in a subset of 50 individuals. We observed substantial daily fluctuations, with intra-individual variations in gut microbiome and metabolism associated with changes in stool moisture and faecal pH, and inter-individual variations accounted for by whole-gut and segmental transit times and pH. Metabolites derived from microbial carbohydrate fermentation correlated negatively with the gut passage time and pH, while proteolytic metabolites and breath methane showed a positive correlation. Finally, we identified associations between segmental transit time/pH and coffee-, diet-, host- and microbial-derived metabolites. Our work suggests that gut physiology and environment are key to understanding the individuality of the human gut microbial composition and metabolism.

초록

인간 장내 미생물 군집은

매우 개인적입니다.

그러나

장내 미생물 군집의 변이에 대한 장 생리학 및 환경의 기여도는

여전히 충분히 연구되지 않았습니다.

본 연구에서는 61명의 건강한 성인을 대상으로

9일 연속 관찰 시험을 수행하여

다중 오믹스 기술을 활용해 미생물 군집 구성과 대사 프로파일을 분석했습니다.

50명의 하위 집단에서 무선 운동성 캡슐을 사용하여

장 환경 요인의 일일 변화를 평가하고,

전체 장 및 분절별 장 통과 시간과 pH를 측정하였다.

상당한 일일 변동성을 관찰했으며,

장내 미생물군집과 대사 활동의 개인 내 변동은

대변 수분 함량 및 분변 pH 변화와 연관되었고,

개인 간 변동은 전체 장 및 분절별 통과 시간과 pH에 의해 설명되었다.

미생물 탄수화물 발효에서 유래한 대사산물은

장 통과 시간 및 pH와 음의 상관관계를 보인 반면,

단백질 분해 대사산물과 호흡 메탄은 양의 상관관계를 나타냈다.

마지막으로,

분절별 통과 시간/pH와 커피, 식이, 숙주 및 미생물 유래 대사산물 간의 연관성을 확인했다.

본 연구는

장 생리학 및 환경이

인간 장내 미생물 구성과 대사의 개체차를 이해하는 핵심 요소임을 시사한다.

Similar content being viewed by others

Intestinal pH: a major driver of human gut microbiota composition and metabolism

Article 02 July 2025

초록 인간의 위장관에서 pH는 장내 미생물 구성과 활동을 형성하는 핵심 요소인 동시에 미생물 대사에 의해 영향을 받는다. pH는 숙주, 식이, 미생물 및 외부 요인의 복합적 작용으로 인해 개인 내 위장관 구간별로, 그리고 개인 간에 상당한 차이를 보인다. 지난 세기 동안 미생물군집 구성과 대사 반응에 대한 pH의 중요성은 광범위하게 연구되어 왔다. 본고에서는 위장관 내 pH에 영향을 미치는 주요 생리적 및 식이적 요인을 탐구하기 위해 문헌을 검토한다. 미생물 생태학적 관점에서 위장관 pH가 미생물군집 구성과 대사에 미치는 영향을 논의한다. pH가 세균의 산 반응 시스템, 유전자 발현 및 건강에 중요한 미생물 대사산물 생성에 영향을 미치는 메커니즘을 탐구한다. 마지막으로, 위장관 pH가 인간 질환에서 수행할 수 있는 잠재적 역할에 관한 문헌을 검토한다. 위장관 pH를 고려함으로써 건강 및 질병 상태에서의 장내 미생물군집에 대한 이해를 심화할 수 있다고 제안한다. pH에 의해 매개되는 장내 미생물 대사 변이가 인간 건강과 관련된 대사 산출물을 예측하고 조절하는 데 매우 중요하다고 주장한다.

Temporal variability in quantitative human gut microbiome profiles and implications for clinical research

Article Open access18 November 2021

Metaproteomic portrait of the healthy human gut microbiota

Article Open access28 June 2024

Main

Diet influences the gut microbial composition and metabolism1,2. However, even with identical dietary intake, the gut microbiome varies3,4, suggesting that other factors in the gut contribute to these variations. Gut transit time accounts for substantial variation in the microbiome composition of healthy populations5,6,7,8, with longer transit time associated with increased microbial protein degradation and methane production9. While short-chain fatty acids (SCFAs), the main microbial products of saccharolysis, are typically considered beneficial10, microbial proteolysis results in metabolites associated with poor health outcomes, including hydrogen sulfide, ammonia, branched-chain fatty acids (BCFAs), p-cresol, indole and phenylacetate11,12.

Changes in pH along the gut are also linked to gut microbial composition and metabolism13. The presence of SCFAs and other organic acids lowers colonic pH13, inhibiting bacteria sensitive to acidic environments14. Yet, little is known about how the gut environment, determined by physiological factors such as transit time and luminal pH, associates with diet–host–microbiota metabolism. Understanding these factors could be crucial for future personalized dietary microbiome-based strategies.

Transit time can be assessed using wireless motility capsules (SmartPills), which measure intraluminal pH, temperature and pressure along the gastrointestinal tract that allow for the determination of segmental transit times. SmartPills are one of the standard clinical methods for directly assessing transit time, along with radio-opaque markers and scintigraphy15, and have been used in previous microbiome studies16,17.

In this Article, we conducted a 9 day observational study with 61 healthy volunteers and assessed whole-gut and segmental transit time and pH by SmartPills. We also collected data on bowel habits, dietary intake, and breath hydrogen and methane levels and profiled urine and faecal samples using multi-omics techniques. This allowed us to follow and link inter-individual and day-to-day changes in the gut environment, gut microbiota and microbiota-derived metabolites.

메인

식이요법은

장내 미생물군집 구성과 대사에 영향을 미친다1,2.

그러나

동일한 식이 섭취량에서도

장내 미생물군집은 다양하게 나타난다3,4.

이는 장내의 다른 요인들이 이러한 변이에 기여함을 시사한다.

장 통과 시간은

건강한 인구의 미생물군집 구성에 상당한 변이를 설명한다5,6,7,8.

통과 시간이 길수록

미생물 단백질 분해와 메탄 생산이 증가한다9.

당분해의 주요 미생물 생성물인 단쇄지방산(SCFAs)은

일반적으로 유익한 것으로 간주되지만10,

미생물 단백질 분해는

황화수소, 암모니아, 분지쇄 지방산(BCFAs), p-크레졸, 인돌 및 페닐아세테이트11,12와 같은

건강 악화와 연관된 대사 산물을 생성한다.

장 내 pH 변화 역시

장내 미생물 구성과 대사 활동과 연관되어 있다13.

SCFA 및 기타 유기산의 존재는 결장 pH를 낮추어13

산성 환경에 민감한 세균의 증식을 억제한다14.

아미노기 : 알칼리성
황산, 인산, 카르복실기 : 산성

그러나

통과 시간 및 관강 pH와 같은 생리적 요인에 의해 결정되는 장 환경이

식이-숙주-미생물군집 대사와 어떻게 연관되는지에 대해서는 거의 알려져 있지 않다.

이러한 요인들을 이해하는 것은

향후 맞춤형 식이 미생물군집 기반 전략 수립에 핵심적일 수 있다.

이동 시간은

무선 운동성 캡슐(SmartPills)을 사용하여 평가할 수 있으며,

이는 위장관을 따라 관강 내 pH, 온도 및 압력을 측정하여 분절별 이동 시간을 결정할 수 있게 합니다.

SmartPills는 방사선 불투과성 마커 및 신티그래피15와 함께

이동 시간을 직접 평가하는 표준 임상 방법 중 하나이며,

이전 미생물군 연구16,17에서도 사용되었습니다.

본 논문에서는

61명의 건강한 자원자를 대상으로 9일간의 관찰 연구를 수행하여

스마트필을 통해 전장 및 분절별 통과 시간과 pH를 평가했습니다.

또한 배변 습관, 식이 섭취량, 호흡 수소 및 메탄 농도 데이터를 수집하고

다중 오믹스 기법을 활용하여 소변 및 대변 샘플을 분석했습니다.

이를 통해 개인 간 및 일일 변동성을 보이는

장 환경, 장내 미생물군집 및 미생물군집 유래 대사물의 변화를 추적하고

연관성을 규명할 수 있었습니다.

Results

Study design and participant characteristics

We enrolled 61 healthy participants (aged 39 ± 13.5 years, with body mass index (BMI) of 23.6 ± 2.8 kg m−2; Table 1 and Extended Data Fig. 1) and asked them to maintain their habitual lifestyle and diet for 9 consecutive days (Fig. 1a). The study included two visits (day 2 and day 9) where fasting blood glucose, insulin and C-peptide, as well as breath hydrogen and methane, were measured (Table 1). On the first visit, participants were given a breakfast that accounted for 25% of their daily energy needs (Supplementary Table 1) to provide a standardized meal before a subset of the volunteers (n = 50) ingested a wireless motility capsule (SmartPill) to measure whole-gut and segmental transit time and pH18. While previous investigations used granola bars (SmartBar) before the monitoring19,20, we used a complex meal similar to a recent study21 to investigate diet–microbiota interactions. Postprandial breath and urine samples were obtained as indicated in Fig. 1a. The participants recorded daily 24 h dietary records (days 1–8) using the myfood24 nutrition platform (https://www.myfood24.org); noted daily bowel habits including defecation time, stool consistency assessed by the Bristol Stool Form Scale (BSS)22 and stool frequency (number of bowel movements per day); and collected daily urine (the first morning sample) and faecal samples (the first bowel movement). The study population had normal bowel habits (Table 1) with a median BSS of type 4 and 1 bowel movement per day. Transit time was also estimated by a self-administered sweet-corn transit time test23 on days 3 and 5 (corn TT). We measured faecal water content (indication of stool moisture, a proxy marker of transit time23), pH and microbial load in all collected faecal samples (n = 484). All collected urine samples (daily spot and postprandial samples, n = 1,154) and a subset of faecal samples (n = 170) were profiled by untargeted liquid chromatography–mass spectrometry (LC-MS)-metabolomics to obtain urine and faecal metabolomes. Finally, we obtained the gut microbiome composition via 16S ribosomal RNA (rRNA) gene sequencing of a subset of faecal samples (n = 362) and assessed both relative microbiome profiles (RMPs) and quantitative microbiome profiles (QMPs) after adjusting for microbial load24.

연구 설계 및 참가자 특성

61명의 건강한 참가자(평균 연령 39±13.5세,

체질량지수(BMI) 23.6±2.8 kg/m−2; 표 1 및 확장 데이터 그림 1)를 모집하여

9일 연속으로 평소 생활습관과 식단을 유지하도록 요청했습니다(그림 1a).

이 연구에는 두 차례의 방문(2일차 및 9일차)이 포함되었으며,

이때 공복 혈당, 인슐린 및 C-펩타이드와 호흡 수소 및 메탄을 측정했습니다(표 1).

첫 방문에서 참가자들은 일일 에너지 요구량의 25%에 해당하는 아침 식사(보충 자료 표 1)를 제공받았으며,

일부 자원자(n = 50)는 무선 운동성 캡슐(SmartPill)을 섭취하여

전장 및 분절별 통과 시간과 pH를 측정하기 전에 표준화된 식사를 제공받았습니다18.

이전 연구에서는 모니터링 전에 그래놀라 바(SmartBar)를 사용했지만19,20, 우리는 식이-미생물군 상호작용을 조사하기 위해 최근 연구21과 유사한 복합 식사를 사용했습니다. 식후 호흡 및 소변 샘플은 그림 1a에 표시된 대로 채취하였다.

참가자들은 myfood24 영양 플랫폼(https://www.myfood24.org)을 사용하여 1~8일차 동안 매일 24시간 식이 기록을 작성하였으며, 배변 시간, 브리스톨 변형 척도(BSS)22로 평가한 변의 농도, 배변 빈도(일일 배변 횟수)를 포함한 일일 배변 습관을 기록하였다. (일일 배변 횟수)를 기록했으며, 매일 아침 첫 소변과 첫 배변 시 대변 샘플을 채취했다.

연구 대상자는 정상적인 배변 습관(표 1)을 보였으며, 중앙값 BSS는 4형, 일일 배변 횟수는 1회였다. 또한 3일차와 5일차에 자가 투여형 스위트콘 통과 시간 검사23(옥수수 TT)을 통해 통과 시간을 추정했다. 수집된 모든 대변 샘플(n = 484)에서 대변 수분 함량(대변 수분의 지표, 통과 시간의 대리 지표23), pH 및 미생물 부하를 측정했습니다 .

수집된 모든 소변 샘플(일일 스팟 및 식후 샘플, n = 1,154)과 일부 대변 샘플(n = 170)은 비표적 액체 크로마토그래피-질량 분석법(LC-MS) 대사체학 분석을 통해 소변 및 대변 대사체 프로파일을 얻었다.

마지막으로, 일부 대변 샘플(n = 362)의 16S 리보솜 RNA(rRNA) 유전자 시퀀싱을 통해 장내 미생물군 구성 정보를 얻었으며, 미생물 부하를 조정한 후 상대적 미생물군 프로파일(RMP)과 정량적 미생물군 프로파일(QMP)을 모두 평가하였다24.

Table 1 Participants’ characteristics (N = 61)

Full size table

Fig. 1: PRIMA study design and variations in gut environmental factors, gut microbiome and metabolomes.

a, PRIMA study design. The study included two site visits, at which fasting blood and breath samples were taken. At visit 1, anthropometric measurements were attained, and all participants were given a standardized breakfast; a subset of 50 volunteers ingested SmartPills immediately after. Postprandial breath hydrogen and methane were measured every 30 min for 6 h, and postprandial urine was collected at 0.5 h and every hour until 24 h as indicated. On days 3 and 5, participants performed a sweet-corn test to measure WGTT. In addition, daily 24 h dietary records (days 1–8), records of bowel habits (stool consistency, stool frequency and time of defecation) and daily urine and faecal samples were obtained. Solid line indicates sample collection on site and dashed line sample collection at home. b, Inter- and intra-individual variations in the gut environmental factors over the 9 consecutive days. The red and blue lines represent median and mean values, respectively. Grey lines represent intra-individual fluctuations over time. Asterisks indicate the statistical significance of mixed-effect models accounting for repeated measures (two-sided ***P < 0.001, **P < 0.01, *P < 0.05; NS, not significant; see Supplementary Table 3 for details; no adjustment for multiple testing was applied). c, Percentage of variation explained by individual and study day in the gut microbiome and urine and faecal metabolomes based on PERMANOVA tests (two-sided *P < 0.05).

Full size image

Gut environment stability varies for each individual

Daily sampling allowed us to eval‎uate the fluctuations in gut environmental factors, faecal and urine metabolomes, gut microbiomes and diets within and between healthy adults over time (Supplementary Fig. 1). First, we observed varying degrees of day-to-day fluctuations within individuals for faecal pH (coefficient of intra-individual variation (CVIntra) 0.3–8.1%), BSS (0–57.8%), stool frequency (0–73.1%), stool moisture (2.2–24%) and microbial load (7.6–72.7%) (Fig. 1b and Supplementary Table 2), suggesting that some individuals have more stable gut environments than others. Most of the gut environmental factors varied within individuals over the 9 days, whereas faecal pH remained relatively stable (Fig. 1b and Supplementary Table 3a). Participant ID explained a significant proportion of variation for day-to-day fluctuations in all of the gut environmental factors (Supplementary Table 3b), indicating that the stability of the gut environment is to some extent personal.

Next, we performed a permutational multivariate analysis of variance (PERMANOVA) on the QMPs, urine and faecal metabolomes and found that the individual explained more than 50% of the inter-individual variations in all three cases (Fig. 1c). The sampling day explained on average 6.7% of the urine metabolome variation but did not explain day-to-day variations in the gut microbiome and faecal metabolome. (Fig. 1c). By inspecting the β-diversities of individual microbiome and metabolome profiles, we observed that some individuals showed less variation over the study period than others (Extended Data Fig. 2).

Stool moisture and pH explain daily gut microbiome fluctuations

To explore what drives the intra-individual fluctuation in the metabolomes and the microbiome, we performed distance-based redundancy analysis (db-RDA). We considered daily dietary macronutrients and fibres, as well as the gut environmental factors. None of the dietary components explained intra-individual fluctuations in the gut microbiome or metabolomes. By contrast, stool moisture, faecal pH, BSS and time of defecation markedly affected the gut microbiome (QMP, genus level; Fig. 2a), explaining 3.5%, 2.5%, 2% and 1.3% of the variations, respectively. Similar results were observed using the RMP data (Extended Data Fig. 3a) and in previous studies6,25. Notably, these explanatory factors are proxies for gut transit time, suggesting that day-to-day variations in transit time are reflected in the gut microbiome variation.

개체마다 장 환경 안정성이 다름

일일 샘플링을 통해 건강한 성인 내·간에서 시간에 따른

장 환경 요인, 대변·소변 대사체, 장내 미생물군집 및 식단의 변동성을 평가할 수 있었다(보충 그림 1).

첫째, 개체 내 대변 pH(개체 내 변동계수 (CVIntra) 0.3–8.1%), BSS (0–57.8%), 배변 빈도 (0–73.1%), 대변 수분 함량 (2.2–24%) 및 미생물 부하 (7.6–72.7%)에서 관찰되었으며(그림 1b 및 보충 자료 표 2), 이는 일부 개인이 다른 개인보다 더 안정적인 장 환경을 가질 수 있음을 시사한다. 대부분의 장 환경 요인은 9일 동안 개인 내에서 변동한 반면, 분변 pH는 상대적으로 안정적으로 유지되었다(그림 1b 및 보충 자료 표 3a). 참가자 ID는 모든 장 환경 요인의 일일 변동에서 상당한 변동 비율을 설명하였다(보충 자료 표 3b). 이는 장 환경의 안정성이 어느 정도 개인에 따라 다르다는 것을 시사한다.

다음으로, QMP, 소변 및 대변 대사체에 대해 순환 다변량 분산분석(PERMANOVA)을 수행한 결과, 세 가지 경우 모두에서 개인이 개인 간 변이의 50% 이상을 설명하는 것으로 나타났다(그림 1c). 채취일은 소변 대사체 변이의 평균 6.7%를 설명했으나, 장내 미생물군집과 대변 대사체의 일일 변동은 설명하지 못했다(그림 1c). 개별 미생물군집 및 대사체 프로파일의 β-다양성을 분석한 결과, 일부 개인은 연구 기간 동안 다른 개인보다 변동성이 적게 나타났다 (확장 데이터 그림 2).

Fig. 2: Intra- and inter-individual variations in gut microbiome and urine metabolome explained by gut environment.

a,b, Contributions of dietary and gut environmental factors on intra-individual variations in gut microbiome (QMP, all days) (a) and urine metabolome (all days, all features) (b). c, Boxplots showing segmental and WGTT measured by the SmartPill (n = 50) at day 2 and mean transit time of sweet corn (n = 61, day 3 and day 5) with each dot representing an individual. d, Boxplots showing pH throughout the gastrointestinal tract measured by the SmartPill (n = 50) and in faeces measured by pH meter at day 2 (n = 61) with each dot representing an individual. e,f, Contributions of clinical variables, dietary components, gut environmental and physiological factors and participant characteristics to inter-individual variations in the gut microbiome (QMP, sample closest to the capsule body exit) (e) and urine metabolome (24 h, day 2, all features) (f). Panels a, b, e and f were quantified by db-RDA with permutation tests using Bray–Curtis distances. Effect sizes are plotted. The asterisks indicate statistical significance after adjustments for multiple testing (*q < 0.1). See Extended Data Fig. 3a,b for RMPs and faecal metabolome. Boxplot centre in c and d represents median, and box represents interquartile range (IQR). Whiskers extend to most extreme data point <1.5 IQR. SB, small bowel; CH4, breath methane; H2, breath hydrogen.

Source data

Full size image

Stool moisture and faecal pH further explained 3.1% and 3%, respectively, of the intra-individual variation in urine metabolomes, despite subtle day-to-day fluctuations (Fig. 2b). This suggests that even small changes in the colonic water content and pH may be associated with the host–microbiota metabolism. However, these observations could also be influenced by daily diet variations. Gut environmental factors did not contribute to the intra-individual fluctuations in the faecal metabolomes (Extended Data Fig. 3b). It should be noted that faecal metabolome data were derived only from three consecutive days, and stool moisture still tended (P = 0.081) to have an effect.

Transit times and pH vary between individuals

SmartPills were used to obtain whole-gut transit time (WGTT), gastric emptying time (GET), small-bowel transit time (SBT), colonic transit time (CTT) and intestinal transit time (ITT; SBT + CTT), as well as pH throughout the gastrointestinal tract (GIT) (Extended Data Fig. 3c). In 8 individuals, the capsule was retained in the stomach for over 8 h, a common event reported in other studies19,26. Therefore, GET and WGTT values from these participants were excluded from our analyses. In addition, we could not determine CTT and WGTT in one participant due to a signal loss.

The median values of transit time were as follows: GET, 4.8 h (range 3.1–6.2 h); WGTT, 23.3 h (12.4–72.3 h); CTT, 13.6 h (2.1–63.5 h); and SBT, 5.1 h (2.5–10.3 h), in agreement with previously reported data on healthy populations27. For comparison, the corn TT showed a median of 23.6 h (10.8–109.7 h) at day 3 and 19.7 h (12.0–84.5 h) at day 5. Furthermore, we found a strong correlation between the two corn TT measurements (Spearman correlation coefficient (SCC) = 0.8, P < 0.001) suggesting consistency within individuals. The median of the mean corn TT across the two days was 21.7 h (11.7–97.1 h) (Fig. 2c), similar to the WGTT obtained by the SmartPill. However, we did not observe any correlation between the WGTT and corn TT (Extended Data Fig. 4a), indicating that despite providing similar results on average, individually, these two methods showed different results.

When exploring the relationships between segmental transit times, corn TT, gut environmental factors and participant characteristics (Extended Data Fig. 4a), we found that the transit times recorded by both methods were negatively correlated to BSS, as also reported previously23,28. We also observed that women had significantly longer CTT compared with men, while there was no effect of menstruation status among the women (Extended Data Fig. 4b).

Large inter-individual variations in the gastrointestinal segmental pH were also observed (Fig. 2d) with the following median pH values in the upper GIT: the stomach (0.9, range 0.5–4.9), duodenum (6.1, 5.0–7.2) and small intestine (7.4, 6.4–8.2). pH in the proximal colon was slightly acidic (6.3, 5.3–7.0) followed by a gradual increase in the distal colon (6.9, 5.0–8.2) and sigmoid colon (7.2, 5.6–8.6). Interestingly, a small decrease in pH was observed from the sigmoid colon to the rectum (7.0, 5.7–8.6) and also in the faecal pH (6.9, 6.6–7.3), indicating that acidifying processes occur after entry into the rectum.

CTT and pH contribute to gut microbiome variations

To quantify how participant characteristics, clinical variables, diet and gut environmental and physiological factors explain inter-individual variations in the gut microbiome and metabolomes, we performed a db-RDA using data derived from faecal and 24 h urine collections on day 2 from all participants (n = 61; Supplementary Table 4). Moreover, we performed the same analysis with whole-gut and segmental transit times and pH derived from the SmartPills on day 2 (n = 50). Stool moisture and distal colon pH were key factors associated with inter-individual variation in QMPs (Fig. 2e), accounting for 5.5% and 5% of the variation, respectively, on day 2 and also other days (Supplementary Table 4). Unlike previously reported data from larger cohorts5, BSS did not explain a significant proportion of the variation in QMP in this population. WGTT, CTT, corn TT and faecal pH explained 9.1%, 6.2%, 4.9% and 5.4%, respectively, of the inter-individual variations in the 24 h urine metabolome, in comparison to age, which explained 4.5% of the variation (Fig. 2f). These contributions were consistent when testing against the urine metabolomes on different days (Supplementary Table 4). By contrast, segmental transit time did not contribute to the inter-individual variation in the faecal metabolomes, whereas pH in the distal colon and fibre intake showed the largest effects explaining 6.8% and 5.9% of the variations, respectively; however, this was not significant after adjusting for multiple testing (Supplementary Table 4).

We also tested the effect of menstruation during the study period for women (non-menstruating, n = 30; menstruating, n = 13), which showed effect sizes of 3.7% (gut microbiome) and 3.6% (urine metabolome), however without statistical significance. Considering the notable effect size of age on urine metabolome and a significant age difference between the two groups of women (P = 0.01), age might contribute to these observed effects. Our results emphasize that the personal gut environment contributes considerably to the inter-individual differences in the gut microbiota and urinary metabolic profiles.

Individual gut microbiota and metabolite profiles are dynamic

We next assessed intra-individual fluctuations in microbial-derived metabolites including breath hydrogen and methane, faecal SCFAs (acetate, propionate, butyrate, valerate and caproate) and BCFAs (isobutyrate, isovalerate and 2-methylbutyrate; Table 1), as well as 16 other microbial-derived metabolites detected in faeces and urine, including the proteolytic markers, p-cresol sulfate (PCS), phenylacetylglutamine (PAGln) and indoxyl sulfate. Substantial day-to-day fluctuations were observed (Fig. 3a and Supplementary Fig. 2). Breath methane and hydrogen had a median CVIntra of 141% and 47%, respectively, with a moderate positive correlation between the two time points for both gases (hydrogen, SCC = 0.42, P < 0.001; methane, SCC = 0.66, P < 0.001). Faecal concentrations of the SCFAs and BCFAs fluctuated considerably from day to day (median CVIntra ranging from 26% to 40%) with valerate varying the least and acetate the most. Similarly, the relative abundances of the proteolytic markers varied substantially from day to day with a median CVIntra of 26%, 42% and 39% for PAGln, indoxyl sulfate and PCS, respectively. These findings suggest that microbial-derived metabolites in breath, faeces and urine fluctuate from day to day on a habitual diet.

Fig. 3: Fluctuations in microbial metabolites and their correlations to gut physiology and environment.

a, Intra-individual fluctuations in microbial metabolites measured in breath, faeces and urine. Boxplots show coefficients of intra-individual variations. Each dot represents an individual (n = 61). Boxplot centre represents median, and box represents IQR. Whiskers extend to most extreme data point <1.5 IQR. b,c, Correlations between microbial metabolites and gut physiology and environmental factors as assessed by repeated measures correlation coefficient (b) or the Spearman correlation coefficient (c). The asterisks indicate statistical significance after adjustment for multiple testing (****q < 0.001, ***q < 0.01, **q < 0.05, *q < 0.1). Blue, brown and yellow bars indicate breath, faecal and urine metabolites, respectively. The black bar in b indicates repeated measure correlations where daily values for each variable have been used (Extended Data Fig. 5), whereas the green bar in c indicates analysis based on data collected on day 2. Postprandial hydrogen and methane were only measured at one time point and therefore were not included in the repeated measure analysis. (f), faecal; TMAO, trimethylamine N-oxide.

Source data

Full size image

Next, we used repeated measures (Fig. 3b, Extended Data Fig. 5 and Supplementary Table 5) and Spearman correlation analysis (Fig. 3c) to find links between specific microbial metabolites and the gut environment along with diet. Faecal SCFAs were negatively correlated to faecal pH with butyrate showing the strongest correlation (r = −0.77, q < 0.001) in line with previous human studies29. Moreover, higher faecal propionate (SCC = −0.25, q < 0.1) was linked to shorter CTT with a similar tendency observed for faecal butyrate (SCC = −0.29, P < 0.05, q = 0.3). Faecal butyrate also tended to negatively correlate with rectal pH (SCC = −0.37, P < 0.05, q = 0.2) but not with pH in other segments of the colon, suggesting that butyrate production may contribute to the reduced pH observed in the rectum and faeces.

Proteolytic markers including urinary PAGln and faecal BCFAs were negatively correlated to stool moisture, and urinary PCS correlated positively with faecal pH (PAGln, r = −0.12, q < 0.1; isobutyrate, r = −0.39, q < 0.05; isovalerate, r = −0.37, q < 0.1; 2-methylbutyrate, r = −0.43, q < 0.05; PCS, r = 0.12, q < 0.1). Similarly, higher urinary levels of PCS were associated with longer CTT and ITT (SCC = 0.48, SCC = 0.44, respectively, q < 0.05), with similar tendencies observed for PAGln and indole-lactic acid. Furthermore, breath methane was linked to lower stool moisture and longer CTT. It is worth noting that none of this was shown for SBT, indicating that CTT determines the abundance of these metabolites and supports the hypothesis that longer passage through the colon is linked to microbial proteolysis possibly due to the depletion of substrates for saccharolytic fermentation9,30. Urinary indoxyl-glucuronide was positively associated with pH in the distal colon (SCC = 0.33, q < 0.1), and a similar trend was found between urinary PCS and rectal pH (SCC = 0.32, P < 0.05). These metabolites did not correlate to pH in the small intestine and the proximal colon, indicating a higher contribution of microbial proteolysis to pH in the distal gut compared with the proximal gut. In summary, CTT and colonic pH, but not SBT and small-intestinal pH, are associated with levels of several microbial metabolites in breath, faeces and urine. In addition, we found several associations between microbial metabolites and dietary components (Supplementary Table 5) with notable inverse correlations between the intake of dietary fibres and faecal BCFAs, urinary PCS and urinary PAGln, respectively.

Faecal and urine metabolomes are linked to transit time and pH

To explore unknown metabolic features related to gut physiology, we used untargeted metabolomics to profile the urine and faecal metabolomes. We applied univariate and multivariate statistical models on all molecular features identified in urine and faeces. We first used sparse partial least squares (SPLS) models on the SmartPill-derived data and urine metabolomes from 24 h postprandial urine collected on day 2 and faecal metabolomes collected closest to the SmartPill egestion. We then performed linear regression models on the same data and further investigated features selected by both models (446 unique features; Supplementary Table 6).

Several metabolic features in urine and faeces were associated with whole-gut and segmental transit time and pH (Fig. 4a,b). To investigate these features in further detail, the corresponding samples were analysed by tandem mass spectrometry (MS2) and by matching with authentic standards when available, resulting in the identification of 33 metabolites (Supplementary Tables 7 and 8).

Fig. 4: Metabolites identified via LC-MS untargeted metabolomics associated with segmental transit time and pH.

a,b, Volcano plots derived from regression models where each dot represents a metabolic feature with blue representing statistically significant associations (FDR-adjusted P < 0.1) in urine (a) and faeces (b). The x axis shows the regression coefficient values (estimate) indicating either positive or negative associations and the y axis represents the −log10 transformation of the P value.

Source data

Full size image

Apart from urinary levels of PCS and PAGln, several other urinary metabolites derived from the breakdown of aromatic amino acids tryptophan and tyrosine by gut microbes were found to be linked with gut transit time and faecal pH. Specifically, 5-hydroxy-2-oxindole sulfate, 3-hydroxy-2-oxindole sulfate and 4-hydroxybenzoic acid sulfate were associated with longer WGTT/CTT, while 3-hydroxy-2-oxindole glucuronide correlated with higher faecal pH. By contrast, faecal tryptophan was negatively linked to faecal pH. In addition, higher faecal proline and urinary picolinoylglycine levels were linked with increased faecal and rectal pH, respectively.

Several dicarboxylic acids in faeces, pimelic, suberic and sebacic acids were positively associated with WGTT and CTT. By contrast, faecal glutaric acid and pipecolic acid were negatively correlated with WGTT/CTT and sigmoid, rectal and faecal pH. Pipecolic acid is highly abundant in plants; however, it can also be produced by the gut microbiota from lysine31. Furthermore, higher urinary levels of citric acid were positively associated with pH in the proximal colon.

Moreover, faecal levels of 2-oxindole-3-acetic acid, previously linked to the New Nordic Diet and Mediterranean diet32,33, were negatively associated with WGTT, CTT and faecal pH. Similarly, faecal pantothenic and nicotinic acids were negatively associated with CTT and faecal pH, respectively. In addition, dihydroferulic acid glucuronide and argininic acid in urine were negatively associated with rectal pH, while p-hydroxyphenyllactic acid in faeces was negatively linked to faecal pH.

4-Hydroxyhippuric acid and several urinary markers of coffee intake, including 1-methyluric acid, 1-methylxanthine, 1,3-dimethyluric acid, 1,7-dimethyluric acid and 1,3,9-trimethyluric acid, were negatively associated with small-intestinal pH. 1-Methylxanthine and 1,3,9-trimethyluric acid in faeces were also negatively associated with WGTT or faecal pH, suggesting a link between coffee consumption and gut function. In addition, a positive correlation was observed between rectal pH and urinary 4-methylcatechol sulfate, a metabolite of quercetin found in plant-based foods34. Urinary taurine and faecal cholic acid were also positively associated with small-intestinal pH, supporting the role of bile acids in neutralizing the acidic chyme coming from the stomach35.

Finally, urinary pseudouridine, a primary constituent of RNA, was found to be inversely associated with CTT and sigmoid colon pH, in line with our previous work9. Pseudouridine was also found in faeces and showed a similar inverse relationship with faecal pH, as did deoxy-xanthosine and xanthine. This suggests a link between increased cell turnover and lower colonic pH.

Altogether, by using untargeted LC-MS metabolomics, we identified several host-, microbial- and food-derived metabolites associated with WGTT, CTT and pH in the distal part of the colon emphasizing an interplay between diet, the gut environment, the host and the microbiota.

Microbial alpha diversity is linked to long passage

To explore potential links between the identified metabolites and the gut microbiota, Spearman correlation analysis was performed (Fig. 5). Strong positive correlations between microbial alpha diversity measures and microbial proteolysis, CTT and ITT were found. On the contrary, alpha diversity correlated negatively with stool moisture and microbial saccharolysis.

Fig. 5: Associations between metabolites, bacterial genera and gut environmental factors.

Hydrogen and methane were measured in breath, and SCFAs and BCFAs were measured in faeces; other metabolites as indicated in brackets. (b), breath; (f), faecal metabolite; (u), urine metabolite. The colour gradient shows the Spearman correlation coefficient, and the asterisks indicate statistical significance after adjustment for multiple testing (***q < 0.001, **q < 0.01, *q < 0.05).

Source data

Full size image

Products of microbial proteolysis and dicarboxylic acids were positively correlated with the absolute abundances of several bacterial genera including Intestimonas, Flavonifractor, Eubacterium, Lachnospira, Clostridium, Oscillibacter, Alistipes, Dialister and Akkermansia. The same genera negatively correlated with faecal levels of tryptophan, oxindole-3-acetic acid and various coffee-derived metabolites. Not surprisingly, these genera were also positively associated with longer ITT and CTT and higher faecal pH, and negatively associated with stool moisture and/or BSS. Conversely, SCFAs-producing genera including Agathobacter, Faecalibacterium and Blautia36,37, along with lactate-producing Streptococcus, were all positively associated with faecal nicotinic acid, pantothenic acid and the coffee-derived metabolites. Notably, Oscillibacter, Alistipes and Akkermansia have repeatedly been found elevated in samples linked to longer transit time and/or constipation8,9,28, whereas butyrate-producing genera including Faecalibacterium and Agathobacter have been associated with shorter transit time28,38. In summary, these observations highlight the interdependency between gut bacteria, metabolites and gut physiology.

Discussion

Gut transit time and pH are important determinants of gut microbiota composition and metabolism7. Here we showed substantial variation in whole-gut and segmental transit time, along with luminal pH among healthy individuals. These variations explained differences in microbiome composition and host–microbiota co-metabolism. As pH and transit time influence microbial growth and enzyme activities39, these factors could play a key role in shaping the gut microbial composition and metabolism along the GIT as well as microbiome responses to foods. A recent study confirmed that microbiome and metabolome compositions differ along the GIT40. Future studies with sampling along the GIT combined with measurements of regional pH and transit time are needed to ultimately disentangle this. Our study emphasizes that person-specific differences in the luminal pH may pose challenges for studies using pH-sensitive ingestible devices40,41. It is worth noting that repeated measurements in human studies are needed as we observed substantial daily fluctuations in microbial-derived metabolites.

We identified several metabolites associated with longer transit time that have been reported to be elevated in patient groups with constipation42,43,44. Among these, dicarboxylic acids including pimelic acid were associated with longer ITT and/or higher pH, and bacteria consistently associated with constipation8,9,28. Pimelic acid, possibly originating from microbial fatty acid metabolism45,46, has been found at higher faecal levels in patients with chronic kidney disease47 and colorectal cancer48, often associated with constipation49,50. A recent study showed an increased abundance of dicarboxylic acids towards the distal gut40, and the authors speculated that it could be due to the catabolism of host epithelial cells. Whether longer ITT might be associated with increased epithelial cell turnover and shedding needs further research.

토론

장 통과 시간과 pH는

장내 미생물 군집 구성과 대사의 중요한 결정 요인이다7.

본 연구에서는

건강한 개인들 사이에서 전체 장 및 구간별 통과 시간과 관강 내 pH에 상당한 변이가 있음을 보여주었다.

이러한 변이는

미생물군 구성과 숙주-미생물 공동 대사의 차이를 설명하였다.

pH와 통과 시간은

미생물 성장과 효소 활성에 영향을 미치므로39,

이러한 요인들은 위장관(GIT)을 따라 장내 미생물 구성과 대사,

그리고 음식에 대한 미생물군 반응을 형성하는 데 핵심적인 역할을 할 수 있습니다.

최근 연구는

미생물군과 대사체 구성이 위장관(GIT)을 따라 다르다는 점을 확인했습니다40.

이를 궁극적으로 규명하기 위해서는

위장관(GIT)을 따라 샘플링과 함께 지역별 pH 및 통과 시간 측정을 결합한 향후 연구가 필요합니다.

본 연구는

인체별 내강 pH 차이가 pH 감지형 섭취형 장치40,41을 활용한 연구에

어려움을 초래할 수 있음을 강조한다.

미생물 유래 대사산물의 상당한 일일 변동성을 관찰했기에

인간 대상 연구에서 반복 측정이 필요하다는 점은 주목할 만하다.

우리는

변비 환자 집단에서 증가된 것으로 보고된42,43,44

장 통과 시간과 연관된 여러 대사 산물을 확인했습니다.

이 중 피멜산(pimelic acid)을 포함한 디카르복실산(dicarboxylic acids)은

더 긴 ITT 및/또는 더 높은 pH와 연관되었으며,

변비와 일관되게 연관된 박테리아8,9,28과도 연관되었습니다.

미생물 지방산 대사에서 유래했을 가능성이 있는 피멜산45,46은

만성 신장 질환47 및 대장암48 환자에서

배변 시 더 높은 농도로 발견되었으며,

이는 종종 변비와 연관되어 있습니다49,50.

최근 연구에서는 원위부 장으로 갈수록

디카르복실산의 풍부도가 증가하는 것으로 나타났으며40,

저자들은 이것이 숙주 상피 세포의 이화 작용 때문일 수 있다고 추측했습니다.

ITT가 길어질수록

상피 세포의 회전율과 탈락이 증가하는지 여부는 추가 연구가 필요합니다.

A negative association between daily fibre intake and several proteolytic markers was also observed. As dietary fibres can regulate microbial tryptophan metabolism51, availability of fibre in the colon may affect microbial protein fermentation associated with negative health outcomes11,12. Further research is needed to understand these mechanisms and to explore the relationship between gut physiology and microbiome under controlled diets, possibly involving dietitians or providing whole diets.

Despite its limited cohort size, our study shows significant associations between intestinal segmental transit time and pH with intra- and inter-individual differences in the gut microbiome composition and metabolism in a healthy population. Potential limitations to consider are the choice of breakfast made before the SmartPill measurement; given that past validation studies and their normative data rely on specific meal/nutrient combinations52,53, any deviations from these could likely influence gut motility and transit time. Furthermore, the sweet-corn test is not a validated tool to assess WGTT despite being cost effective. Finally, the introduction of corn and the meals during the first visit constitute small dietary changes, but they could possibly have impacted gut physiology. While this study included a rather homogenous group of healthy volunteers, it provides valuable insights into longitudinal variations in gut microbial metabolism and pH over more than 1 week. Our results highlight the important role of transit time and pH for the gut microbiome composition and levels of microbial-derived metabolites, emphasizing the importance of considering gut physiology and environment in human microbiome studies. This may be key for understanding the healthy gut microbiome and for disentangling personal microbiome responses to foods and other lifestyle factors.

일일 섬유질 섭취량과 여러 단백질 분해 마커 사이의 음의 상관관계도 관찰되었다. 식이 섬유는 미생물 트립토판 대사를 조절할 수 있으므로51, 결장에서 섬유질의 가용성은 건강에 부정적인 결과를 초래하는 미생물 단백질 발효에 영향을 미칠 수 있다11,12. 이러한 메커니즘을 이해하고 통제된 식이 조건 하에서 장 생리학과 미생물군집 간의 관계를 탐구하기 위해서는 영양사 참여 또는 전체 식이 제공을 포함한 추가 연구가 필요하다.

코호트 규모가 제한적임에도 불구하고, 본 연구는 건강한 인구 집단에서 장 구간 통과 시간 및 pH와 장내 미생물군집 구성 및 대사에서의 개인 내·간 차이에 대한 유의한 연관성을 보여줍니다. 고려해야 할 잠재적 한계로는 스마트필 측정 전 선택한 아침 식사가 있습니다. 과거 검증 연구 및 그 기준 데이터가 특정 식사/영양소 조합에 의존한다는 점을 감안할 때52,53, 이를 벗어난 경우 장 운동성과 통과 시간에 영향을 미칠 수 있습니다. 또한, 스위트콘 테스트는 비용 효율적임에도 불구하고 WGTT 평가를 위한 검증된 도구가 아닙니다. 마지막으로, 첫 방문 시 옥수수 및 식사 도입은 작은 식이 변화이지만 장 생리학에 영향을 미쳤을 가능성이 있습니다. 본 연구는 비교적 동질적인 건강한 자원자 집단을 포함했지만, 1주 이상에 걸친 장내 미생물 대사 및 pH의 종적 변동에 대한 귀중한 통찰력을 제공합니다. 본 연구 결과는 장내 미생물군 구성과 미생물 유래 대사산물 수준에 있어 통과 시간과 pH의 중요성을 부각시키며, 인간 미생물군 연구에서 장 생리학적 환경을 고려할 필요성을 강조한다. 이는 건강한 장내 미생물군을 이해하고, 식품 및 기타 생활습관 요인에 대한 개인별 미생물군 반응을 규명하는 데 핵심적일 수 있다.

Methods

Study participants

A 9 day human observational trial (PRIMA, toward Personalized dietary Recommendations based on the Interaction between diet, Microbiome and Abiotic conditions in the gut) among healthy participants was conducted at the Department of Nutrition, Exercise and Sports at the University of Copenhagen in Denmark from April to December 2021. The research protocol was approved by the Municipal Ethical Committee of the Capital Region of Denmark (H-20074067), and all participants provided written informed consent to participate according to Case Report (CARE) guidelines and in compliance with the principles of the Declaration of Helsinki. The study was registered at ClinicalTrails.gov (ID NCT04804319).

Out of the initially anticipated 85 individuals, 63 healthy participants living in Denmark were enrolled, and 61 completed the study (43 women and 18 men; Extended Data Fig. 1). The two drop-outs were excluded due to illness and antibiotic administration. Among the 61 study participants, a subset of 50 volunteers (37 women and 13 men) underwent wireless motility capsule monitoring at visit 1 as anticipated. Volunteers were compensated with gift cards (500 or 800 DKK) but received no direct financial compensation. The criteria for inclusion in the study specified for participants who were healthy by self-report (did not suffer from inflammatory bowel syndrome, small-intestinal overgrowth, inflammatory bowel disease, chronic or infectious disease, diabetes or cancer), ages 18–75 years old with a BMI between 18.5 and 29.9 kg m−2, with no intake of medication, except for mild antidepressants and contraceptive pills. Intake of antibiotics, diarrhoea inhibitors and laxatives 1 month before the trial was not allowed. Furthermore, pregnant or lactating women were not included in the trial.

The PRIMA study was an explorative study. The primary outcome was to investigate associations between faecal pH and gut microbial saccharolytic/proteolytic metabolism (assessed by targeted and untargeted metabolomics). The secondary outcomes included to explore relationships between the gut environmental factors (small-intestinal and colonic pH and transit time assessed by wireless motility capsules and various transit time proxy markers measured in faeces) and gut microbiome and metabolome assessed by 16S rRNA sequencing and metabolomics.

Experimental design and sample collection

Seven days before the study, the participants were asked not to consume any sweet corn as two self-administered sweet-corn tests to eval‎uate the WGTT were part of the study. Before both visits, the participants were asked to abstain from alcohol intake, smoking and strenuous exercise.

The participants were asked to maintain their habitual diet and register their food intake online via the Myfood24 tool (myfood24.org) with nutritional values based on the Danish food composition database FRIDA version 4.1 (frida.fooddata.dk) for eight consecutive days during the study. During the trial, the participants collected daily stool samples (first bowel movement whenever possible), stored the samples in their domestic freezers and transported them to the laboratory while being kept cold. Moreover, the participants self-reported daily their defecation patterns including time of defecations, stool consistency assessed by the BSS and stool frequency, intake of dietary supplements and medication (limited to pain killers for a few participants), and their gastrointestinal symptoms. The gastrointestinal symptoms were assessed based on a 10 cm visual analogue scale (0, no symptoms; 10, the most severe symptoms) in regard to stomach ache, bloating, constipation, diarrhoea and overall comfort. Women were asked to note down whether they had menstruation during the study period (yes/no). Furthermore, the participants collected seven daily spot morning urine samples (days 1, 2, 4, 5, 6, 7, 8; the first morning sample) and two 24 h urine samples (days 2–3 and days 8–9) during the study period. The collected urine samples were stored in the participants’ domestic freezers, transported to the study site in a cooling bag and stored at −20 °C overnight. After thawing at 5 °C, aliquots of 1 ml were taken and stored at −80 °C until further use. In addition, the participants consumed 100 g of sweet corn before their evening meal on days 3 and 5 and recorded the time of the corn egestion23.

At both visits (days 2 and 9), fasting blood and breath samples were collected. During the first visit, anthropometric measurements (height, body weight and BMI) were obtained. Furthermore, the first visit also included a standardized meal test for all participants (n = 61). The test meal consisted of rye bread (with butter and jam), a boiled egg, a portion of natural yoghurt along with nuts, walnuts, blueberries and a glass of water (100 ml) with 250 mg of dissolved paracetamol (Table S1), which was used as a marker of postprandial gastric emptying of liquids54. The meal portion size was calculated as 25% of the daily energy demand of each participant based on the Harris–Benedict equation22. Postprandial urine samples (at 30 min, 60 min, 120 min, 180 min, 240 min, 300 min and 360 min and between 6 and 8 h, between 8 and 10 h and between 10 and 24 h) and postprandial breath exhalations (at 30 min, 60 min, 90 min, 120 min, 150 min, 180 min, 210 min, 240 min, 270 min, 300 min, 330 min and 360 min) were collected. A subset of participants (n = 50) ingested a SmartPill capsule immediately after the meal with a bit of additional water if needed. All participants drank 150 ml of water at 2 h and 4 h after the meal, respectively. At 6 h, all participants received a sandwich and 500 ml of water and left the study site.

SmartPill data collection and analysis

The SmartPill capsule is a single-use wireless gastrointestinal capsule (26.8 mm × 13 mm) that transmits data on luminal pH, temperature and pressure to a portable receiver, which was worn by the participants from ingestion to egestion and thereafter returned to the study personnel. The capsule measures a pH range of 1–9, with an accuracy of ±0.5 pH units, pressure at a range of 0–350 mmHg (±5 mmHg) and temperature ranging between 20 °C and 40 °C (±1 °C)23. Upon receiving the portable receiver, the study personnel downloaded the raw data from the receiver to the manufacturer’s software via a docking station (Motility GI v 3.1). Intestinal segmental transit times were determined based on landmark changes in the pH values as follows: gastric emptying was defined as the time point with an abrupt increase of ≥3 pH units indicating passage from the stomach into the duodenum. The passage from the small intestine into the ileocaecal junction was defined as the first time point with a decrease of at least one pH unit. The body exit of the capsule was defined as the time point with a decrease in temperature and/or a loss of data. The time of capsule residence in each of the gastrointestinal segments corresponds to GET, small-intestinal transit time, CTT and combined WGTT. Regional pH and pressure profiles were also obtained, and the median values were determined. The segmental transit time and pH values in the colon were further segmented into proximal, distal and recto-sigmoid, respectively. The proximal colon pH and transit time were estimated as median values of the first 32.3% of the total CTT, the distal colon pH were median values of the next 32.6% and the recto-sigmoid pH were median values of the last 35.4%; this was based on previously reported data, which determined the percentages of total CTT according to the location of radio-opaque markers (visualized by X-rays) in the different segments of the colon16. In addition, the median pH value measured during the last 10 min before the capsule egestion was registered as rectal pH.

Dietary records

Detailed 24 h weighted food intakes were recorded for 8 consecutive days by the participants via the online Myfood24 tool (myfood24.org) with nutritional values based on the Danish food composition database FRIDA version 4.1 (frida.fooddata.dk). The collected data included information about the intake of macronutrients (carbohydrate, protein, fat) and dietary fibre (AOACFIB), in addition to information about more than 80 nutrients. Under-reporting was identified by calculating the reported caloric intake divided by the average daily energy demand for each person with a cut-off value of 0.8 (ref. 55). Accordingly, approximately 25% of the daily dietary records were under-reported, and the data were removed in the subsequent analyses in this study (this essentially affected 10 participants who under-reported more than 4 out of the 8 days, while the other participants occasionally under-reported daily intakes). By contrast, no over-reporters (cut-off >2.5) were detected. The total dietary profiles (all macro- and micronutrients available in Myfood24) were used in the principal component analysis, whereas macronutrient and fibre intake were used in the redundancy analyses. The daily intake was used for intra-individual analysis, whereas mean intake across the 8 days was used for the inter-individual analysis.

Breath exhalations measurements

Fasting and postprandial levels of hydrogen and methane were measured in all breath samples by the M.E.C. Lactotest 202 Xtend device (M.E.C. R&D sprl).

Biochemical analysis of blood

Blood samples were immediately put on ice upon collection until they were centrifuged for precipitation of blood cells and stored at −80 °C. Glucose was measured in plasma samples by using Pentra ABX 400 (HORIBA ABX) with a detection limit of 0.11 mmol l−1. Serum insulin and C-peptide levels were measured by using Immulite 2000 XPi (Siemens Healthcare Diagnostics) with the detection limit of 14.4 pmol l−1 and 27 pmol l−1, respectively. Before the analyses, both instruments’ performances were validated using external and internal insulin, c-peptide and glucose controls. Three participants arrived for the second visit in a postprandial state; the blood was collected and analysed accordingly, but the glucose, insulin and c-peptide values were not included in the data analysis.

Faecal measurements

Upon receipt, faecal samples were stored at −20 °C overnight, thawed and homogenized in sterile water with a sample-to-water ratio of 1:1 (w/v) (faecal slurry). Subsequently, pH was measured in the faecal slurry using a digital pH meter (Mettler Toledo). The homogenized samples were subsequently aliquoted to cryotubes and stored at −80 °C until further analyses. Stool moisture was determined by evaporating the water of one aliquot (approximately 1 ml) using a vacuum concentrator (Speed-Vac, Christ RVC 2-25) and by calculating the faecal weight difference before and after evaporation.

Faecal SCFAs and BCFAs were quantified by LC-MS in samples collected between day 2 and day 5 (n = 170) as previously described32. In brief, the samples were thawed, mixed with ethanol and purified by filtration (0.2 µm filter). Subsequently, the samples were derivatized with 3-nitrophenylhydrazine, and labelled internal SCFA standards were added. Dilution series of external SCFA standards were spiked with internal SCFA standards, and all derivatized samples were analysed on ultra-performance liquid chromatography (UPLC)-quadrupole time-of-flight mass spectrometry (QTOF-MS) (Synapt G2, Waters) in negative ionization mode (cone voltage 3.0 kV) with an ACQUITY BEH C18 guard column (2.1 × 5 mm, 1.7 µm, Waters) coupled to an ACQUITY BEH C18 column (2.1 × 100 mm, 1.7 µm, Waters) and with collision energy of 6.0 eV. The faecal concentrations of SCFAs and BCFAs were determined using vendor software (Quanlynx, Waters).

Bacterial load in faeces was determined using approximately 500 µl of frozen faecal slurry (238–816 mg) and diluting it 400,000 times in physiological saline (8.5 g l−1 NaCl; VWR International). Next, 1 ml of the microbial cell suspension obtained was stained with 1 μl SYBR Green I (1:100 dilution in dimethylsulfoxide; shaded during 20 min incubation at 37 °C; 10,000 concentrate, Thermo Fisher Scientific). The flow cytometry analysis of the bacterial cells present in the suspension was performed using a Cytoflex flow cytometer (CytoFLEX 3; Beckman) as previously described (Supplementary Fig. 3)24. The final microbial load was calculated per gram of faeces.

Microbiome profiling

DNA was extracted in random order from the faecal slurries (n = 484) using DNeasy PowerLyzer PowerSoil kit (Qiagen, 12855-100), and the V3 region of the 16S rRNA gene was PCR amplified using 0.2 µl Phusion High-Fidelity DNA polymerase (ThermoFisher Scientific, F-553L), 4 µl high-fidelity buffer, 0.4 µl dNTP (10 mM of each base), 1 µM forward primer (primer binding upstream); 5′-A-adapter-TCAG-barcode-CCTACGGGAGGCAGCAG-3′), 1 µM reverse primer (primer binding reverse); 5′-trP1-adapter-ATTACCGCGGCTGCTGG-3′) and 0.05–5 ng faecal DNA in 20 µl total reaction volume. Both primers (TAG Copenhagen A/S) were linked to sequencing adaptors, and the forward primer additionally contained a unique 10 bp barcode (Ion Xpress Barcode Adapters) for each sample. The PCR program consisted of an initial denaturation for 30 s at 98 °C, followed by 24 cycles of 98 °C for 15 s and 72 °C for 30 s and a final extension at 72 °C for 5 min. The PCR products were purified by the HighPrep PCR clean-up system (AC-60500 Magbio) according to the manufacturer’s protocol. The resulting DNA concentrations were determined by Qubit HS assay and libraries constructed with mixing equimolar amounts of each PCR product. Partial 16S rRNA gene sequencing was performed on an Ion S5 System (ThermoFisher Scientific) using OneTouch 2 Ion5: 520/530 kit-OT2 400 bp and an Ion 520 Chip. The raw data were pre-processed into an amplicon sequence variant (ASV) table using our in-house pipeline56 based on the DADA2 algorithm and settings recommended for IonTorrent reads57, with taxonomy assigned to the ASVs using the Ribosomal Database Project (RDP, v18). The resulting ASV table, taxonomy and ASV sequences were merged into a phyloseq object for further analysis. For quantitative microbiome profiling analyses, the relative abundances derived from the pre-processed 16S rRNA sequencing analysis were adjusted for the bacterial loads as previously published58. In brief, samples with <10,000 reads were removed (n = 362) and downsized to even sampling depth, defined as the ratio between sample size (16S rRNA gene copy number corrected sequencing depth) and bacterial load. 16S rRNA gene copy numbers were retrieved from the rRNA operon copy number database rrnDB73 (ref. 59). The copy-number-corrected sequencing depth of each sample was rarefied to the level necessary to equate the minimum observed sampling depth in the cohort while assuring a minimum number of 10,000 reads in each sample and optimizing the chosen sampling depth to exclude as few samples as possible. In case of no copy number correction, an average copy number of 3.88 was used6.

Metabolic profilingPreparation of urine and faecal samples

Untargeted urine and faecal metabolomics were performed as previously published32. All urine samples were thawed on ice, centrifuged at 10,000 g at 4 °C for 2 min and transferred to a new tube to remove solid particles. The urine samples were kept cold on ice during preparation. Samples were randomized and pipetted into 15 plates (96-well). All urine samples from the same individual were placed on the same 96-well plate. Subsequently, they were diluted to 1:5 with an internal standard mixture (l-adenine-8-13C (Cambridge Isotope Lab), l-phenyl-d5-alanine-2,3,3-d3 (Cambridge Isotope Lab), caffeic acid 13C3 (Toronto Research Chemicals), caffeine 13C3 (Toronto Research Chemicals), l-tyrosine 13C9 (Sigma Aldrich), para-aminobenzoic acid (Sigma Aldrich), l-tryptophan-(indole-d5) (Sigma Aldrich), hippuric acid-[13C6] (IsoSciences), cortisone-d8 (Sigma Aldrich) and glycocholic acid-[2H4] (IsoSciences)). Quality control (QC) samples were obtained by mixing 20 µl of each urine sample in each plate (plate pools) and by mixing 20 µl of each plate pool to create the global pool. The QC samples, blank assays (0.1% formic acid) and mixtures of known standards (including 33 microbial-derived compounds) were included in each plate. The plates were sealed and stored at 4 °C until analysis (24 h maximum, otherwise stored at −80 °C). If the plate was frozen and thawed again before analysis, the plate was gently mixed by vortex stirring for 30 min immediately before analysis.

Faecal homogenates collected between day 2 and day 5 (n = 170) were thawed at room temperature for 30 min and vortexed. Approximately 50 ±5 mg (≈50 µl) of the homogenates were mixed with 96% ethanol and internal standard mixture (l-adenine-8-13C (Cambridge Isotope Lab), l-phenyl-d5-alanine-2,3,3-d3 (Cambridge Isotope Lab), caffeic acid 13C3 (Toronto Research Chemicals), caffeine 13C3 (Toronto Research Chemicals), l-tyrosine 13C9 (Sigma Aldrich), lysophosphatidylcholine (17:1d7) (Avanti Polar Lipids), l-tryptophan-(indole-d5) (Sigma Aldrich), hippuric acid-[13C6] (IsoSciences), cortisone-d8 (Sigma Aldrich) and glycocholic acid-[2H4] (IsoSciences)) resulting in a 1:60 dilution. The diluted samples were vortexed for 30 s and subsequently mixed at 60 °C for 2 min in a Thermo mixer at 1,400 r.p.m., before being centrifuged at 20,000 g (Eppendorf centrifuge 5417R) at 4 °C for 2 min. The supernatants were filtered through a 0.2 µm filter, and 200 µl of each faecal suspension was transferred to a 96-well plate, evaporated using a cooled vacuum centrifuge and re-dissolved in 200 µl 0.1% formic acid before the UPLC-MS. All faecal samples from the same individual were placed on the same 96-well plate, and QC samples were prepared in the same way as for the urine samples. In addition, each 96-well plate contained blank assays (96% ethanol) and mixtures of known standards (including 33 microbial-derived compounds).

UPLC-electrospray ionization-QTOF-MS analysis

Both urine and faecal samples were profiled by UPLC coupled with a QTOF mass spectrometer equipped with electrospray ionization (Synapt G2, Waters) in both positive and negative ionization modes32. Blank samples (0.1% formic acid), assay blanks, standard mixtures and QC samples were injected regularly to eval‎uate LC-MS system stability, possible contamination and/or loss of metabolites. The injected samples (5 µl) were separated on a reversed-phase column (ACQUITY HSS T3 C18 column, 2.1 × 100 mm, 1.8 µm) coupled with a pre-column (ACQUITY VanGuard HSS T3 C18 column, 2.1 × 5 mm, 1.8 µm). The mobile phases consisted of 0.1% formic acid in water (solvent A) and 0.1% formic acid in 70:30 acetonitrile/methanol (solvent B). The duration of the analytical run was 7 min with the following flow rate: start condition (0.5 ml min−1), 1 min (0.5 ml min−1), 2 min (0.6 ml min−1), 3 min (0.7 ml min−1), 4 min (0.8 ml min−1), 4.5 min (1.0 ml min−1), 6.4 min (1.1 ml min−1), 6.6 min (1.0 ml min−1), 6.8 min (0.5 ml min−1), 7.0 min (0.5 ml min−1), and the following gradient: start condition (5% B), 1 min (8% B), 2 min (15% B), 3 min (40% B), 4 min (70% B), 4.5 min (100% B), 6.6 min (5% B) and 7 min (5% B). Mass spectrometry data were acquired in full scan mode with a scan range of 50–1,000 mass/charge (m/z). Data-dependent acquisition was performed on the top three most abundant ions on QC samples (only urine) to provide MS2 data. Electrospray settings were the following: the cone voltage was 2.5 kV and 3.2 kV; the collision energy was 6.0 and 4.0 eV; and the temperature of the ion source and desolvation nitrogen gas temperature were 120 °C and 400 °C for positive and negative ionization mode, respectively.

Metabolomics data processing

The raw data obtained by UPLC-MS were converted to mzML format by publicly available msConvert (ProteoWizard Toolkit)60. The converted data were pre-processed using the open-source R package XCMS (v3.18) using the centWave algorithm (requiring three consecutive scans with an intensity of over ten counts)61. The pre-processing steps included noise filtering, peak picking, retention time alignment and feature grouping across samples, and filling of missing features, which were done separately for the urine and faecal samples (and for positive and negative mode), respectively. The detailed pre-processing parameter settings can be found in Supplementary Table 9. Noise filtering settings included that features should be detected in a minimum of 10% of all samples. Features with a retention time below 0.5 min or above 6.8 min were excluded. Data tables were generated comprising mass-to-charge ratio (m/z), retention time (rt) and intensity (peak area) for each variable in every sample. Each detected peak is represented by a feature defined by a rt and a m/z. The obtained data were corrected for within- and between-batch intensity drift using the locally estimated scatterplot smoothing correction method62. The processed data were normalized by the probabilistic quotient normalization63 method to correct for variations in urine and faecal concentrations within and between batches. Upon analyses of 15 plates with urine samples, QC samples clustered closely together in the principal component analysis score plots, confirm‎ing a stable UPLC system during the course of analysis with the exception of two plates in the negative mode and one plate in the positive mode, which had to be removed from further statistical analyses (Supplementary Fig. 4).

Moreover, features with high variability after normalization across the pooled QC samples were filtered out (coefficient of variation, CV% >50%). Finally, the CAMERA package64 (v1.52) was used to group features together based on retention time (tolerance = 0.1 s) and to annotate possible adducts and isotopes. Upon pre-processing, 641 and 651 molecular features were detected in the urine in positive and negative modes, respectively, whereas 453 and 445 molecular features were detected in faeces in positive and negative modes, respectively. MzMine 3 (ref. 65) and MassLynx (Waters) were used for data visualization.

Metabolite identification and structure elucidation

MS2 analyses were performed by an ultra-high performance LC system coupled to a Vion IMS QTOF mass spectrometer (Waters) for obtaining spectra with higher mass accuracy. The samples were separated on a reversed-phase column (ACQUITY HSS T3 C18 column, 2.1 × 100 mm, 1.8 µm) coupled with a pre-column (ACQUITY VanGuard HSS T3 C18 column, 2.1 × 5 mm, 1.8 µm) at a temperature of 50 °C. The mobile phases consisted of 0.1% formic acid in water (solvent A), methanol (solvent B), 0.1% formic acid in 70:30 acetonitrile/methanol (solvent C) and isopropanol (solvent D). The duration of the analytical run was 10 min with the following flow rate: start condition (0.4 ml min−1), 0.75 min (0.4 ml min−1), 6 min (0.5 ml min−1), 6.5 min (0.5 ml min−1), 8 min (0.6 ml min−1), 8.1 min (0.4 ml min−1), 9 min (0.4 mlmin−1), 10 min (0.4 ml min−1), and the following gradient: start condition (100% A), 0.75 min (100% A), 6 min (100% B), 6.5 min (70% C, 30% D), 8 min (70% C, 30% D), 8.1 min (70% C, 30% D), 9 min (100% A) and 10 min (100% A). Full scan acquisition was performed on selected urine samples with a scan range of 50–1,500 m/z. Data-dependent acquisition was performed on a selected list of precursors at three different collision dissociation energies, 10, 30 and 50 eV.

Mass spectra were manually interpreted, and metabolites were identified by matching the precursor ion and fragmentation patterns with databases such as Human Metabolome Database (https://hmdb.ca/), Metlin (https://metlin.scripps.edu/), mzCloud (https://www.mzcloud.org/), combinatorial database of bile acid conjugates66 (http://melolab.org/smilib/) and an in-house database (https://gitlab.com/stanstrup_R_packages/mscurate and https://gitlab.com/stanstrup_R_packages/xcms-annotator). In addition, we used several software annotations including GNPS (v30)67 (https://gnps.ucsd.edu/), microbeMASST (v2024.08.26)68 (https://masst.gnps2.org/microbemasst/) and SIRIUS (v6.0.5)69 (https://bio.informatik.uni-jena.de/software/sirius/), without obtaining additional plausible matches. Furthermore, authentic standards were run together with the samples with the highest intensity on the same batch and platform. If needed, the authentic standards were sulfated or glucuronidated with either biomimetic synthesis70 or chemical synthesis32. The identification level of metabolites that were identified was classified according to Sumner et al.71 as level I (confirmed by matching to a standard with two orthogonal measures (rt, m/z), level II (matching MS2 fragmentation to a spectral library), level III (compound classification) or level IV (unknown)25. See Supplementary Tables 7 and 8 for further details. 3-Hydroxy-2-oxindole, 5-hydroxyoxindole, 2-picolinic acid, 4-methylcatechol, xanthine, 2-oxindole-3-acetic acid, pantothenic acid, nicotinic acid, tryptophan, sebacic acid, pipecolic acid, glutaric acid, citric acid, psedouridine, taurine, 1,3-dimethyluric acid, suberic acid and 1,3,7-trimethyluric acid were purchased from Sigma-Aldrich. 4-Hydroxyhippuric acid, 1-methylxanthine and 1-methyluric acid were purchased from Toronto Research Chemicals.

Statistical analysis

Statistical analyses were conducted in R (v 4.2). The area under the curves for hydrogen and methane concentrations during the postprandial period was calculated using the trapezoid rule in GraphPad Prism (v 9.2.0). The normality of data was assessed with the Gaussian distribution and Shapiro–Wilk test procedure.

Mixed-effects linear regression models were used to examine the day-to-day fluctuations and inter-individual variation in gut environmental factors using data from all 9 days. The models were generated using the lme4 R package (v 1.1-31) as lmer (gut environmental factor ~ factor(day) + (1 | Participant ID)); moreover, ranova function from the lmerTest package (v 3.1-3) was used to perform the random effects-likelihood ratio tests to infer whether Participant ID significantly contributes to explaining the variation in the gut environmental factors. P value of < 0.05 was considered statistically significant. Coefficients of intra-individual variation were calculated as CVIntra = (s.d.intra/Meanintra) × 100 where mean and s.d. were based on all measurements from a single individual over the 9 days.

Gut microbiome beta-diversity analysis using Bray–Curtis distances as well as metabolome and diet beta-diversity analyses using Euclidian distances were performed with the phyloseq package (v 1.42.0) and PERMANOVA tests by adonis2 function from the vegan package (v 2.6) with 999 permutations and strata = Participant ID when testing the day-to-day fluctuations.

Single time point correlations were calculated using standard Spearman’s rank correlation, as implemented in the Hmisc R package (v 4.7), and heat maps were generated by the corrplot package (v0.92). Repeated measure correlations were performed using the rmcorr (v 0.5)72.

db-RDA was performed to quantify the effect sizes of gut environmental factors and other variables on the intra-individual and inter-individual variation in the gut microbiome (both relative and quantitative profiles at genus level) and faecal and urine metabolomes (untargeted data, all features). The analyses were performed with Bray–Curtis dissimilarity using the capscale function as implemented in the vegan package (v 2.6). With regards to intra-individual analyses, data available from all samples (day 1 to day 9) and strata = Participant ID were used. For the inter-individual analyses, data collected on day 2 (visit 1) were used separately for all participants (n = 61) and for the SmartPill subgroup (n = 50). The statistical significance was determined by permutation test with 9,999 random permutations (anova.cca function), and P values were adjusted for multiple testing by false discovery rate (Benjamin–Hochberg)73. An adjusted P value (q value) below 0.1 was considered significant.

For the untargeted metabolomics data, the area of each m/z feature was log-transformed, and missing values were imputed and replaced by values reflecting half of the minimum intensity of the given m/z feature. Linear regression models and SPLS models were performed to examine the relationship between the m/z features and the variables of interest (that is, segmental transit time and pH). The modelling was performed using the SmartPill-derived data and the 24 h postprandial urine metabolome collected at day 2 as well as the faecal metabolome closest to the time of the SmartPill egestion. The linear mixed models were performed with the lme4 R package (v 1.1-31). The multivariate SPLS models were performed with the caret R package (v 6.0-92). P values were corrected for multiple testing by the Benjamin–Hochberg false discovery rate (q value). Features with q < 0.1 were considered to be statistically significant, and only features selected by both the linear regression and SPLS were further submitted for identification including the MS2.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data availability

All sequencing data have been submitted to the National Center for Biotechnology Information Sequence Read Archive. BioProject ID, PRJNA1044006. MS2 data of global urine and faecal pool samples are deposited at MassIVE MSV000095466. Individual-level personally identifiable MS2 data from the participants cannot be made freely available to protect the privacy of the participants, in accordance with the Danish Data Protection Act and European Regulation 2016/679 of the European Parliament and of the council (GDPR) that prohibits open distribution even in pseudoanonymized form. Metabolomics data and data tables can be shared upon request. For data inquiries, please contact the principal investigator, H.M.R., via email. Access will be eval‎uated and granted upon signing a data processing agreement between the governing legal entities. Source data are provided with this paper.

Code availability

No custom code was generated for this work.

References

1. David, L. A. et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature 505, 559–563 (2014).

   CAS PubMed Google Scholar 

2. Wu, G. D. et al. Linking long-term dietary patterns with gut microbial enterotypes. Science 334, 105–108 (2011).

   CAS PubMed PubMed Central Google Scholar 

3. Johnson, A. J. et al. Daily sampling reveals personalized diet–microbiome associations in humans. Cell Host Microbe 25, 789–802.e5 (2019).

   CAS PubMed Google Scholar 

4. Guthrie, L. et al. Impact of a 7-day homogeneous diet on interpersonal variation in human gut microbiomes and metabolomes. Cell Host Microbe 30, 863–874 (2021).

   Google Scholar 

5. Falony, G. et al. Population-level analysis of gut microbiome variation. Science 352, 560–564 (2016).

   CAS PubMed Google Scholar 

6. Vandeputte, D. et al. Temporal variability in quantitative human gut microbiome profiles and implications for clinical research. Nat. Commun. 12, 6740 (2021).

   CAS PubMed PubMed Central Google Scholar 

7. Procházková, N. et al. Advancing human gut microbiota research by considering gut transit time. Gut https://doi.org/10.1136/gutjnl-2022-328166 (2022).

   Article PubMed Google Scholar 

8. Vandeputte, D. et al. Stool consistency is strongly associated with gut microbiota richness and composition, enterotypes and bacterial growth rates. Gut 65, 57–62 (2016).

   CAS PubMed Google Scholar 

9. Roager, H. M. et al. Colonic transit time is related to bacterial metabolism and mucosal turnover in the gut. Nat. Microbiol. 1, 16093 (2016).

   CAS PubMed Google Scholar 

10. Gill, P. A., van Zelm, M. C., Muir, J. G. & Gibson, P. R. Review article: short chain fatty acids as potential therapeutic agents in human gastrointestinal and inflammatory disorders. Aliment. Pharmacol. Ther. 48, 15–34 (2018).

    CAS PubMed Google Scholar 

11. Edamatsu, T., Fujieda, A., Ezawa, A. & Itoh, Y. Classification of five uremic solutes according to their effects on renal tubular cells. Int. J. Nephrol. 2014, 512178 (2014).

    PubMed PubMed Central Google Scholar 

12. O’Keefe, S. J. D. Diet, microorganisms and their metabolites, and colon cancer. Nat. Rev. Gastroenterol. Hepatol. 13, 691–706 (2016).

    PubMed PubMed Central Google Scholar 

13. Macfarlane, G. T., Gibson, G. R. & Cummings, J. H. C