Re: 장내 미생물 변동 요인 '대변통과시간, 대변 수분, pH) 2024년 nature review

[문형철]

대사체학 덕분에 인체의 엄청난 대사 결과를 볼 수 있는 논문이 나오기 시작!!

이 논문은 20-50대 젊은 사람 연구

동일한 식사를 할 때.. 장통과시간, pH, 대변 수분'이 

장내미생물 변동에 영향을 주더라!!

장내미생물이 분비하는 SCFAs 등의 변동치도 분석...

언젠가... 조만간에 혹은 이미 나왔는지도 모르겠는데...

케토제닉 식사, 지중해 식단이 ... 장내 미생물, 장내미생물 분비에 대한 논문이 나오고..

언젠가.. 건생병사 4과제 후 변화에 대한 논문을 쉽게 쓰는 날이 올듯!!

1. 개인 내(일일) 변동성
개인별로 장내 미생물 조성과 대사물질 농도가 매일 크게 변동.
변동은 주로 대변 수분 함량과 분변 pH 변화와 관련됨.

2. 개인 간 차이
장 통과 시간과 pH 수준의 차이가 개인 간 장내 미생물 조성 차이의 주요 원인.
통과 시간이 길고 pH가 높을수록 단백질 분해성 대사와 메탄 생성이 증가.

3. 대사물질 특징
탄수화물 발효로 생성된 **짧은 사슬 지방산(SCFA)**은 통과 시간이 짧고 pH가 낮을수록 증가.
반면, 단백질 분해 대사물질(branched-chain fatty acids, p-cresol 등)과 호기 메탄 농도는 통과 시간이 길고 pH가 높을 때 증가.

4. 식이 요인
커피, 섬유질 등 식이 요인이 장 통과 시간과 pH에 영향을 줌.
특히 식이 섬유 섭취량이 많을수록 단백질 분해성 유해 대사물질이 감소 → 장 건강에 보호 효과.

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

•  
•  
•  
•  

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와 커피, 식이, 호스트, 미생물 유래 대사산물 간의 연관성을 확인했습니다.

본 연구는

장 생리학과 환경이 인간 장 미생물 구성과 대사 다양성의 개인 차이를 이해하는 데

핵심적 역할을 한다는 것을 시사합니다.

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.

