Paired plasma cfRNA and fecal metagenomics jointly resolve the two-compartment bacterial landscape of SCZ

기술 요약: 쌍을 이룬 혈장 cfRNA와 분변 메타게놈을 통한 조현병(SCZ)의 2개 구획 박테리아 경관 공동 규명

문제 제기
조현병(SCZ)은 유전적 감수성이 미생물군집을 포함한 환경적 요인과 상호작용하는 이질적인 질환이다. 장내 미생물군이 '제2의 게놈'으로서 미생물-장-뇌 축에 영향을 미친다는 점은 잘 알려져 있으나, 기존 연구들은 코호트의 이질성, 혼란 변수(식단, 약물), 생물정보학적 불일치로 인해 재현성 측면에서 상당한 어려움을 겪고 있다. 또한, 분변 메타게놈은 장내 미생물 저장소를 특징화하지만, 이는 실제 숙체의 전신 생리에 관여하는 분자 신호를 전달하는 간접적인 대리 지표 역할에 그친다. 장에서 방출되어 장벽을 통과하고 순환계에서 전신적으로 검출 가능한 특정 미생물 분자 집단은 정신 의학 연구에서 여전히 미개척 영역으로 남아 있다. 따라서 잔류 장내 미생물군과 숙체에 노출된 순환 미생물 신호를 구분하여, 조현병에서 나타나는 전신적 미생물 생태계 재편을 이해할 수 있는 기준이 절실히 필요하다.

연구 방법론
본 단면 연구는 동일한 부모 코호트에서 모집된 178명의 참가자(조현병 환자 89명, 건강 대조군 89명)를 대상으로 쌍을 이룬 혈장 세포 유리 RNA(cfRNA) 시퀀싱과 분변 샷건 메타게놈 분석을 통합하였다.

• 시료 처리 및 시퀀싱: 엄격한 공복 상태에서 혈장을 채취하였으며, 혈장 cfRNA를 추출하여 샘플당 평균 약 11.38 Gb의 깊이로 시퀀싱하였다. 분변 DNA를 추출하여 샷건 메타게놈 시퀀싱을 수행하였다.
• 오염 제어: 엄격한 오염 인지 파이프라인을 구현하였다. 여기에는 추출 블랭크(extraction blanks), 피부 스왑, 순수 배양액 등 27개의 프로세스 컨트롤이 생물학적 샘플과 함께 처리되었다. 두 가지 직교 분류기(MetaPhlAn 및 Kraken2)를 통해 특정 풍부도 및 통계적 농축 기준을 충족하며 추출 블랭크에서 검출된 종은 고신뢰 오염원으로 지정되었다. 총 47종의 오염 종이 혈장 매트릭스에서 제거되었다.
• 생물정보학적 분석:
  • 분류학적 프로파일링: 1차 분류 할당을 위해 MetaPhlAn v4를 사용하였으며, Kraken2를 사용하여 직교 검증 및 오염 평가를 수행하였다.
  • 기능적 프로파일링: 혈장 cfRNA 리드로부터 효소(EC) 및 오솔로지(KO) 풍부도를 추론하기 위해 HUMAnN v4를 사용하였으며, 오염이 조정된 풍부도 계산법을 적용하였다.
  • 통계 모델링: 연령, 성별, BMI를 보정하여 분석하였다. 차별적 풍부도는 MaAsLin2를 사용하여 평가하였다. 알파 및 베타 다양성은 Shannon/Simpson 지수와 Bray-Curtis 비유사도(PERMANOVA)를 통해 평가하였다.
  • 구획 간 분석: 쌍을 이룬 혈장 및 분변 프로파일 간의 Spearman 상관관계를 계산하였다. 그룹 수준의 일치성과 개별 수준의 쌍 형성을 구분하기 위해 10,000회의 순열 테스트를 수행하였다.
  • 임상적 연관성: 미생물 특징과 PANSS(양성 및 음성 증후군 척도) 하위 척도 점수 간의 연관성을 강건한 표준 오차를 이용한 선형 회귀를 통해 테스트하였다.
  • 분류: 진단 성능(AUC)을 평가하기 위해 반복 5회 교차 검증을 사용하여 혈장 전용, 분변 전용, 결합된 특징 세트로 학습된 16개의 머신러닝 분류기를 훈련시켰다.
