Population-scale genomics reveals divergent pathogenicity of variant classes across paralogous collagen IV genes

문제 정의
알포트 증후군(Alport Syndrome, AS)은 제4형 콜라겐의 $\alpha3(IV)$, $\alpha4(IV)$, $\alpha5(IV)$ 사슬을 인코딩하는 상동 유전자인 COL4A3, COL4A4, COL4A5의 병원성 변이에 의해 발생하는 유전성 신장 질환이다. 이형 접합(biallelic) 또는 X-연관(X-linked) 변이의 임상적 중증도는 잘 규명되어 있는 반면, COL4A3 및 COL4A4의 이형 접합 자동상동(heterozygous autosomal) 변이에 대한 병원성은 덜 알려져 있다. 현재의 임상 가이드라인과 변이 분류 체계는 동일한 변이 클래스(예: 글리신 치환, 절단 변이, 비콜라겐 도메인[NC1] 미스센스 변이)에 대해 동등한 병원성을 가정하며, 이 두 상동 유전자를 기능적으로 상호 교환 가능한 것으로 취급하는 경향이 있다. 그러나 이러한 가정은 인구 규모 수준에서 엄격하게 검증된 적이 없으며, 이는 잠재적으로 부정확한 위험 계층화 및 유전 상담으로 이어질 수 있다.

연구 방법론
저자들은 영국 바이오뱅크(UK Biobank; UKB; n ≈ 500,000)의 전전체 게놈 시퀀싱(WGS) 데이터를 활용하였고, All of Us 연구 프로그램(n ≈ 414,000)에서 그 결과를 재현하였으며, 분석 대상을 유럽 혈통의 개인으로 제한하였다. 연구에는 다각적인 접근 방식이 사용되었다:

1. 표현형 정의: 사구체 혈뇨, 단백뇨(ICD-10 코드 및 요중 알부민-크레아티닌 비율 사용), 그리고 중증 만성 신장 질환(CKD)에 대한 사례를 정의하였으며, 비사구체 원인을 엄격히 배제하였다.
2. 변이별 연관성: 개별 희귀 코딩 변이와 신장 표현형 사이의 연관성을 평가하기 위해 피셔 정확 검정(Fisher's Exact Test)을 사용하였다.
3. 유전자 기반 결합(Collapsing) 및 PheWAS: 희귀 예측 병원성 변이를 세 가지 기능적 범주, 즉 글리신 치환(Gly), 절단 변이(Trunc), NC1 도메인 미스센스 변이로 분류하였다. 이들은 COL4A3와 COL4A4에 대해 각각 별도로 이진(binary) 보유 상태 변수로 결합되었다. 1,363개의 ICD-10 표현형에 대한 연관성을 테스트하기 위해 SAIGE를 사용하여 표현형-전체 연관 연구(PheWAS)를 수행하였다.
4. 효과 크기 비교: 각 변이 클래스에 대해 COL4A3와 COL4A4 사이의 유전적 효과 크기($\beta$) 차이를 두 꼬리 Z-검정을 사용하여 공식적으로 테스트하였다.
5. 전전체 연관 연구(GWAS): 혈뇨와 관련된 로커스(loci)를 식별하기 위해 GWAS를 수행하였으며, 신호가 희귀 변이에 의한 것인지 아니면 독립적인 공통 변이/조절 인자에 의한 것인지 결정하기 위해 희귀 COL4A3/A4 변이의 가중 부하 점수(weighted burden score, WBS)를 조정한 이차 분석을 실시하였다.

주요 결과

• 변이 연관성: 본 연구는 혈뇨 또는 단백뇨와 유의미하게 연관된 64개의 COL4A3 및 92개의 COL4A4 희귀 변이를 식별하였다. 글리신 치환은 모든 표현형에서 가장 빈번한 클래스였다.
• 유전자 특이적 발현 차이: 변이 클래스에 따른 두 유전자 간의 병원성 발현 차이가 핵심적인 발견이었다:
  • 글리신 치환: COL4A3와 COL4A4 모두에서 유사한 위험을 부여하였다.
  • 절단 및 NC1 미스센스 변이: 이 변이들은 COL4A4 보유자에서 혈뇨 및 단백뇨와 강력하게 연관되었다. 반면, COL4A3 보유자의 경우, 유사한 보유 빈도와 예측된 병원성 점수(CADD/REVEL)에도 불구하고 절단 또는 NC1 변이에 의한 연관성이 실질적으로 약화되거나 나타나지 않았다.
  • 이러한 차이는 영국 창시자 변이인 p.Ser969Ter를 제외한 후에도 지속되었으며, All of Us 코호트에서도 재현되었다.
• 보유율: COL4A4 절단 및 NC1 미스센스 보유자는 COL4A3 보유자(단백뇨의 경우 각각 17.1% 및 12.1%)에 비해 단백뇨 사례율(각각 52.7% 및 37.2%)이 유의하게 높았다.
• 유전적 구조: COL4A3/A4 로커스는 혈뇨의 지배적인 유전적 결정 요인으로 식별되었다. 연관 신호는 전적으로 희귀 코딩 변이의 집합적 효과에 기인하였다. 가중 부하 점수를 조정했을 때 로커스 신호가 사라졌으며, 독립적인 공통 변이 또는 트랜스-작용(trans-acting) 조절 인자의 증거는 없었다. HLA-B 및 THSD4 로커스는 혈뇨와 연관된 독립적인 신호로 식별되었다.

의의 및 주장
본 논문은 COL4A3와 COL4A4가 절단 및 NC1 도메인 변이의 병원성에 있어 기능적으로 상호 교환 가능하지 않다는 것을 보여주는 첫 번째 인구 규모의 증거를 제공한다고 주장한다. 저자들은 다음과 같이 결론짓는다:

1. 유전자 정체성의 중요성: 콜라겐 IV 신병증에 대한 위험 계층화 및 변이 해석은 단순히 변이 클래스가 아닌 유전자의 정체성을 고려해야 한다. COL4A4 절단 및 NC1 변이는 COL4A3의 대응 변이보다 실질적으로 더 강한 병원성 효과를 부여한다.
2. 임상적 함의: 이러한 발견은 두 상동 유전자가 동일한 병원성을 가질 것이라는 가정에 도전하며, 현재의 가이드라인이 COL4A3 절단/NC1 보유자의 위험은 과대평가하고 COL4A4 보유자의 위험은 과소평가할 수 있음을 시사한다.
3. 광범위한 원칙: 본 연구는 동일한 헤테로머 복합체를 인코딩하는 상동 유전자들이 기능 상실 변이에 대한 내성에서 상당한 차이를 보일 수 있음을 입증한다. 이는 임상 유전학에서 '변이 클래스 인식'이 아닌 '유전자 인식' 프레임워크가 필요함을 시사하며, 이 원칙은 다른 헤테로다이머(heteromultimer) 인코딩 유전자 가족에도 적용될 수 있다.
4. 기전적 통찰: 결과는 $\alpha4(IV)$ 사슬이 $\alpha3(IV)$에 비해 사슬 선택 또는 프로토머 안정화에서 불균형적으로 중요한 역할을 할 수 있음을 시사하지만, 구체적인 생화학적 기전은 여전히 규명되어야 할 과제로 남아 있다.

저자들은 유럽 혈통으로의 제한과 무증상 미세 혈뇨의 유병률을 과소평가할 수 있는 병원 기반 ICD-10 코드에 대한 의존성을 포함한 한계점을 언급하였다. 따라서 보고된 보유자 사례율은 침투율(penetrance)에 대한 보수적인 추정치로 해석된다.