Urinary metabolomic signatures of muscle injury and recovery in elite football players

기술 요약: 엘리트 축구 선수의 근육 부상 및 회복에 나타나는 소변 대사체학적 특징

문제 제기
근육 부상은 엘리트 축구에서 시간 손실(time-loss)과 경기력 저하를 야기하는 주요 원인이다. 현재의 선수 모니터링 전략은 전자 퍼포먼스 트래킹 시스템(EPTS)에서 도출된 외부 부하 지표(예: 달린 거리, 속도)와 주관적인 내부 부하 지표(예: 인지된 노력도, 피로도)에 크게 의존하고 있다. 그러나 이러한 방법들은 주어진 부하에 대한 생물학적 반응이 매우 개별적이며 복잡한 대사 및 염증 네트워크의 영향을 받기 때문에, 개별적인 생리학적 스트레스와 생물학적 비용을 완전히 포착하지 못하는 경우가 많다. 따라서 구조적 영상 진단 및 기존 모니터링을 보완할 수 있는, 근육 부상 및 조직 치유와 관련된 대사 변화를 반영하는 객관적이고 비침습적인 바이오마커가 필요하다.

연구 방법론
본 관찰 종단 연구는 2두 시즌 동안 121명의 엘리트 남성 축구 선수로부터 수집된 287개의 첫 기상 소변 샘플을 분석하였다. 데이터셋은 임상적으로 두 그룹으로 분류되었다: 10명의 서로 다른 선수로부터 수집된, 임상적으로 진단된 근육 부상 기간의 샘플 30개와 부상이 없는 대조 기간의 샘플 257개이다. 근육 부상은 UEFA 합의 권고안(시간 손실 정의)에 따라 정의되었으며, 임상 평가 및 영상(MRI/초음파)을 통해 확인되었다.

• 분석 접근 방식: 표적형 초고성능 액체 크로마토그래피-질량 분석법(UPLC-MS)을 사용하여 특정 대사체 클래스(AccQ-Tag 유도체를 통한 아미노산 및 트립토판 관련 대사체)를 정량화하였다.
• 통계 분석:
  • 단변량 분석: 부상 상태와 대조 상태 사이에서 유의미한 차이를 보이는 개별 대사체를 식별하기 위해 Welch의 t-검정, 폴드 변화(fold-change) 평가 및 Hedge's g 효과 크기 추정을 사용하였다.
  • 경로 분석: 대사 경로의 농축을 확인하기 위해 KEGG 데이터베이스와 mummichog 알고리즘을 이용한 과잉 표현 분석(Over-representation analysis, ORA)을 수행하였다.
  • 다변량 분석: 조정된 대사 패턴을 식별하기 위해 부분 최소 제곱 판별 분석(PLS-DA)을 활용하였다. 부상 샘플 수가 적었기 때문에(n=30), 학습 데이터를 극대화하기 위해 Leave-One-Out 교차 검증(LOOCV)을 적용하였다. 모델의 유의성을 검증하기 위해 순열 검정(300회 반복)을 실시하였다. 또한, 개별 간 변동성에 대한 모델의 견고성을 평가하기 위해 Leave-One-Player-Out(LOPO) 교차 검증을 이용한 민감도 분석을 수행하였다.
  • 특징 선택: 변수 중요도(Variable Importance in Projection, VIP) 점수에 기반한 후방 단계적 제거법(Backward stepwise elimination)을 사용하여 식별력을 가진 감소된 대사체 하위 집합을 식별하였다.
  • 종단적 추적: 부상당한 선수들의 대사 점수를 시계열로 추적하여 재활 과정 중의 정상화 경향을 관찰하였다.

주요 결과

• 대사체 서명(Metabolic Signatures): 근육 부상은 네 가지 주요 대사 영역에서의 조정된 변화와 관련이 있었다:
  1. 아미노산 회전(Amino Acid Turnover): 알라닌, β-알라닌, 프롤린, 오르니틴, 라이신의 유의미한 변화.
  2. 에너지 스트레스: 퓨린 대사의 변화, 특히 세포 에너지 고갈을 나타내는 하이포잔틴(hypoxanthine)의 상승.
  3. 트립토판 이화작용: 크이뉴레닌 경로 대사체(xanthurenic acid, kynurenic acid, tryptamine 포함)의 유의미한 상승.
  4. 메틸화: S-아데노실메티오닌(S-adenosylmethionine)의 변화.