식사 후 ~4시간 위장 통과(gastric emptying time

            ~6시간 소장 통과(small intestine transit time)

            ~3일 대장 통과(large intestine transit time)

단쇄 지방산(SCFAs)은

당분 분해의 주요 미생물 대사산물로 일반적으로 유익한 것으로 간주됩니다10.

반면 미생물 단백질 분해는

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

건강에 해로운 대사물을 생성합니다11,12.

microbial proteolysis 

hydrogen sulfide

ammonia,

branched-chain fatty acids (BCFAs), 

p-cresol,

indole and phenylacetate

장내 미생물에 의한 단백질 발효는 대장 내 대사물 풀에 크게 기여하며, 호스트의 아미노산 균형에 영향을 미칠 수 있습니다. 그러나 장내 및 전신 순환계에서 단백질 분해 대사물의 역할에 대한 이해는 아직 제한적입니다. 최근 연구 결과와 이전 배양 기반 실험 결과를 종합한 검토는 미생물 단백질 발효가 장내 미생물군집을 변화시키고 다양한 생물활성 분자를 생성하는 데 중요한 역할을 한다는 것을 시사합니다. 이러한 대사산물은 장 내 염증 반응, 조직 투과성, 대장염의 심각도를 증가시키는 것으로 나타났습니다. 또한 비만, 당뇨병, 비알코올성 지방간 질환(NAFLD)과 같은 대사 질환의 발병과도 연관되어 있습니다. 단백질 분해 발효의 특정 제품인 수소 황화물, 암모니아, p-크레솔은 대장암 발병에 기여할 수 있습니다. 이러한 결과는 트립토판 대사물이 장 장벽 기능을 개선하고 다발성 경화증 모델에서 중증도를 완화한다는 다른 연구 결과와 충돌합니다. 장 내 단백질 분해 발효를 조사하는 추가 연구는 미생물과 호스트 대사 상호작용이 건강에 미치는 영향을 이해하는 데 핵심이 될 수 있습니다.

장 내 pH 변화는

장 미생물 구성과 대사 과정과도 연관되어 있습니다13.

SCFAs 및 기타 유기산은

대장 pH를 낮추며13,

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

젖산이 Escherichia coli O157:H7, Pseudomonas aeruginosa 및 Salmonella enterica serovar Typhimurium의 외막 투과성에 미치는 영향을 형광 탐침 흡수 분석법과 세균 용해에 대한 감수성 증가를 통해 연구하였다. 대조 실험을 위해, 킬레이트 작용을 통해 투과성을 증가시키는 EDTA와 염산(lactic acid가 생성하는 pH 값과 동일한 pH 값에서)을 사용한 유사한 실험이 수행되었으며, KCN이 존재하는 조건에서도 실험이 진행되었습니다. 이미 5 mM(pH 4.0) 젖산은 각 종에서 현저한 투과성을 유발했으며, 형광 분석에서의 효과는 EDTA나 HCl보다 강했습니다. KCN 존재 시에도 유사한 결과가 관찰되었으나, P. aeruginosa의 경우 KCN 존재 시 HCl의 효과가 증가하는 현상이 관찰되었습니다. 젖산과 염산에 의한 투과성 증가는 MgCl2에 의해 부분적으로 억제되었습니다. 젖산은 E. coli와 세로형 Typhimurium을 나트륨 도데실 설페이트(SDS)의 용해 작용에 대해 HCl보다 더 효율적으로 감작시켰으며, 두 산 모두 P. aeruginosa를 SDS와 트리톤 X-100에 대해 감작시켰습니다. P. aeruginosa는 젖산과 HCl에 의해 리소자임에 대한 감수성이 효과적으로 증가했습니다. 이 산들은 세로형 Typhimurium에서 상당량의 리포폴리사카라이드를 방출시켰으며, 전기영동 및 지방산 분석을 통해 방출된 물질의 분석 결과, 젖산은 EDTA나 HCl보다 외막에서 리포폴리사카라이드를 방출시키는 데 더 활성적임을 보여주었습니다. 따라서 젖산은 pH를 낮추어 항균 효과를 발휘하는 것 외에도 그람 음성 세균의 외막을 투과성으로 만드는 역할을 하며, 다른 항균 물질의 효과를 증강시키는 작용을 할 수 있습니다.

그러나

생리적 요인(통과 시간, 장 내 pH 등)에 의해 결정되는 장 환경이

식이-호스트-미생물 대사 과정과 어떻게 연관되는지에 대한 연구는 아직 부족합니다.

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

미래 맞춤형 식이 미생물군 기반 전략 개발에 결정적 역할을 할 수 있습니다.

통과 시간은 무선 운동 캡슐(SmartPills)을 사용하여 측정할 수 있습니다. 이 캡슐은 위장관 전체의 내강 pH, 온도 및 압력을 측정하여 세분화된 통과 시간을 산정할 수 있습니다. SmartPills는 방사선 불투과성 표지자 및 스캔티그래피와 함께 통과 시간을 직접 측정하는 표준 임상 방법 중 하나이며15, 이전 미생물군 연구에서도 사용되었습니다16,17.

본 연구에서는 61명의 건강한 자원자를 대상으로 9일간의 관찰 연구를 수행했으며,

SmartPills를 사용하여

전체 장 및 세분화 장 통과 시간과 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²; 표 1 및 확장 데이터 그림 1)을 모집하고,

9일 동안 평소 생활 방식과 식습관을 유지하도록 요청했습니다(그림 1a).