• 주요 결과
• 혈장 cfRNA 경관: 엄격한 오염 필터링 후, 혈장의 미생물 알파 다양성은 조현병 환자와 건강 대조군 사이에 유의미한 차이를 보이지 않았다. 그러나 베타 다양성은 완만하지만 통계적으로 유의미한 변이를 보였다 (PERMANOVA P = 0.0316, R² = 0.0141). 탐색적 임계값(q < 0.10)에서 Flavonifractor plautii 및 Aeribacillus pallidus를 포함한 12종의 박테리아 종이 조현병 혈장에서 차별적으로 풍부하게 나타났다.
• 분변 미생물군: 혈장과 대조적으로, 분변 군집은 광범위한 재구조화를 보였다. 알파 다양성은 변화가 없었으나, 베타 다양성은 매우 유의미한 그룹 분리를 나타냈다 (PERMANOVA P = 0.0001). 분변에서는 81종의 종이 차별적으로 풍부하게 나타났다 (조현병에서 52종 증가, 29종 감소).
• 구획 간 일치성: 86개의 공유된 종 중에서 3종(Blautia obeum, Streptococcus salivarius, Evtepia gabavorous)이 두 구획 모두에서 차별적인 풍부도를 보였다. 특히, B. obeum과 S. salivarius는 조현병에서 분변에서는 감소했으나 혈장에서는 상승했는데, 이는 단순한 수동적 유출이 아닌 구획 특이적 역학을 시사한다. 그룹 수준의 혈장-분변 상관관계는 대조군보다 조현병 환자에서 유의미하게 강했으나, 개별 수준의 쌍 상관관계는 우연한 기대치를 초과하지 못했다. 이는 개별적인 전위(translocation) 사건보다는 그룹 수준의 조절이 강화되었음을 나타낸다.
• 기능적 연관성: 분류학적 연관성은 PANSS 점수와 약하게 연결되었으나, 혈장 cfRNA 유래 효소 기능은 증상 심각도와 더 강력한 연관성을 보였다. 구체적으로, 25개의 고유한 EC 반응이 PANSS 하위 척도와 연관되었다 (q < 0.10). 음의 계수가 반복되는 패턴은 특정 미생물 효소(예: 글루탐산 탈탄산효소 [EC 4.1.1.15])의 풍부도 감소가 증상 심각도(양성, 일반, 전체)의 증가와 상관관계가 있음을 시사한다.
• 진단적 분류: 분변 전용 모델이 혈장 전용 모델보다 높은 성능을 보였다 (최적 AUC: 0.964 vs. 0.873). 그러나 혈장 cfRNA 특징을 분변 메타게놈 특징과 통합했을 때 모든 16개 분류기에서 진단 성능이 소폭 개선되었다. 최적의 결합 모델(Ridge 회귀)은 0.976의 내부 탐색적 AUC를 달-성하였다. 특징 중요도 분석 결과, 분변 종이 주요 특징을 차지했으나, 두 가지 혈장 종(Aeribacillus pallidus 및 Flavonifractor plautii)이 결합 모델의 성능에 지속적으로 기여함을 확인하였다.

의의 및 주장
저자들은 본 연구가 조현병에서 최초로 오염을 고려한, 종 수준의 해상도를 가진 혈장 박테리아 cfRNA 경관을 제공한다고 주장한다. 주요 의의는 혈장 cfRNA가 단순히 분변 저장소를 반영하는 것이 아니라, 비중복적이고 보완적인 미생물 신호를 포착한다는 점을 입증했다는 데 있다.

주요 주장은 다음과 같다:

1. 전신적 재편: 본 연구 결과는 조현병에서 미생물 생태계가 전신적으로 재편된다는 모델을 뒷받로하며, 이는 개별적으로 검출 가능한 빈번한 전위 사건보다는 (잠재적으로 장벽 무결성 변화에 의해 유도된) 그룹 수준의 변화에 의해 발생한다.
2. 기능적 관련성: 혈장 cfRNA의 기능적 주석(특히 효소 수준)은 분류학적 조성보다 증상 심각도와 더 직접적으로 연결되어 있어, 향후 검증을 위한 잠재적 기능적 서명(signature)을 제공한다.
3. 진단적 유용성: 분변 메타게놈이 여전히 더 강력한 예측 인자이지만, 혈장 cfRNA 특징을 통합함으로써 진단적 식별력을 소폭 개선하였으며, 이는 정신 질환에 대한 '액체 생검' 접근법의 타당성을 입증한다.

저자들은 단면 설계로 인한 인과관계 추론의 한계, 식단 및 약물로 인한 잔류 혼란 변수의 가능성, 그리고 관찰된 AUC 값과 미생물 서명을 확인하기 위한 독립적인 외부 검증의 필요성을 명시하며 연구의 한계를 인정하였다. 또한, 보고된 분류 성능은 내부적이며 탐색적임을 강조하였다.