• 경로 농축: 경로 분석 결과 "아르기닌 및 프롤린 대사"가 가장 유의미하게 농축된 경로로 식별되었으며, 이는 대사 변화가 조직 복구 및 콜라겐 합성 과정으로 수렴됨을 시사한다.
• 분류 성능: PLS-DA 모델은 통계적으로 유의미한 분류 성능을 입증하였다 (순열 검정 기준 p < 0.003).
  • 교차 검증된 AUROC: 0.88.
  • 민감도(Sensitivity): 73%; 특이도(Specificity): 88%.
  • 음성 예측도(NPV): 97%.
  • 참고 (모델 견고성): 샘플 수준의 LOOCV는 0.88의 AUROC를 나타냈으나, 더 엄격한 선수 수준의 LOPO-CV 결과 AUROC는 0.72로 감소하였다. 저자들은 이를 데이터셋의 불균형(10명의 선수로부터 얻은 30개의 부상 샘물)과 개별 간 대사 변동성을 처리하기 위한 샘플 수준의 보정 필요성 때문이라고 설명하였다.
• 주요 판별 대사체: 다변량 모델은 β-아미노이소부티르산(BAIBA), 하이포잔틴, 잔투레닉산(xanthurenic acid), β-알라닌, 알라닌을 가장 강력한 샘플 판별 기여자로 식별하였다.
• 종단적 역학: 부상당한 선수들에 대한 종단적 분석 결과, 재활 과정 중 대사 프로파일이 점진적으로 정상화되어, 선수들이 경기 복귀(RTP) 시점에 도달함에 따라 대사 점수가 기본 대조군 값에 근접하는 것이 관찰되었다. 이러한 경향은 흉터 성숙 및 부종 해소에 대한 구조적 MRI 증거와 함께 관찰되었다.

의의 및 주장
본 논문은 소변 대사체학이 엘리트 축구 선수의 근육 부상과 관련된 조정된 대사 변화를 식별할 수 있으며, 근육 부상 및 조직 손상과 회복에 대한 "내부 부하"와 생리학적 스트레스를 들여다볼 수 있는 비침습적인 창을 제공한다고 주장한다. 본 연구는 대사체 프로파일링이 외부 부하 모니터링 및 구조적 영상 진단을 보완하는 도구로서 유용하다는 점을 뒷받침한다.

구체적으로, 저자들은 다음과 같이 상정한다:

1. 보완적 모니터링: 대사체 서명은 외부 부하 지표나 거시적 구조 영상만으로는 완전히 반영되지 않을 수 있는 생리학적 스트레스(예: 에너지 고갈, 염증 조절)를 포착함으로써 부상에 대한 전신적 관점을 제공한다.
2. 회복 평가: 재활 중 대사 프로파일의 점진적 정상화는 소변 대사체학이 회복 역학을 모니터링하고 잠재적으로 경기 복귀 결정을 내리는 데 객적인 지표로 사용될 수 있음을 시사한다.
3. 중개 가능성(Translational Potential): 제한된 판별 대사체 패널(예: 하이포잔틴, 잔투레닉산, BAIBA)의 식별은 개인별 선수 모니터링을 위한 표적화된 일상적 임상 분석법 개발의 타당성을 뒷받침한다.

저자가 인정한 한계점
저자들은 다음과 같은 몇 가지 제약 사항을 명시하였다:

• 교란 요인: 영양 섭취, 급성 훈련 부하, 약물 복용 등이 기록되지 않았으며, 이는 소변 대사체에 영향을 미칠 수 있다.
• 표본 크기 및 일반화 가능성: 확인된 부상 사례가 제한적이었고(n=30), 단일 센터 코호트(10명의 선수)에 국한되어 있어 통계적 일반화 능력이 제한된다.
• 전신적 특성: 소변 대사체학은 골격근 부상과 반드시 직접적으로 관련되지 않은 과정을 포함할 수 있는 전신적 반응을 반영한다.
• 인과관계: 본 연구는 관찰 연구이므로 발견된 연관성에 대한 인과적 해석에는 한계가 있다.

본 논문은 이러한 결과가 유망하지만, 개인 맞춤형 스포츠 의학을 실현하기 위해서는 더 큰 규모의 독립적인 다기관 코호트에서의 검증과 구조적 영상 및 임상 평가와의 통합이 필요하다고 결론짓는다.