연구는 두 차례의 방문(2일차와 9일차)으로 구성되었으며, 이 기간 동안 공복 혈당, 인슐린, C-펩타이드, 호흡 수소 및 메탄을 측정했습니다(표 1). 첫 번째 방문 시 참가자들은 일일 에너지 필요량의 25%를 차지하는 아침 식사(보충 자료 표 1)를 제공받았으며, 자원자 중 일부(n = 50)는 전체 장 및 부위별 통과 시간과 pH를 측정하기 위해 무선 운동 캡슐(SmartPill)을 복용했습니다18. 이전 연구에서는 모니터링 전 그라놀라 바(SmartBar)를 사용했습니다19,20, 우리는 최근 연구21와 유사한 복합 식사를 사용하여 식이-미생물군 상호작용을 조사했습니다. 식후 호흡 및 소변 샘플은 그림 1a에 표시된 대로 수집되었습니다. 참가자들은 myfood24 영양 플랫폼(https://www.myfood24.org)을 사용하여 1일부터 8일까지 매일 24시간 식이 기록을 작성했으며, 배변 습관(배변 시간, 브리스톨 배변 형태 척도(BSS)22로 평가된 배변 일관성, 배변 빈도(하루 배변 횟수))을 기록하고, 매일 소변(첫 아침 샘플)과 대변 샘플(첫 배변 시)을 수집했습니다. 연구 대상자는 정상적인 배변 습관(표 1)을 보였으며, BSS 중간값은 4형이었고 하루 배변 횟수는 1회였습니다. 3일차와 5일차에 자기 보고형 스위트콘 통과 시간 검사(corn TT)를 통해 통과 시간을 추정했습니다. 모든 수집된 대변 샘플(n = 484)에서 대변 수분 함량(대변 습도의 지표로, 통과 시간의 대리 지표로 사용됨23), pH 및 미생물 부하를 측정했습니다. 모든 수집된 소변 샘플(일일 채취 및 식후 샘플, n = 1,154)과 대변 샘플의 일부(n = 170)는 비표적 액체 크로마토그래피-질량 분석법(LC-MS) 대사체학 분석을 통해 소변 및 대변 대사체 프로파일을 얻었습니다. 마지막으로, 대변 샘플의 일부(n = 362)에 대한 16S 리보솜 RNA(rRNA) 유전자 시퀀싱을 통해 장 미생물군 구성체를 얻었으며, 미생물 부하를 조정24한 후 상대적 미생물군 프로파일(RMPs)과 정량적 미생물군 프로파일(QMPs)을 평가했습니다.

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).

분변 수분 함량과 pH가 장 미생물군집의 일일 변동을 설명합니다

개체 내 대사체와 미생물군집의 변동을 유발하는 요인을 탐색하기 위해 거리 기반 중복 분석(db-RDA)을 수행했습니다. 일일 식이 대영양소와 식이 섬유, 장 환경 요인을 고려했습니다. 식이 성분 중 어느 것도 장 미생물군집이나 대사체의 개체 내 변동을 설명하지 않았습니다. 반면, 대변 수분, 대변 pH, BSS 및 배변 시간은 장 미생물군집에 유의미한 영향을 미쳤습니다(QMP, 속 수준; 그림 2a), 각각 변동의 3.5%, 2.5%, 2% 및 1.3%를 설명했습니다. RMP 데이터(확장 데이터 그림 3a) 및 이전 연구6,25에서도 유사한 결과가 관찰되었습니다. 특히, 이러한 설명 변수들은 장 통과 시간의 대리 지표로 작용하며, 이는 장 통과 시간의 일일 변동이 장 미생물군집 변동에 반영된다는 것을 시사합니다.

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.

변의 수분 함량과 분변 pH는

각각 3.1%와 3%를 차지하며,

소변 대사체 내 개인 내 변이성을 설명했습니다.

이는 하루하루의 미묘한 변동에도 불구하고(그림 2b)

이러한 요인들이 영향을 미쳤음을 시사합니다.

이는 대장 내 수분 함량과 pH의 작은 변화가

호스트-미생물군집 대사 과정과 연관될 수 있음을 의미합니다.

그러나 이러한 관찰 결과는 일일 식이 변동에도 영향을 받을 수 있습니다. 장 환경 요인은 분변 대사체 내 개인 내 변동성에 기여하지 않았습니다(확장 데이터 그림 3b). 주의할 점은 분변 대사체 데이터가 3일 연속 수집된 데이터에서만 도출되었으며, 분변 수분 함량은 여전히 유의미한 영향을 미칠 경향(P = 0.081)을 보였습니다.

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.

통과 시간과 pH는 개인 간에 차이가 있습니다

