Rare variants drive high variance in human ancestral fitness at mutation-selection-drift balance

문제
인간 유전학의 핵심 질문은 무조건적으로 해로운 돌연변이(unconditionally deleterious mutations)의 유전적 부하(genetic load) 변이가 인간의 질병과 적응도(fitness)의 가변성에 실질적으로 기여하는가 하는 점이다. 흔한 변이(common variants)에 기반한 다유전자 위험 점수(polygenic risk scores, PRS)는 표준적인 방법이지만, 가계 연구에서 관찰되는 "누락된 유전력(missing heritability)"을 포착하는 데는 종종 실패한다. 기존의 돌연변이 부하 모델은 대개 "지체 부하(lag load)"(평균 집단 적응도와 가상의 돌연변이 없는 개체 간의 차이)에 초점을 맞추었으며, 이는 터무니없이 높은 적응도 비용을 초래하는 "돌연변이 부하 역설(mutation load paradox)"로 이어졌다. 또한, 서열 제약(sequence constraint)으로부터 도출된 표준적인 유해 돌연변이율($U_{obs}$) 추정치는 중립적인 비율로 고정되는 약한 해로운 돌연변이들을 반영하지 못하기 때문에 실제 돌연변이율($U$)을 과소평가할 가능성이 높다. 게다가 많은 모델은 유전적 부하의 분산, 특히 사이트가 단형성(monomorphic)일 확률을 무시하는데, 이는 새로운 돌연변이가 유효 집단 크기에 비해 드문 인간 게놈에서 매우 중요하다.

방법론
저자들은 단순한 돌연변이-선택 균형(mutation-selection balance)이 아닌, 돌연변이-선택-부동 균형(mutation-selection-drift balance, MSD)을 해결하기 위한 수학적 프레임워크를 개발하였다.

1. 수정된 $U$ 추정: 저자들은 관찰된 유해 돌연변이율($U_{obs}$)을 교정하여 전체 인간 유해 돌연변이율($U$)을 추정한다. 이들은 $U_{obs}$가 고정되는 약한 해로운 돌연변이를 놓친다는 점을 고려한다. 인간 비동형(non-synonymous) 돌연변이(특히 ESP 샘플) 및 초보존 비부호화 영역(ultra-conserved non-coding regions)으로부터 추론된 적응도 효과 분포(distribution of fitness effects, DFE)를 사용하여, 선택이 고정을 방지할 확률을 계산한다. 이를 통해 DFE 파라미터에 따라 $U \approx 1.34 U_{obs}$에서 $1.68 U_{obs}$에 이르는 교정 계수를 산출한다.
2. 분산 모델링: 저자들은 게놈 전체의 로그 상대 적응도($\sigma_r^2$)를 모델링한다. 사이트가 항상 분리(segregating)되어 있다고 가정하는 기존의 확산 근사법(diffusion approximations)과 달리, 이들은 대기 이론(queueing theory)(돌연변이 도착에 대한 포아송 과정 및 평균 분리 시간)을 사용하여 단형성 확률을 포함시킨다.
3. 파라미터: 모델은 현실적인 인간 1인당 디플로이드 사이트 돌연변이율($\mu = 1.515 \times 10^{-8}$), 유효 집단 크기($N_e = 11,300$), 그리고 다양한 평균 선택 계수($\bar{s}$, $0.00064$에서$0.01$ 사이)와 변동 계수(CV)가 2.43인 감마 분포 형태의 DFE를 활용한다.
4. 지표: 지체 부하 대신, 저자들은 75백분위수 개인과 25백분위수 개인 간의 적응도 이점(Q75-Q25 advantage)과 생물학적 타당성을 보장하기 위한 최적의 개체의 절대 적응도($W_b$)를 계산한다.

주요 기여

• 돌연변이율의 교정: 본 연구는 제약 기반 추정치보다 전체 인간 유해 돌연변이율이 유의미하게 높음을 입증하며, 약한 해로운 돌연변이의 포함으로 인해 $U > 3.8$ (상한선은 13에 근접)로 추정한다.
• 부하 분산의 정량화: 저자들은 MSD 균형 하에서 무조건적으로 해로운 돌연변이에 기인하는 조상 인간 적응도의 분산을 엄밀하게 계산한다.
• 희귀 변이의 역할: 모델은 대립유전자 빈도별로 부하의 분산을 구분하며, 희귀 변이(빈도 < 1%)와 초희귀 변이(빈도 < 0.01%)가 적응도 분산의 주요 동인임을 보여준다.
• 부동(Drift)의 영향: 논문은 부동이 해로운 대립유전자의 빈도를 높일 수 있지만, 동시에 대립유전자 소실(단형성)을 증가시켜 결정론적 모델에 비해 전체적인 적응도 분산을 감소시킨다는 점을 명확히 한다.

결과

• 적응도 분산: $U$의 불확실성을 고려할 때 두 전형적인 인간은 조상 적응도에서 17–33% 차이가 나며, 낮은 하한값($\bar{s}$)을 포함한 광범위한 대안적 DFE를 고려할 경우 **6–49%**까지 차이가 난다.
• 분산의 원천: 평균 선택 계수에 따라 희귀 변이(빈도 1% 미만)가 로그 적응도 분산의 약 **83–99%**를 설명한다. 초희귀 변이(<0.01%) 또한 상당한 부분을 기여한다.
• 생물학적 타당성: 모델은 집단 내 최적의 개체가 기대하는 생존 자손 수($W_b$)가 2.2에서 3.0 미만임을 예측하며, 이는 인간의 생애사적 제약(상한선 < 6.5)과 일치하여 돌연변이 부하 역설을 해결한다.
• 교차 종 적용 가능성: 높은 평균 선택 계수와 높은 돌연변이율을 가진 다른 포유류(예: 생쥐)에 대해서는 높은 부하 변이가 예측되지만, 평균 선택 계수가 훨씬 작은 종(예: 초파리)에 대해서는 그렇지 않다.
• 부동의 영향: 일부 기존 주장과 달리, 부동(및 단형성 확률)을 포함하면 결정론적 돌연변이-선택 균형 모델에 비해 부하의 변동이 완만하게 감소한다.

의의
본 논문은 희귀하고 무조건적으로 해로운 돌연변이가 인간의 적응도와 질병 감수성의 주요 원인이라고 주장한다. 이러한 높은 부하 분산은 왜 가계 연구에서 관찰되는 유전력이 흔한 변이에 의존하는 전장 유전체 연관 분석(GWAS)에서 자주 누락되는지를 설명하는 데 도움이 된다. 저자들은 질병 위험이 높은 개인을 식별하기 위해서는 흔한 변이로부터 유도된 다유전자 위험 점수뿐만 아니라, 전장 유전체 시퀀싱을 통한 변이 효과 예측과 같이 희귀 및 초희귀 변이를 점수화할 수 있는 방법이 필요하다고 제언한다. 이 결과는 주변부 유전자의 돌연변이가 질병 위험에 확산적으로 기여한다는 "옴니제닉(omnigenic)" 모델을 뒷받침한다. 또한 저자들은 유전적 부하 개념이 차별적인 정책을 위한 역사적 오용에 사용되는 것을 명시적으로 비난하며, 공평한 사회를 증진하기 위해 유전적 부하를 이해하는 것의 윤리적 중요성을 강조한다.