Seven replicated genomic associations of myalgic encephalomyelitis/chronic fatigue syndrome: a biobank study

문제 정의
근통성 뇌척수염/만성 피로 증후군(ME/CFS)은 여성 편향적인 경향을 보이는 쇠약성 질환이지만, 진단용 바이오마커, 효과적인 치료법, 또는 명확히 이해된 병인론이 결여되어 있다. 기존의 유전학 연구는 소규모의 이질적인 코호트에 의존하고 단일 진단 코드에 의존해 왔기 때문에 재현 가능한 결과를 도출하는 데 어려움을 겪어 왔으며, 이로 인해 재현 가능한 발견이 부족했다. 이전에는 독립적인 코호트 전반에 걸쳐 재현되는 ME/CFS의 게놈 위험 로커스(loci)가 없었다. 또한, 사회경제적 지위(SES)와 유전자-성별 상호작용이 위험을 수정하는 역할에 대해서는 선행 연구들의 상충하는 보고와 방법론적 한계로 인해 논쟁의 여지가 있다.

방법론
저자들은 UK Biobank(UKB)와 All of Us(AoU) 연구 프로그램을 활용하여 전전체 연관 분석(GWAS)을 수행하였다.

• 표현형 전략: 단일 진단 코드에 의존하는 대신, 사례(cases)와 대조군(controls)을 정의하기 위해 다중 증거 표현형 전략을 채택하였다. 사례군은 자가 보고된 임상 진단, 불량/보통 수준의 전반적 건강 상태 등급, 그리고 ME/CFS에 관한 특정 긍정 응답을 모두 갖추어야 했다. 대조군은 양호/우수 수준의 건강 등급을 가져야 하며, 의료 기록이나 설문 조사에서 ME/CFS의 증거가 없어야 했다.
• 코호트: 세 개의 서로 다른 모집단이 설정되었다:
  1. 발견 코호트 (UKB1): 1,268 사례 및 113,132 대조군.
  2. 재현 코호트 1 (UKB2): 319 사례 및 28,471 대조군.
  3. 재현 코호트 2 (AoU): 371 사례 및 33,140 대조군.
• 통계적 프레임워크: 본 연구는 표적 학습(Targeted Learning)에 기반한 준매개변수적, 이중 강건 추정 워크플로우인 TarGene(Targeted Genomic Estimation)를 활용하였다. 이 접근 방식은 기존의 매개변수적 GWAS 방법이 가진 한계를 해결하면서, 모델 오지정(misspecification)과 집단 층화로부터 발생하는 편향을 최소 जबकि 통계적 검정력을 극대화하기 위해 선택되었다.
• 분석 파이프라인:
  • 발견 단계: UKB1 내 670,243개의 유전형 SNP에 대해 GWAS를 수행하였으며, 허위 발견율(FDR)을 5%로 제어하였다.
  • 재현 단계: 발견 단계에서 식별된 176개의 FDR 유의적 변이들을 두 개의 재현 코호트에서 테스트하였다.
  • 미세 매핑(Fine-mapping) 및 공국소화(Colocalisation): 재현된 로커스에 대해 SuSiE를 사용하여 미세 매핑을 수행하여 인과적 변이의 신뢰 집합(credible sets)을 생성하였다. 이후 이를 GTEx v10 데이터베이스의 발현 양적 형질 위치(eQTL)와 공국소화하였다.
  • 상호작용 분석: 유전자-성별 및 유전자-결핍(Townsend 결핍 지수 및 다중 결핍 지수 사용) 상호작용을 TarGene을 사용하여 테스트하였다.
  • 검증: 엄격한 임상 기준(캐나다 컨센서스 및 IoM 2015)을 적용하여 운동 후 악화(PEM)를 요구하는 연구인 DecodeME 코호트(2,143 사례)를 대상으로 룩업(look-up)을 수행하였다.

주요 결과

• 발견: 분석 결과 UKB1 발견 코호트에서 ME/CFS 위험과 유의하게 연관된 176개의 변이가 식별되었다 (FDR < 5%).
• 재현: 7개의 게놈 위험 로커스가 두 개의 서로 다른 재현 코호트(UKB2 또는 AoU) 중 적어도 하나에서 재현되었다. 그러나 단 하나의 변이도 다중 검정 보정 후 세 코호트 모두에서 유의성을 확보하지 못했다.
  • 재현된 7개의 변이는 다음과 같다: rs115186419 (ETV5 근처), rs73175505 (CSMD1 근처), rs261902 (BICD1 내), rs117553493 (CLYBL 근처), rs72741654 (RORA 근처), rs74963073 (GRIN2A 내), rs76847656 (SBK1 근처).
• 미세 매핑 및 공국소화:
  • 염색체 13 로커스에서, 미세 매핑은 피각(putamen, 기저핵)에서의 CLYBL 발현 감소와 공국소화되는 신뢰 집합을 해결하였다. 위험 대립유전자는 낮은 CLYBL 발현과 연관되어 있다 (사후 확률 H4 = 0.994).
  • 주목할 점은, 알려진 기능 상실(loss-of-function) 변이인 CLYBL Arg259 정지-획득(stop-gain) 변이는 ME/CFS 위험과 연관되지 않았다. 이는 해당 신호가 단백질 절단보다는 발현 용량을 조절하는 조절 변이에 의해 주도됨을 시사한다.
• 상호작용: 유전자-성별 또는 유전자-결핍 상호작용 중 다중 검정 보정을 통과한 것은 없었다. 그러나 rs261902 및 rs76847656에서 성별 차등 효과, 그리고 rs73175505에서 결핍 차등 효과에 대한 명목적(nominal) 신호가 관찰되었다.
• DecodeME 검증: 재현된 4개의 테스트 가능한 변이 중 어느 것도 DecodeME 코호트에서 연관성을 보이지 않았다. 저자들은 이러한 불일치를 표현형 정의(DecodeME는 PEM을 요구함), 연령 분포, 그리고 모집 전략(지역사회 기반 vs 바이오뱅크)의 차이 때문으로 돌렸다.

의의 및 주장
본 논문은 독립적인 코호트에서 재현되는 유전적 연관성을 보고한 최초의 ME/CFS 연구라고 주장하며, 7개의 후보 로커스를 식별하였다. 저자들은 다중 증거 표현형과 TarGene 통계 프레임워크의 사용이 기존의 단일 코드 연구들이 놓쳤던 신호들을 탐지할 수 있게 해주었다고 강조한다.

연구는 비록 모든 변이가 세 코호트 모두에서 재현되지는 않았으나(이는 표현형의 이질성, 인구 집단 차이, 그리고 작은 규모의 재현 코호트에서의 "승자의 저주" 때문으로 보임), 7개의 재현된 로커스가 향후 생물학적 조사를 위한 후보 타겟을 제공한다고 겸허하게 결론짓는다. 염색체 13 신호와 피각에서의 CLYBL 발현 사이의 공국소화는 미토콘드리아 비타민 B12 대사와 ME/CFS 병태생리 사이의 잠재적 연결을 강조한다. 저자들은 이 결과가 인과적 증거를 구성하는 것이 아니라, ME/CFS의 생물학적 메커니즘에 대한 후속 연구를 위한 토대를 제공한다고 명시적으로 밝혔다. 또한, DecodeME와의 불일치는 질병이 어떻게 정의되느냐(예: PEM의 포함 여부)에 따라 유전적 구조가 달라질 수 있음을 시사한다고 주의를 주었다.