SmartPills를 사용하여 전체 장 통과 시간(WGTT), 위 배출 시간(GET), 소장 통과 시간(SBT), 대장 통과 시간(CTT), 장 통과 시간(ITT; SBT + CTT) 및 위장관(GIT) 전체의 pH를 측정했습니다(확장 데이터 그림 3c). 8명의 개인에서 캡슐이 위에서 8시간 이상 머무는 현상이 관찰되었으며, 이는 다른 연구에서도 보고된 일반적인 현상입니다19,26. 따라서 이 참가자들의 GET 및 WGTT 값은 분석에서 제외되었습니다. 또한, 신호 손실로 인해 1명의 참가자에서 CTT 및 WGTT를 측정할 수 없었습니다.

통과 시간의 중간값은 다음과 같습니다: GET, 4.8시간 (범위 3.1–6.2시간); WGTT, 23. 3시간 (12.4–72.3시간); CTT, 13.6시간 (2.1–63.5시간); SBT, 5.1시간 (2.5–10.3시간)으로, 건강한 인구 집단에 대한 이전 보고 데이터와 일치했습니다27. 비교를 위해 옥수수 TT는 3일차에 중간값 23.6시간 (10.8–109.7시간)과 5일차에 19.7시간 (12.0–84.5시간)을 나타냈습니다. 또한 두 옥수수 TT 측정값 사이에 강한 상관관계가 관찰되었습니다(스피어먼 상관 계수(SCC) = 0.8, P < 0.001), 이는 개인 내 일관성을 시사합니다. 두 날짜에 걸친 옥수수 TT의 평균 중간값은 21.7시간 (11.7–97.1시간) (그림 2c)로, SmartPill로 측정된 WGTT와 유사했습니다. 그러나 WGTT와 옥수수 TT 사이에는 상관관계가 관찰되지 않았습니다(확장 데이터 그림 4a), 이는 평균적으로는 유사한 결과를 제공하지만, 개별적으로는 두 방법이 다른 결과를 보였음을 나타냅니다.

세분화 위장관 통과 시간, 옥수수 TT, 위장 환경 요인 및 참가자 특성 간의 관계를 탐구할 때(확장 데이터 그림 4a), 두 방법에서 기록된 통과 시간이 BSS와 음의 상관관계를 보였으며, 이는 이전 연구에서도 보고된 바 있습니다23,28. 또한 여성은 남성보다 CTT가 유의미하게 더 길었으며, 여성 내에서 월경 상태의 영향은 관찰되지 않았습니다(확장 데이터 그림 4b).

위장관 세분화 구간별 pH의 큰 개인 간 변동성도 관찰되었습니다(그림 2d). 상부 위장관(GIT)의 중간 pH 값은 위(0.9, 범위 0.5–4.9), 십이지장(6.1, 5.0–7.2), 소장(7.4, 6.4–8.2) 순이었습니다. 근위 결장의 pH는 약간 산성(6.3, 5.3–7.0)이었으며, 원위 결장(6.9, 5.0–8.2)과 S자 결장(7.2, 5.6–8.6)으로 갈수록 점차 증가했습니다. 흥미롭게도, S자 결장에서 직장까지 pH가 약간 감소했으며(7.0, 5.7–8.6), 대변 pH(6.9, 6.6–7.3)에서도 감소가 관찰되어 직장에 진입한 후 산성화 과정이 발생함을 나타냈습니다.

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.

CTT와 pH는 장 미생물군집 변이에 기여합니다

