장내미생물 ‘하위종(OSU)’로 본 대장암 예측의 새 지평

🧭 왜 ‘하위종’인가? — 종 수준의 한계를 넘어

기존 메타지놈 도구는 주로 ‘종’ 수준에서 풍부도를 계산합니다. 그러나 같은 종 내에서도 유전자 구성이 크게 달라 기능이 달라질 수 있어, 종 단위만 보면 질환 연관성이나 기전을 놓치기 쉽습니다. 하위종(형제 하위종 간 기능 차이가 뚜렷한 가장 좁은 공식 분류 단위)을 쓰면 개인·코호트 간 이질성은 줄이면서, 기능 차이는 살릴 수 있습니다.

---

🏗️ HuMSub는 어떻게 만들었나 — OSU와 panhashome의 아이디어

연구진은 인간 장내 세균 유전체(품질 필터링 후)를 바탕으로 코딩서열을 스케치(sketch)하고, 예측 강도(predictive strength)로 군집화하여 운영 하위종 단위(OSU, operational subspecies unit)를 정의했습니다. 그 결과 225,918개 유전체(3,483 종)에서 977 종에 속한 5,361개 OSU를 얻었고, 종의 28%에서 무시되던 하위종 변이가 드러났습니다.

 

또한 각 하위종에 고유한 해시값 모음(panhashome)을 만들어 직접 정량이 가능하도록 했습니다. 모의 데이터에서 F1-score 중앙값 0.794, L2 거리 0.144로 높은 성능을 보였고, 종 수준 대표 도구(MetaPhlAn4) 대비 계산 자원도 줄었습니다.

 

HuMSub 카탈로그의 하위종은 세균역 전반에 광범위하게 존재하며 직접 정량 가능함

 

---

🌍 전 세계에서 통하는가? — 생지리 분포와 일반화 가능성

국가 정보가 있는 5,272개 공공 메타지놈 샘플을 검색하니, 정량된 OSU의 62%가 4개 대륙 중 최소 3곳에서 발견되었고, 0% 유병인 OSU는 2.5%에 불과했습니다. 지리적 특이성(Geographical Enrichment Score, GES)을 정의해 GES=0.4를 임계값으로 삼았을 때, 1,869 OSU는 광범위 공유, 338 OSU는 지역 특이로 분류되었습니다. 이는 종·균주 수준 분석에서 놓치는 정보입니다.

 

하위종의 전 지구적 생지리 분포

 

---

🧪 대장암 메타분석 — 하위종은 ‘형제 중 누구’가 문제인지 보여준다

7개 연구(대장암 555, 대조 530)에서 하위종 정량을 수행하니, 2,800 OSU 중 218개가 유의하게 연관되었고, 그중 104개는 ‘어떤 하위종만’ 연관(형제는 비연관)인 경우였습니다. 예컨대 Fusobacterium animalis는 두 하위종 중 OSU 001002만 대장암에서 증가했습니다. 반대로 부모 종이 비연관인데 특정 하위종만 연관된 28예도 관찰되어, 종 수준에선 보이지 않는 신호를 하위종이 드러냄을 확인했습니다.

 

하위종에는 종 수준에서 보이지 않는 고유 정보가 내재함

 

---

🤖 머신러닝 성능 — 임상적 예측력의 실질적 향상

같은 구조의 모델로 종 vs 하위종 입력을 비교했습니다. LODO(Leave-One-Dataset-Out) 설정에서 하위종 기반 AUROC 중앙값 0.838(일부 데이터셋 0.89+)이 종 기반 0.785를 모든 6개 테스트 중 5개에서 앞섰습니다. 또한 대변잠혈(FOBT, fecal occult blood test) 결과를 가중치로 포함하자 AUROC 0.893까지 상승했습니다. 특징 중요도 분석을 활용하면 상위 64개 OSU만으로도 성능 고원(plateau)에 도달해, 해석 가능하고 임상에 적합한 시그니처 구성이 가능했습니다.

 

하위종은 대장암(CRC) 예측에서 더 우수한 지표임

 

---

🔬 무엇이 달라지는가 — 비타민 B12·아르기닌·헴 대사로 본 기전 단서

하위종 간 유전자 차이를 들여다보니, 높은 영향력의 변이(∣ΔΔG∣>1)가 OSU 특이 보조유전자보다 더 흔했고, 핵산·아미노산 생합성 경로가 특히 영향 받았습니다. 예를 들어,

• Ruthenibacterium lactatiformans의 CRC-연관 OSU 001003과 비연관 001001/001002를 비교하면, B12(코발라민) 생합성 관련 유전자들에 파괴적 변이가 비연관 쪽에 집중되어 B12 생성 능력 상실이 시사되었습니다.
• Enterocloster 한 하위종에선 아르기닌 생합성(argH)의 불안정화 변이가 비연관 OSU에 있어 대사 수준 차이가 추정됩니다.
• Porphyromonas에서는 헴 생합성(hemH, gltX) 유전자 불안정화가 비연관 OSU에 위치했습니다. 헴은 상피 과증식·염증과 연결되어 대장암에 유리한 환경을 만들 수 있습니다.

이처럼 “같은 종 안의 누구(하위종)가 무엇을(기능 유전자) 어떻게(고영향 변이·유전자 보유/결손) 바꾸는가”가 보이기 시작합니다.

 

형제 하위종은 숙주 생리에 서로 다르게 영향을 미침

 

---

⚠️ 한계와 주의

조립 단편(Contig) 가장자리의 유사유전자(pseudogene) 오분류 가능성, 품질점수의 정량적 표준 부재 등 방법론적 한계가 있습니다. 다만 실제 데이터 기준으로 가장자리 유전자 비율을 감안하면 영향은 매우 제한적으로 추정됩니다.

---

✍️ 한줄평

“하위종 해상도는 ‘누가’ 문제인지와 ‘무엇을’ 바꾸는지를 한 번에 보여주는, 마이크로바이옴 예측·기전 연구의 실용적 분기점입니다.”

 

참고문헌 : DOI: 10.1016/j.chom.2025.07.015