Microbes on the Mind: Multi-modal Neuroimaging Reveals Gut-Brain Axes in Infants

기술 요약: 마음 위의 미생물: 다중 모드 신경 영상 기술을 통한 영아의 장-뇌 축 규명

문제 제기
장내 미생물 군집은 생후 첫 3년 동안 급격한 성숙 과정을 거치며, 이 시기는 뇌 발달의 결정적 시기와 일치한다. 미생물-장-뇌 축(microbiome-gut-brain axis)은 면역, 신경내분비, 미주신경 경로를 통한 상호 신호 전달을 포함하는 것으로 알려져 있으나, 초기 생애의 미생물 구성과 동시기에 발생하는 뇌 발달 사이의 구체적인 관계는 여전히 명확히 정의되지 않았다. 기존 연구들은 작은 표본 크기, 낮은 분류학적 해상도(예: 16S rRNA), 제한된 대사 특성 규명, 그리고 불완전한 신경 영상 범위 등의 한계를 지니고 있었다. 본 연구는 전형적으로 발달 중인 대규모 영아 코호트를 대상으로 장내 미생물과 다중 모드 뇌 영상 표현형 사이의 다변량 공변성(multivariate covariance)을 규명하는 것을 목표로 한다.

연구 방법론
본 연구는 UNC/UMN Baby Connectome Project의 가속 종단적 코호트를 활용하였으며, 여기에는 0~3세 사이의 전형적으로 발달 중인 아동 223명이 포함되었다. 데이터셋에는 각각 ±30일 이내에 수집된 분변 샷건 메타게놈(fecal shotgun metagenomics)과 다중 모드 MRI(구조, 확산, 휴지기 기능 MRI)의 489개 쌍 데이터가 포함되었다.

• 데이터 처리:

  • 미생물군: 샷건 메타게놈 데이터는 MetaPhl�lAn 3.0.2 및 HUMAnL을 사용하여 분류학적 프로파일(종, 속)과 기능적 프로파일(554개 MetaCyc 경로, 44개 장-뇌 모듈[GBM], 99개 장 대사 모듈[GMM])을 생성하였다.
  • 신경 영상: 구조적 MRI(T1w/T2w)는 iBEAT V2.0을 사용하여 148개 영역에 걸친 피질 형태(부피, 표면적, 두께)를 도출하였다. 확산 MRI는 34개 신경로에 걸친 미세구조 지표(FA, MD, MO)를 제공하였다. 자연 수면 중에 획득된 휴지기 fMRI는 7개의 표준 네트워크 내 기능적 연결성을 추출하도록 처리되었다.
  • 전처리: 신경 영상 표현형에서 연령 관련 추세를 제거하기 위해 참여자 무작위 절편을 포함한 일반화 가법 혼합 모델(GAMM)을 사용하였다. 미생물 특징은 중심 로그 비(CLR) 변환을 거쳤다.

• 통계 분석:

  • 미생물군 전 범위 연관 연구(MWAS): 선형 혼합 모델(MaAsLin2)을 사용하여 사이트, 연령, 성별, 출산 방식, 수유 방식 및 시퀀싱 깊이를 조정한 후 6가지 미생물 특징 세트와 신경 영상 측정치 간의 연관성을 테스트하였다.
  • 다변량 공변성: 107명의 완전한 데이터를 대상으로 리지 규제 정규화 정준 상관 분석(rCCA)을 적용하여 미생물 특징과 신경 영상 특징 사이의 주요 공변 패턴을 식별하였다.
  • 종단적 분석: 참여자 내 발달 기울기를 미생물 및 신경 영상 특징에 대해 추정하였다. 자기회귀 교차 지연 모델(autoregressive cross-lag models)을 통해 1~12개월 창(window) 동안의 정준 점수(canonical scores) 간 시간적 순서를 테스트하였다.
  • 검증: rCCA 패턴이 다른 다변량 목적 함수 하에서도 회복 가능한지 확인하기 위해 희소 정준 부분 최소 제곱법(Sparse canonical PLS)을 사용하였다. 참여자 수준의 재표본 추출을 통해 안정성을 평가하였다.

주요 결과