참가자 특성, 임상 변수, 식이 요인 및 장 환경 및 생리적 요인이 장 미생물군집과 대사체군집의 개인 간 변이를 설명하는 정도를 정량화하기 위해, 모든 참가자(n = 61; 보충 표 4)의 2일차 분변 및 24시간 소변 수집 데이터로부터 db-RDA를 수행했습니다. 또한, 2일차에 SmartPills에서 수집된 전체 장 및 세분화 장 통과 시간과 pH를 사용하여 동일한 분석을 수행했습니다(n = 50). 대변 수분 함량과 원위 결장 pH는 QMP의 개인 간 변이에 관련된 주요 요인으로, 2일차 및 기타 날짜에서 각각 변이의 5.5%와 5%를 설명했습니다(보완 표 4). 이전 대규모 코호트 연구 결과와 달리5, BSS는 이 인구 집단에서 QMP 변이의 유의미한 부분을 설명하지 않았습니다. WGTT, CTT, 옥수수 TT 및 대변 pH는 24시간 소변 대사체 변이의 9.1%, 6.2%, 4.9% 및 5.4%를 각각 설명했으며, 이는 연령이 변이의 4.5%를 설명한 것과 대비됩니다(그림 2f). 이 기여도는 다른 날의 소변 대사체와 비교할 때도 일관되었습니다(보충 표 4). 반면, 분절 통과 시간은 분변 대사체의 개인 간 변이에 기여하지 않았으며, 원위 결장의 pH와 식이 섬유 섭취량이 각각 변이의 6.8%와 5.9%를 설명하는 가장 큰 영향을 보였으나, 다중 검정 조정 후 유의미하지 않았습니다(보충 표 4).

여성 참가자(월경 중이 아닌 그룹, n = 30; 월경 중 그룹, n = 13)를 대상으로 연구 기간 중 월경의 영향을 분석한 결과, 장 미생물군집에 3.7%, 소변 대사체에 3.6%의 효과 크기가 관찰되었으나 통계적 유의미성은 없었습니다. 연령이 소변 대사체에 미치는 유의미한 효과 크기 및 두 여성 그룹 간의 연령 차이(P = 0.01)를 고려할 때, 연령이 이러한 관찰된 효과에 기여할 수 있습니다. 본 연구 결과는 개인의 장 환경이 장 미생물군과 소변 대사 프로파일의 개인 간 차이에 크게 기여함을 강조합니다.

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.

개인의 장 미생물군과 대사체 프로파일은 동적입니다

다음으로,

호흡 수소 및 메탄,

분변 SCFAs(아세테이트, 프로피오네이트, 부티레이트, 발레레이트 및 카프로에이트) 및

BCFAs(이소부티레이트, 이소발레레이트 및 2-메틸부티레이트; 표 1)을 포함한 미생물 유래 대사물,

그리고 분변과 소변에서 검출된 16개의 다른 미생물 유래 대사물

(단백질 분해 표지물질인 p-크레졸 황산염(PCS), 페닐아세틸글루타민(PAGln) 및 인돌황산염)을 평가했습니다.

상당한 일일 변동성이

관찰되었습니다(그림 3a 및 보충 그림 2).

호흡 메탄과 수소의 CVIntra 중간값은 각각 141%와 47%였으며,

두 가스 모두 두 시점 간에 중간 정도의 양의 상관관계가 관찰되었습니다

(수소, SCC = 0.42, P < 0.001; 메탄, SCC = 0.66, P < 0.001).

분변 내 SCFA 및 BCFA 농도는

하루하루 크게 변동했으며(중앙값 CVIntra 26%~40%),

발레이트가 가장 적게 변동하고

아세테이트가 가장 크게 변동했습니다.

同様に, 단백질 분해 표지물질의 상대적 풍부도는 날마다 크게 변동했으며,

PAGln, 인돌설페이트, PCS의 중간 CVIntra는 각각 26%, 42%, 39%였습니다.

이러한 결과는

습관적인 식이에서 호흡, 분변, 소변에 존재하는 미생물 유래 대사물이

날마다 변동함을 시사합니다.

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.

다음으로, 반복 측정 분석(그림 3b, 확장 데이터 그림 5 및 보충 표 5)과 스피어맨 상관 분석(그림 3c)을 사용하여 특정 미생물 대사산물과 장 환경 및 식이 요인 간의 연관성을 조사했습니다. 분변 SCFAs는 분변 pH와 음의 상관관계를 보였으며, 특히 부티레이트가 가장 강한 상관관계를 나타냈습니다(r = −0.77, q < 0.001), 이는 이전 인간 연구29와 일치합니다. 또한, 분변 프로피오네이트 수치(SCC = −0.25, q < 0.1)가 높을수록 CTT가 짧았으며, 분변 부티레이트에서도 유사한 경향이 관찰되었습니다(SCC = −0.29, P < 0.05, q = 0.3). 분변 부티레이트는 직장 pH와 음의 상관관계를 보였으나(SCC = −0.37, P < 0.05, q = 0.2), 대장의 다른 부위 pH와는 상관관계가 없었습니다. 이는 부티레이트 생산이 직장 및 분변에서 관찰된 pH 감소에 기여할 수 있음을 시사합니다.

단백질 분해 지표인 요중 PAGln과 분변 BCFAs는 분변 수분과 음의 상관관계를 보였으며, 요중 PCS는 분변 pH와 양의 상관관계를 보였습니다(PAGln, r = −0.12, q < 0.1; 이소부티레이트, r = −0.39, q < 0. 05; 이소발레이트, r = −0.37, q < 0.1; 2-메틸부티레이트, r = −0.43, q < 0.05; PCS, r = 0.12, q < 0.1). 同様に, 소변 내 PCS 수치가 높을수록 CTT 및 ITT가 길어지는 경향이 관찰되었으며 (SCC = 0.48, SCC = 0.44, 각각, q < 0.05), PAGln 및 인돌-젖산에서도 유사한 경향이 관찰되었습니다. 또한, 호흡 메탄은 대변 수분 함량 감소 및 CTT 연장과의 연관성이 확인되었습니다. 참고로, 이러한 결과는 SBT에서는 관찰되지 않았으며, 이는 CTT가 이러한 대사물의 풍부도를 결정하며, 대장 통과 시간이 길수록 미생물 단백질 분해와 연관될 수 있음을 시사합니다. 이는 사카로리틱 발효의 기질 고갈 때문일 수 있습니다9,30. 요중 인독실-글루쿠로니드는 대장 말단부의 pH와 양의 상관관계를 보였습니다(SCC = 0.33, q < 0.1), 그리고 요중 PCS와 직장 pH 사이에도 유사한 경향이 관찰되었습니다(SCC = 0.32, P < 0.05). 이 대사산물들은 소장 및 근위 대장의 pH와 상관관계를 보이지 않았으며, 이는 원위 장에서 미생물 단백질 분해가 근위 장에 비해 pH에 더 큰 기여를 한다는 것을 시사합니다. 요약하면, CTT와 대장 pH는 호흡, 분변 및 소변의 여러 미생물 대사산물 수준과 연관되어 있지만, SBT와 소장 pH는 그렇지 않습니다. 또한, 미생물 대사산물과 식이 성분 사이의 여러 연관성을 발견했으며(보충 표 5), 식이 섬유 섭취량과 대변 BCFAs, 요중 PCS, 요중 PAGln 사이의 역상관 관계가 눈에 띄었습니다.

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).

대변과 요 대사체는 통과 시간과 pH와 연관되어 있습니다

장 생리학과 관련된 알려지지 않은 대사적 특성을 탐구하기 위해, 우리는 소변과 분변의 대사체를 프로파일링하기 위해 비표적 대사체학을 사용했습니다. 소변과 분변에서 식별된 모든 분자적 특성에 대해 단변량 및 다변량 통계 모델을 적용했습니다. 먼저 SmartPill에서 파생된 데이터와 식사 후 24시간 동안 수집된 소변 대사체, SmartPill 배출 시점에 가장 가까운 시점에서 수집된 분변 대사체에 대해 희소 부분 최소 제곱(SPLS) 모델을 적용했습니다. 동일한 데이터에 선형 회귀 모델을 수행한 후 두 모델에서 선택된 특징(446개 고유 특징; 보충 표 6)을 추가로 조사했습니다.