1. 발달 궤적:

   • 미생물군: 알파 다양성은 연령에 따라 증가하였다. 베타 다양성은 연령($R^2 = 0.107$)과 가장 강력하게 상관관계가 있었으며, 그 다음으로 어머니의 교육 수준, 수유 방식, 출산 방식 순이었다. Bifidobacterium 종들이 초기 프로파일을 지배하다가, 첫 1년 이후 Bacteroidetes와 Firmicutes로 변화하였다.
   • 신경 영상: 구조적 측정치는 첫 1년 동안 급격히 상승하였다. 기능적 네트워크는 뚜렷한 궤적을 보였다: 감각운동 연결성은 감소하였고, 시각 연결성은 초기에 정점에 도달한 후 감소하였으며, 고차 네트워크(디폴트 모드, 실행 제어)는 점진적으로 강화되었다.

2. 특징 수준 연관성 (MWAS):
   엄격한 유의성 임계값($q < 0.05$)을 충족하는 5가지 연관성이 발견되었다:

   • 총 뇌 부피는 퓨린 데옥시뉴클레오사이드 분해 및 퓨린 뉴클레오타이드 구제 경로와 양의 상관관계가 있었다.
   • Actinomyces sp. HMSC035G02의 풍부도는 피질 면적 및 부피 PC3 점수와 음의 상관관계가 있었다.
   • Gordonibacter pamelaeae의 풍부도는 실행 제어 기능적 연결성과 양의 상관관계가 있었다.

3. 다변량 공변 축 (rCCA):
   두 가지 주요 모드가 식별되었다:

   • 구조-대사 축 (CCA1): 피질 형태(특히 청각-시각 지역 대비를 나타내는 PC3: 더 큰 청각/측두엽 영역 vs 더 작은 시각 영역)를 미생물 대사 기능(멜리비오스 및 프럭탄 분해)과 연결하였다.
   • 기능-분류학 축 (CCA2): 종 구성(예: Bifidobacterium pseudocatenulatum) 및 GBM을 기능적 연결성, 특히 복측 주의 네트워크와 연결하였다.

4. 종단적 역학:

   • 기울기 공변성: 참여자 내 정준 점수의 성장률은 유의미한 상관관계를 보이지 않았다.
   • 시간적 순서: 자기회귀 교차 지연 모델 결과, 신경 영상 CCA1 점수가 13개월($\beta = 0.35$) 및 36개월($\beta = 0.31$) 후의 미생물 CCA1 점수를 예측하였다. 역방향 연관성은 보정 후 살아남지 못했다.

5. 방법론적 검증:
   희소 PLS 분석은 CCA1의 구조적 대비와 CCA2의 기능적 네트워크 패턴을 재현하였다. 그러나 별도의 테스트 상관관계는 완만하였는데, 이는 패턴이 샘플 내에서는 견고하지만 외부 분할로의 일반화는 제한적임을 시사한다.

의의 및 주장
저자들은 본 연구가 샷건 메타게놈과 포괄적인 신경 영상을 사용하여 전형적으로 발달하는 영아의 장-뇌 축을 최초로 다중 모드 및 종단적으로 규명했다고 주장한다. 연구 결과는 다음과 같다:

• 교차 영역 공변성: 특정 미생물 대사 및 분류학적 패턴이 뇌 구조(청각-시각 대비) 및 기능(복측 주의 연결성)의 뚜렷한 측면과 공변한다.
• 시간적 선행성: 측정된 창 내에서, 신경 영상 표현형(특히 구조-대사 축)이 이후의 미생물 상태를 예측하며, 이는 뇌 발달이 미생물 구성을 영향을 미치거나 공유된 미측정 요인이 이러한 순서를 유도할 수 있다는 시간적 순서를 시사한다.
• 방법론적 엄밀성: 본 연구는 단순한 단변량 상관관계를 넘어기 위해 엄격한 다변량 기법(rCCA, sparse PLS)과 종단적 모델링을 채택하였다.

본 논문은 이러한 결과가 조정된 변동과 시간적 순서를 나타내지만, 인과관계를 입증하지는 않는다고 명시적으로 밝히고 있다. 저자들은 관찰적 설계, 시변 혼란 변수(예: 질병, 항생제)의 부재, 그리고 메타게놈 데이터의 특성(활성 플럭스가 아닌 유전적 잠재력)이 기전적 결론을 내리는 것을 방해한다고 강조한다. 식별된 연관성은 방향성 있는 기전을 추론하기 위해 대사체학 및 실험 모델을 포함하는 향후 연구를 위한 토대로 기능한다.