소변과 분변의 여러 대사 특징이 전체 장 및 부위별 통과 시간과 pH와 연관되었습니다(그림 4a,b). 이러한 특징을 더욱 자세히 조사하기 위해 해당 샘플을 tandem mass spectrometry (MS2)로 분석하고, 가능한 경우 표준 물질과 일치시켜 33개의 대사체를 식별했습니다(보충 표 7 및 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.

PCS 및 PAGln의 요중 농도 외에도,

장 내 미생물에 의해 분해되는 아로마틱 아미노산 트립토판과 티로신에서 유래한 여러 요중 대사산물이

장 통과 시간과 분변 pH와 연관되어 있음을 발견했습니다.

구체적으로, 5-하이드록시-2-옥신돌 황산염, 3-하이드록시-2-옥신돌 황산염 및 4-하이드록시벤조산 황산염은 더 긴 WGTT/CTT와 연관되었으며, 3-하이드록시-2-옥신돌 글루쿠로니드는 분변 pH가 높은 것과 관련되었습니다. 반면, 분변 트립토판은 분변 pH와 음의 상관관계를 보였습니다. 또한, 분변 프로린과 요중 피콜리노일글리신 수치는 각각 분변 및 직장 pH 증가와 연관되었습니다.

분변 내 여러 이산화탄소산인 피멜릭산, 수베릭산 및 세바시산은 WGTT 및 CTT와 양의 상관관계를 보였습니다. 반면, 분변 글루타르산과 피페콜릭산은 WGTT/CTT 및 시그모이드, 직장 및 분변 pH와 음의 상관관계를 보였습니다. 피페콜릭산은 식물에서 풍부하게 존재하지만, 장 미생물에 의해 라이신에서 생성될 수 있습니다31. 또한, 요중 시트르산 수치가 높을수록 근위 결장의 pH와 양의 상관관계를 보였습니다.

또한, 이전에 뉴 노르딕 다이어트와 지중해 다이어트와 연관된 것으로 알려진 분변 2-옥신돌-3-아세트산32,33은 WGTT, CTT 및 분변 pH와 음의 상관관계를 보였습니다. 마찬가지로, 분변 판토텐산과 니코틴산은 각각 CTT 및 분변 pH와 음의 상관관계를 보였습니다. 또한, 소변 내 디히드로페룰산 글루쿠로니드와 아르기닌산은 직장 pH와 음의 상관관계를 보였으며, 대변 내 p-하이드록시페닐락틱산은 대변 pH와 음의 상관관계를 보였습니다.

4-하이드록시히푸르산과 커피 섭취를 나타내는 소변 내 여러 지표(1-메틸우르산, 1-메틸잔틴, 1,3-디메틸우르산, 1,7-디메틸우르산, 1,3,9-트리메틸우르산)는 소장 pH와 음의 상관관계를 보였습니다. 분변 내 1-메틸잔틴과 1,3,9-트리메틸우르산은 WGTT 또는 분변 pH와 음의 상관관계를 보였으며, 이는 커피 섭취와 장 기능 간의 연관성을 시사합니다. 또한, 직장 pH와 식물성 식품에 함유된 케르세틴의 대사산물인 요중 4-메틸카테콜 설페이트 사이에 양의 상관관계가 관찰되었습니다34. 요중 타우린과 분변 콜리산은 소장 pH와 양의 상관관계를 보였으며, 이는 담즙산이 위에서 오는 산성 chyme을 중화하는 역할을 한다는 것을 지지합니다35.

마지막으로, RNA의 주요 구성 성분인 요중 가짜우리딘은 CTT 및 S자 결장 pH와 역상관 관계를 보였으며, 이는 이전 연구 결과와 일치합니다9. 가짜우리딘은 분변에서도 검출되었으며, 분변 pH와 유사한 역상관 관계를 보였으며, 데옥시-잔토신과 잔토신도 마찬가지였습니다. 이는 세포 회전율 증가와 결장 pH 저하 사이의 연관성을 시사합니다.

전체적으로, 비표적 LC-MS 대사체학을 통해 대장 후부에서 WGTT, CTT 및 pH와 연관된 호스트, 미생물, 식품 유래 대사체들을 식별했으며, 이는 식이, 장 환경, 호스트, 미생물군집 간의 상호작용을 강조합니다.

미생물 알파 다양성은 장 통과 시간과 연관됩니다

식별된 대사산물과 장 미생물군집 간의 잠재적 연관성을 탐색하기 위해 스피어맨 상관 분석을 수행했습니다(그림 5). 미생물 알파 다양성 지표와 미생물 단백질 분해, CTT 및 ITT 사이에 강한 양의 상관관계가 발견되었습니다. 반면, 알파 다양성은 대변 수분 함량과 미생물 당분해와 음의 상관관계를 보였습니다.

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.

미생물 단백질 분해산물과 이산화탄소산은 Intestimonas, Flavonifractor, Eubacterium, Lachnospira, Clostridium, Oscillibacter, Alistipes, Dialister 및 Akkermansia를 포함한 여러 세균 속의 절대 풍부도와 양의 상관관계를 나타냈습니다.

같은 속들은 분변 내 트립토판, 옥시인돌-3-아세트산 및 다양한 커피 유래 대사산물 수준과 음의 상관관계를 보였습니다.

예상대로, 이 속들은 ITT 및 CTT가 더 길고 분변 pH가 높은 것과 양의 상관관계를 보였으며, 분변 수분 함량 및/또는 BSS와 음의 상관관계를 보였습니다. 반면, SCFAs를 생성하는 속인 Agathobacter, Faecalibacterium 및 Blautia36,37과 젖산을 생성하는 Streptococcus는 모두 분변 니코틴산, 판토텐산 및 커피 유래 대사산물과 양의 상관관계를 보였습니다. 특히, Oscillibacter, Alistipes 및 Akkermansia는 장 통과 시간과/또는 변비와 관련된 샘플에서 반복적으로 증가된 것으로 확인되었습니다8,9,28, 반면 부티레이트 생성 속인 Faecalibacterium 및 Agathobacter는 더 짧은 장 통과 시간과 연관되었습니다28,38. 요약하면, 이러한 관찰 결과는 장 내 세균, 대사산물 및 장 생리학 간의 상호 의존성을 강조합니다.

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.

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.

논의

장 통과 시간과 pH는

장 미생물군집의 구성과 대사 과정에 중요한 요인입니다7.

본 연구에서는

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

이러한 변이는

미생물군집 구성과 호스트-미생물군집 공대사 차이를 설명했습니다.

pH와 통과 시간은 미생물 성장과 효소 활성에 영향을 미치기 때문에39,

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

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

최근 연구는 GIT를 따라 미생물군집과 대사체 구성에 차이가 있음을 확인했습니다40.

GIT를 따라 샘플링과 지역별 pH 및 통과 시간 측정을 결합한 미래 연구가 필요하며,

이를 통해 이러한 관계를 최종적으로 규명할 수 있을 것입니다.

본 연구는 개인별 장 내 pH 차이로 인해 pH에 민감한 섭취 가능 장치40,41를 활용한 연구에 도전 과제가 될 수 있음을 강조합니다. 인간 연구에서 미생물 유래 대사물의 상당한 일일 변동성을 관찰했기 때문에 반복 측정이 필요하다는 점도 주목할 만합니다.

우리는 장 통과 시간이 길과 연관된 여러 대사물을 식별했으며, 이는 변비 환자 그룹에서 증가된 것으로 보고된 대사물입니다42,43,44. 이 중 디카르복실산(pimelic acid 포함)은 더 긴 ITT 및/또는 높은 pH와 연관되었으며, 변비와 일관되게 연관된 세균과도 관련되었습니다8,9,28. pimelic acid는 미생물 지방산 대사에서 유래했을 가능성이 있으며45,46, 만성 신장 질환 환자에서 변비와 자주 연관된 높은 분변 수준에서 발견되었습니다47,48. 최근 연구에서는 장의 말단부에서 이산화탄소산 함량이 증가하는 것이 관찰되었으며40, 저자들은 이는 호스트 상피 세포의 분해 때문일 수 있다고 추측했습니다. ITT가 길어지는 것이 상피 세포의 회전율과 탈락과 연관될 수 있는지 추가 연구가 필요합니다.

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

제한된 코호트 규모에도 불구하고, 본 연구는 건강한 인구에서 장 세분화 통과 시간과 pH가 장 미생물군 구성 및 대사와의 유의미한 연관성을 보여주었습니다. 고려해야 할 잠재적 한계는 SmartPill 측정 전 아침 식사 선택입니다. 과거 검증 연구 및 기준 데이터는 특정 식사/영양소 조합에 의존하기 때문에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. 1. 1. 1.